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1.0  INTRODUCTION 


The  Gas  Diagnostics  Division  of  The  University  of  Tennessee  Space 
Institute  has  constructed  and  delivered  to  the  U.  S.  Army  Ballistics  Research 
Laboratory  (BRL)  a  Particle  Sizing  Interferometer  (PSI)  designed  to  measure 
droplet  sizes  over  ranges  from  0.3  to  greater  than  1000  micrometres.  The 
PSI  is  perhaps  the  only  nonintrusive  diagnostic  instrument  which  can  obtain 
time  histories  of  particle  or  droplet  size  distributions  in  hostile  environ¬ 
ments.  While  the  PSI  is  simple  in  concept,  its  successful  utilization  and 
application  requires  a  good  understanding  of  many  subtle  points  pertinent  to 
the  particles  being  measured  and  their  effects  on  instrument  operation. 
Therefore,  the  purpose  of  this  report  is  to  1)  describe  the  expected  operation 
of  the  system  in  realistic  injector  sprays  (a  subject  for  study  at  BRL), 

2)  develop  in  greater  detail  theoretical  references  alluded  to  in  previous 
reports,  and  3)  provide  suggestions  for  extending  the  application  to  other 
potential  applications  such  as  burning  powders. 


1.1  Theoretical  Synopsis 

/ 

Consider  the  optical  system  shown  in  figure  *.]*  A  beam  from  a  laser  is 
split  into  two  equal  intensity  beams.  A  lens  system  makes  the  beams  cross 
and  focus  at  a  common  origin  called  the  geometric  center.  Define  the  included 
angle  between  the  beams  as  a.  Near  the  geometric  center  the  wavefronts  are 
planar.  A  Huygen’s  diagram  of  the  wavefronts  shows  that  planar  interference 
fringes  are  formed  which  are  perpendicular  to  the  plane  defined  by  the  beam 
centerlines  and  are  parallel  to  the  bisector  between  the  beams.  The  distance 
6  between  the  periodic  fringes  is  given  by 

6  =  A/2sin(a/2)  (1) 

where  X  is  the  wavelength  of  the  coherent  light.  At  this  point  it  is 
convenient  to  divide  the  particles  which  may  interact  with  the  interference 
fringes  into  two  classes.  Class  I  particles  are  "small”  (diameters  comparable 
to  an  optical  wavelength)  scattering  particles  with  scattering  characteristics 
which  can  be  described  by  Fraunhofer  diffraction  or  Lorenz-Mie  scattering 
theory.  The  light  scattered  by  these  particles  is  assumed  to  be  observed 
with  a  scattered  light  collection  lens  centered  on  the  bisector  between  the 
beam  in  the  forward  direction.  Class  II  particles  are  "large"  (diameters  much 
greater  than  a  wavelength)  refracting  or  reflecting  particles  with  scattering 
characteristics  which  can  be  described  using  geometric  optics.  Light 
scattered  from  these  particles  is  observed  at  some  angle  off  the  bisector 
between  the  beams  with  the  scattered  light  collection  lens  centered  in  a  plane 
perpendicular  to  the  plane  of  the  beams.  The  following  discussion  presents 
two  simple  models  for  understanding  the  signal  generated  by  the  two  classes 
of  particles. 


1.2  Class  I  Particles 

When  a  particle  (assumed  spherical)  much  less  than  5  in  diameter  crosses 
the  fringe  pattern,  it  can  be  assumed  to  be  uniformly  illuminated  across  its 


7 


PiOal  V0ILJM1 


LU 

PC  ui 

3  -> 


*-2£ 

w  t: 

o>- 

DC*^ 

SI  3 

ttO 


2> 


iuN 


8 


FIGURE  1.  PRINCIPLE  OF  THE  PARTICLE  SIZING  INTERFEROMETER 


diameter  at  all  points  along  its  path  through  the  fringe  pattern.  The 
light  which  is  scattered  by  the  particle  is  proportional  to  the  observable 
flux  illuminating  it.  Thus,  measurement  of  the  time  period,  X,  of  the 
scattered  light  is  related  to  the  velocity,  v,  of  the  particle  through  the 
relationship 


v  =  6/t  .  (2) 

As  the  size  of  the  scattering  particle  increases  relative  to  6,  the 
illumination  of  the  particle  is  no  longer  uniform  and  must  be  averaged 
over  the  cross-sectional  area  of  the  particle.  The  non-uniform  illumination 
of  the  particle  results  in  a  reduction  of  the  contrast  or  visibility  of  the 
scattered  light  signal.  Let  Imax  be  the  maximum  value  in  intensity  in  a 
period  of  the  scattered  light  from  a  particle  and  Imin  the  next  successive 
minimum.  The  visibility,  V,  can  then  be  defined  as 

I  -  I  . 

_  max  min 

v  "  i  +  x  *  ^ 

max  min 


It  is  straightforward  to  show  that  V  is  fully  equivalent  to  the  ratio  of 
AC  amplitude  divided  by  the  DC  amplitude  of  the  scattered  light  signal. 

The  high  frequency  "Doppler"  portion  of  the  signal  is  defined  as  the  "AC". 

It  usually  has  many  cycles  of  information  relative  to  that  of  the  DC 
component  (the  DC  component  refers  to  the  Gaussian  shaped  low  frequency  terra 
describing  the  signal).  Analytically,  the  visibility  may  be  written  as 


V  ~ 


/ 

A 

J£L 


IqCOs  (27Ty/6)dAp 


/  I  dA 

J  op 

A 

P 


(4) 


where  Ap  is  the  cross-sectional  area  of  the  particle,  IQ  is  intensity  distri¬ 
bution  across  one  of  the  illuminating  beams,  and  y  is  the  coordinate  normal 
to  the  fringe  planes.  When  IQ  is  a  Gaussian  function  (TEMQOq*  laser  beam  - 
the  subscripts  indicate  the  order  of  the  Hermite  polynomial  used  to  describe 
the  oscillation  mode  of  the  laser  cavity),  it  can  be  shown  that  Equation  4 
is  an  accurate  approximation  over  a  depth  of  field,  £,  given  by 

H=0.8b  /a  (5) 

o 

-2 

where  bQ  is  the  radius  of  the  e  intensity  point  in  the  illumination  beam. 

For  depths  of  field  greater  than  £,  the  signal  visibility  is  a  function  of 
particle  size  and  position  in  the  illumination.  In  order  to  simplify  Equation 
4  for  Gaussian  beams  and  still  maintain  accuracy,  it  is  required  that  the 
particle  diameter,  D,  satisfy  the  relationship 

D  <  0.2b  (6) 

—  o 
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and  for  6  to  satisfy 


6<0.2bQ  .  (7) 

Under  these  conditions,  V  for  a  sphere  can  be  written  as 

V  «  2J1(7TD/6)/(7TD/6)  (8) 

where  Jj  is  a  Bessel  function  of  the  first  kind.  For  a  cylinder,  V  can  be 

written  as 

V  =  sin(TrL/6)/ (ttLc/6)  (9) 

where  Lc  is  the  length  of  a  cylinder  of  diameter  D  «  6.  Equations  8  and  9 
are  plotted  in  figure  2  to  illustrate  the  salient  features  of  the  visibility 
in  particle  size  measurement. 

Figure  .2  shows  that  for  spherical  particles,  the  visibility  function  is 
not  mono  tonic  for  sizes  greater  than  about  1,05  6.  This  value  fixes  the 
upper  limit  of  the  PSI  size  range.  The  lower  limit  results  from  the  acceptable 
error  in  the  particle  size  measurements.  Since  particle  size  is  a  non-linear 
function  of  visibility,  a  1%  uncertainty  in  visibility  gives  a  1%  uncertainty 
in  D  when  D/6  is  about  1,  but  +  30%  when  D/6  is  about  0.1.  Acceptable  sizing 
uncertainty  and  signal  processor  accuracy  thus  limit  the  low  end  of  the  PSI 
size  range  to  about  0.05  D/6.  Hence,  the  PSI  can  cover  about  20:1  size  range 
for  a  given  6,  Figure  3  shows  the  uncertainty  in  particle  diameter  as  a 
function  of  D/5  for  fixed  uncertainty  values  of  V. 


1.3  Class  II  Particles 


These  particles  may  be  imagined  to  be  small  spherical  lenses  or  mirrors. 

As  the  particle  enters  the  fringe  pattern,  it  causes  the  incident  beams  to 
focus  to  two  point  sources  immediately  in  front  of  the  particle.  The  separation 
of  the  point  sources  is  a  function  of  the  particle  size  and  index  of  refraction. 
The  light  from  the  two  point  sources  generates  a  set  of  interference  fringes 
which  follow  hyperbolic  contours  in  the  plane  of  the  observation  lens.  By 
keeping  the  observation  lens  in  the  plane  perpendicular  to  the  plane  of  the 
beams  and  intersecting  the  bisector,  these  fringes  will  be  nearly  constant 
across  the  lens  aperture.  The  observed  signal  visibility  will  thus  be  a 
function  of  the  separation  of  the  point  sources  and  the  aperture  used  to  detect 
the  refracted  or  reflected  light  from  the  particle.  Analysis  shows  (see 
Appendix  I)  that  the  resulting  visibility  may  be  written  as 


2J1(TTD/6e) 

ttD/6 

e 


(10) 


where  5e  is  the  "effective"  fringe  period  of  the  system.  For  example,  for 
scatter  in  the  near  forward  direction 
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FIGURE  2.  VISIBILITY  AS  A  FUNCTION  OF  PARTICLE  SIZE 
FOR  A  CYLINDER  AND  SPHERE. 
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UNCERTAINTY  IN  PARTICLE  DIAMETER/ UN IT  DIAMETER. 


FIGURE  3.  UNCERTAINTY  IN  SPHERICAL  PARTICLE  DIAMETERS  FOR 
VARIOUS  VALUES  OF  UNCERTAINTY  IN  THE  VISIBILITY. 
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6  *=  4F6  [l  +  tan^B]  ~ — r~ 

e  L  J  m' 

where  F  is  the  F  number  subtended  by  the  scattered  light  receiver,  B  is 
the  observation  angle,  and  m1  is  the  index  of  refraction.  Thus,  the  same 
visibility  function  as  for  Class  I  particles  still  applies.  However,  the 
fringe  period  is  scaled  as  shown  in  Equation  11.  Equation  11  shows  that 
the  size  range  which  can  be  covered  for  a  given  fringe  period  and  particle 
index  of  refraction  can  be  adjusted  through  adjustment  of  the  scattered 
light  collector  F  number. 


2.0  FUEL  SPRAY  CHARACTERISTICS  EXPECTED  TO  AFFECT  PSI  PERFORMANCE 

The  purpose  of  this  section  is  to  develop  simple  parametric  equations 
which  can  be  used  to  illustrate  the  expected  performance  of  a  PSI  system 
under  possible  operational  conditions.  The  parameters  associated  with  a  fuel 
spray  which  may  be  expected  to  have  a  significant  impact  on  PSI  performance 
are: 

f- 

1.  Droplet  number  density 

2.  Droplet  size  range 

3.  Droplet  size  distribution 

4.  Droplet  index  of  refraction 

5.  Material  density  of  the  fuel. 

These  factors  are  discussed  separately  in  the  following  subsections. 

2.1  Droplet  Number  Density 

Consider  figure  4  which  shows  the  geometry  used  to  estimate  droplet 
number  density  as  a  function  of  position  in  the  spray.  Assume  the  spray 
has  a  divergence  angle  given  by  0.  The  mass  flow  rate,  M,  crossing  a  plane 
in  the  spray  is  given  by 


M 


7T  • 

T  P 

O  O  3 


(12) 


where  pQ  =  material  density  of  fuel,  ^3  =  third  moment  of  size  distributions, 
and  N  =  number  of  particles  per  second  crossing  the  plane.  The  droplet  number 
density  can  be  written  as 


=  N 

PN  V  a  At 
~z  s 


(13) 


where  ag  =  cross-sectional  area  of  the  spray.  V2  =  spray  droplet  velocity 
normal  to  ag.  At  =  observation  time  during  which  p^j  is  detected,  and  N  = 
number  of  particles  crossing  ag  during  At.  If  pN  is  constant  during  At 
then  N  =  N/At.  For  some  axial  position,  Z,  in  the  spray  as  is  given  by 
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Fig,  4.  Assumed  spray  nozzle  geometry. 
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TTZ2tan29 


a 
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(14) 


where  0  is  spray  cone  half  angle.  Using  Equations  12  and  14  in  Equation  13 
yields 


M 


N 


7T  2  2 

X  P  UQV  Z  tan  0 
6  o3*s,z 


(15) 


Equation  15  shows  that  can  be  expected  to  vary  inversely  as  Z  .  This 
implies  that  as  Z  decreases,  p^  will  eventually  reach  a  point  where  the  PSI 
cannot  function.  The  factors  affecting  this  limitation  will  be  discussed 
in  detail  in  a  later  section.  It  should  be  clearly  understood  that  Equation 
15  is  an  approximate  relationship  which  assumes  that  the  droplets  have  all 
assumed  the  same  Vz  value.  In  general,  this  will  not  be  the  case.  A  velocity 
distribution  will  exist  which  depends  on  the  droplets  angular  trajectory  and 
variations  in  background  air  and  nozzle  pressure  during  the  sample  interval. 


2.2  Droplet  Size  Range 

A  PSI  can  accommodate  a  size  range  of  approximately  10-20:1  for  any  one 
optical  setting.  If  the  nozzle  produces  particles  over  a  size  range  greater 
than  this  then  it  is  necessary  to  recognize  and  correct  for  the  fact  that  a 
number  of  overlapping  size  ranges  may  be  necessary  to  properly  describe  the 
size  distribution. 


2.3  Droplet  Size  Distribution 

The  shape  of  the  droplet  size  distribution  can  directly  affect  the 
operating  characteristics  of  the  PSI  system.  To  illustrate  these  effects, 
consider  the  following  possible  size  distributions. 

1.  Monodisperse  size  distribution 

2.  Single  mode  log-normal  size  distribution 

3.  Bimodal  log-normal  size  distribution. 


It  is  assumed  for  simplicity  that  the  range  of  the  log-normal  distributions 
is  0  to  °°.  These  distributions  are  chosen  for  simplicity  and  because  they 
are  known  to  characterize  many  droplet  and  particle  processes.  Before 
proceeding,  it  will  be  convenient  to  define  the  moments  of  the  distributions. 
The  kth  sample  moment,  jj^,  of  a  droplet  distribution  expressed  as  a  size 
histogram  is  defined  by 


.1  l 

Nt  i-i 


fiDi 


(16) 


where  N  is  the  total  number  of  measurements  for  the  histograms  and  it  is 
assumed  that  there  are  n  size  increments  with  f^  measurements  in  the  ith  size 
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interval.  is  assumed  to  be  the  largest  diameter  in  the  ith  increment. 

Note  that 


n 

N  =  Z  f 
i=l 


Cl  7) 


Log— normal  distributions  are  characterized  by  two  parameters:  1)  geometric 
mean  diameter,  Dg,  and  2)  the  geometric  standard  deviation  a„.  D_  and  Oc 
are  defined  according  to  ° 


g 


n 


Dg  =  eXp{  N  Z  filnDi} 
S  t  i=l 


(18) 


n 


aK  =  exp{  I  1  f4ln2(D ±/D  )}  . 
B  t  i=l  1  s 


(19) 


A  log-normal  probability  density  distribution  may  be  written  as 

2 


P(D) 


exp 


V2tt  lna  \  LV2  lna 

g  ' 


/  rln(D/vi 

V  Lv2  lna  -I 


(20) 


The  first  four  moments  of  the  distribution  may  be  written  as* 


1  2 

=  Dg  exp (y  In  a  ) 

(21) 

=  D  ^  exp  (2  ln^a  ) 

B  B 

(22) 

=  D  ^  exp (4. 5  ln2a  ) 
g  g 

(23) 

=  D  ^  exp (8  ln2a  ) 

s  s 

(24) 

The  first  four  moments  are  used  to  express  various  mean  diameters  which  have 
appeared  in  the  literature  and  have  been  used  to  characterize  nozzles.  For 
example,  the  volume  to  area  diameter  often  called  the  "Sauter"  mean  diameter, 


Ds,  may  be  expressed  as 


DS  =  y3/y2 


Ds  =  Dg  exP(2*5  ln  ag) 


(25) 

(26) 
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(27) 


The  volumetric  mean  diameter  Dv  is  defined  as 


Dv  =  Dg  exp (1.5  ln2ag)  .  (28) 

Also  found  is  the  mass  mean  diameter  D  given  by 

mm 


D™  ■  Ms  (29) 

D .  =  D  exp (3. 5  ln20  )  .  (30) 

mm  g  g  x  * 

It  is  often  found  that  bimodal  distributions  exist  In  smokes  and  in 
some  sprays.  In  this  case,  the  total  distribution  may  be  viewed  as  the  sum 
of  two  weighted  distributions.  One  distribution  contributes  a  fraction  fN 
to  the  total  number  of  particles  while  a  second  contributes  (l-f)N.  If  such 
a  distribution  is  considered  to  be  the  sum  of  two  log-normal  distributions, 
then  the  resulting  probability  density  will  be  of  the  form 


P(D) 


ln(D/Dgl)  2 
2  lna  , 


+  a-f) 

n/2tt  lna 


In  (D/D  )  2\ 
- (31) 

2  lna  2  J 


where  Dg^,  i«l ,  2  indicates  the  geometric  mean  and  standard  deviation  for 
distributions  1  and  2.  Now  consider  how  the  special  cases  of  interest  affect 
the  expected  number  density. 


2.3.1  Monodisperse  Size  Distribution  of  Diameter  D 

3 

In  this  case,  y^  =  D  •  Equation  15  becomes 

PN(0)  =  2 - - -  (32) 

f  p  D3V  Z2tan20 
6  Ko  -z 

where  PN(0)  Is  the  number  density  for  a  monodisperse  size  distribution.  The 
following  parameters  are  typically  found  in  a  fuel  spray  nozzle. 

D  =  50  pm 

V  =  1524  cm/sec. 

~z 

6  =  70° 

M  *=  667  lb/hour 
M  =  84.11  gm/ sec. 
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Assume  Pq  =  0.9  gm/cc,  then  equation  32  yields 


rns  _  j  ;  10  -i 

2  cc 


(33) 


2.3.2  Single  Mode  Log-Normal  Size  Distribution 

In  this  case,  Equation  15  when  used  with  Equation  23  yields 


-4. 5 In  a 


p„u>  -  _  2 


Me 


g 


(34) 


\  p  D  Z^tan^0 
6  o  g  ~z 


where  p^(l)  is  understood  to  mean  the  number  density  resulting  from  a  single 
mode  log-normal  size  distribution.  If  Dg  is  taken  to  be  50  ym  and  the  parameters 
used  to  compute  p^j(0)  are  used  in  Equation  34,  then 


PBU)  = 


.  -4. 51n  cr 

104  e  8 


cc 


-1 


(35) 


Equation  35  shows  the  result  that  as  a g  increases,  the  number  density  decreases. 
This  should  be  expected  since  a  increasing  corresponds  to  a  broader  distri¬ 
bution  containing  increased  numbers  of  large  particles  carrying  significant 
fractions  of  the  total  mass. 


2.3.3  Bimodal  Log-Normal  Distribution 

For  this  illustration,  it  will  be  assumed  that  the  following  parameters 
apply  to  the  distribution  function  described  by  Equation  31: 


Dgl  =  10  ym 

(.36) 

=  100  ym 

(.37) 

a  .  =  1.64  (lna  .  =  0.5) 
gl  gl 

(38) 

a  0  =  1.105  (lna  0  =  0.1)  . 
g2  g2 

(39) 

To  examine  a  case  that  might  cause  particular  difficulty  with  a  PSI,  consider 
the  situation  where  the  largest  numbers  of  particles  are  in  the  small  size  mode 
but  the  greatest  portion  of  the  mass  is  in  the  large  size  mode.  The  parameters 
in  Equation  31  have  been  chosen  such  that  numerically  92%  of  the  total  number 
of  particles  in  the  distribution  reside  in  the  small  mode.  However,  this  mode 
contributes  only  about  2%  of  the  total  mass.  Let  f  represent  the  numeric 
fraction- of  particles  in  the  first  size  distribution.  Using  Equation  13  and 
32  and  the  parameters  in  Equations  36-39,  there  results  for  the  third  moment 
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of  the  bimodal  size  distribution 


gl  +  Cl-f)Dg23  e 


g2 


(40) 


y3  =  0. 92 (10;  e 


3  4.5(0.25) 


+  0.08(10  )  e 


2.3  4.5(0.01) 


(41) 


8.39  •  104  ym3 


(42) 


=  43.78  ym 


(43) 


Using  Equation  42  along  with  the  parameters  used  in  section  2.3.1  (Equation  32) 
yields  for  the  number  density 


where  y3  has  been  substituted  for  D^.  Equation  44  shows  that  roughly  a  50% 
increase  in  results  from  the  parameters  chosen  for  the  bimodal  distribution 
even  though  the  volumetric  mean  diameter  is  nearly  the  same  as  that  for  the 
monodisperse  case.  Figure  5  plots  the  numeric  distribution  for  the  parameters 
assumed  to  obtain  equation  44.  As  the  figure  shows,  the  PSI  would  need  to 
count  large  numbers  of  particles  very  accurately  if  accurate  pN  and  hence 
mass  concentrations  are  to  be  obtained. 


2.4  Droplet  Index  of  Refraction 

There  are  two  observation  modes  which  have  been  found  to  be  useful  with 
the  PSI.  These  are  called  the  "on  axis"  (for  Class  I  particles)  and  "off  axis" 
(for  Class  II  particles)  viewing  modes.  The  "on  axis"  viewing  mode  requires 
that  the  scattered  light  collection  system  be  centered  on  the  bisector  between 
the  beams  forming  the  PSI  sample  volume.  (A  discussion  of  the  theoretical 
interpretation  and  operating  principles  of  the  PSI  in  this  mode  is  given  in 
reference  5.)  Figure  6  illustrates  this  viewing  mode.  In  this  configuration, 
it  has  been  found  that  the  PSI  size  parameter  is  independent  of  the  droplet 
index  of  refraction  when  the  particle  absorbs  a  small  fraction  of  the  incident 
light.  In  this  case,  the  visibility  function  is  monotonic  when  the  particle 
diameter  is  less  than  about  1  fringe  period,  6,  in  diameter.  Hence,  the  PSI 
measurement  is  independent  of  index  of  refraction  and  covers  a  size  range 
between  roughly  0.05-1.0  6  when  used  in  the  axis  viewing  mode. 

In  the  off  axis  viewing  mode  shown  in  figure  7,  the  PSI  signal  is 
obtained  by  collecting  the  scattered  light  at  some  angle  off  the  bisector 
between  the  beams.  Appendix  I  discusses  how  index  of  refraction  affects 
the  PSI  signal  and  presents  a  simple  model  for  understanding  why  this  is  so. 

In  the  off  axis  viewing  mode,  the  particle  size  range  is  shifted.  A  10-1 
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PCD) 


Fig.  5.  Probability  density  distribution  for  the 
assumed  bimodal  log-normal  distribution. 
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Fig,  7.  Example  of  possible  PSI  off  axis  viewing  directions. 
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or  20-1  size  range  is  still  covered  but  the  location  of  the  size  range  will 
depend  on  the  droplet  index  of  refraction  and  the  F  number  of  the  scattered 
light  collector.  When  droplets  are  opaque  the  discussion  in  reference  5  applies. 


2.5  Material  Density  of  the  Fuel 

The  PSI  produces  sufficient  information  to  compute  the  time  average 
mass  concentration  <C>  or  mass  flow  rate  M  at  a  particular  point  in  the 
spray.  <C>  can  be  written  as 


C  ’  f  P„  «5> 

where  <U3Pn>  is  the  time  averaged  value  of  the  third  moment  and  number  density 
of  the  size  distribution.  pQ  is  the  material  density  of  the  fuel.  Equation 
45  shows  that  if  absolute  values  of  <C>  are  computed,  then  pQ  must  be  known. 
Equation  15  also  shows  that  if  M  is  computed  from  PSI  data  that  pQ  is  a 
fundamental  parameter  which  must  be  known. 


3.0  PSI  DESIGN  LIMITATIONS  IMPOSED  BY  FUEL  SPRAY  REQUIREMENTS 

Section  2.0  discussed  the  spray  characteristics  which  will  affect  the 
performance  and  accuracy  of  a  PSI.  In  this  section,  the  results  of  section 
2.0  are  supplied  to  the  response  parameters  which  define  the  operational 
characteristics  of  the  PSI. 


3.1  Number  Density  Limitations 


3*1.1  Data  Acquisition  Limitations 

Determination  of  number  density  with  any  optical  instrument  is  difficult 
and  subject  to  numerous  errors.  This  is  especially  true  for  a  PSI.  In 
estimating  number  density  with  a  PSI,  at  least  two  different  approaches  have 
been  developed.  The  first  approach  is  to  attempt  to  make  the  sample  volume 
so  small  that  virtually  no  signals  are  ever  seen  from  more  than  one  particle. 
Then  counts  of  individual  particles  can  be  correlated  with  the  total  volume 
of  fluid  sampled.  This  approach  can  only  be  valid  when  particle  number  density 
is  so  low  that  the  average  number  of  particles  in  the  sample  volume  is  much 
less  than  1.  In  sprays  where  the  number  density  may  fluctuate,  erroneous 
counts  may  result  from  the  presence  of  more  than  one  particle  in  the  sample 
volume.  The  second  approach  recognizes  this  fact  and  attempts  to  make  use 
of  the  multiple  particle  detection  capability  of  the  PSI  and  statistical 
models  to  estimate  particle  number  density.  Appendix  II  presents  a  discussion 
of  the  model  employed  in  the  second  approach  for  estimating  number  density. 

This  is  the  method  used  by  UTSI  to  estimate  number  density  and  it  has  been 
found  to  yield  surprisingly  good  results.  The  PSI  parameter  used  to  estimate 
number  density  is  called  the  acceptance  ratio  A.  A  is  the  ratio  of  the  number 
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of  single  particle  signals  (with  aperiodicity  less  than  50%)  detected  to 
the  total  number  of  signals  detected.  Appendix  II  shows  that  A  can  be 
written  as 


A  =  1 


-[»- 


V73**2  V7V 

exp  - - - )  (i_eXp 


.  1/3.  2  1/2 

"PN  Az  1 
■>  a-“p  — — > 


<46) 


where  5  is  the  fringe  period  of  the  system  and  Ax,  Ay,  Az  are  the  rectangular 
spatial  measures  of  the  PSI  sample  volume.  For  purpose  of  illustration,  we 
will  consider  two  extremes  in  the  optical  viewing  geometry.  These  two  cases 
are  illustrated  in  figure  8.  In  the  first  case,  we  consider  an  on  axis 
viewing  geometry  in  which  Ax  =  Ay  and  AxAy  is  the  sample  volume  cross-section 
normal  to  the  bisector  between  the  beams.  In  this  geometry  Az  >>  Ax,  Ay. 
Equation  46  then  becomes 


A  1/3.  2 
rpN  ^ x 

A  =  exp  ( - ^ - ) 


(47) 


If  it  is  assumed  that  the  particles  generate  a  certain  number  of  cycles 
in  the  signal,  N^,  before  they  are  measured,  then  Ax  can  be  written  as 


Ax  -  N  6  . 

JL 


Using  Equation  48  in  Equation  47  yields 


A  = 


tW'\\ 

exp  ( ) 


TT 


(48) 


(49) 


In  the  second  case,  it  is  assumed  that  the  scattered  light  receiver  is 
positioned  90°  off  the  bisector  between  the  beams.  This  case  would  correspond 
to  roughly  the  smallest  possible  sample  volume.  In  this  case,  Ax  =  Ay  =  Az 
and  use  of  Equation  48  reduces  equation  46  to 


A  =  1 


,,  r(W‘  3/2 

(1  -  exp  ( - - - )) 


(50) 


Equations  49  and  50  indicate  the  extremes  in  A  that  could  be  obtained  for 
two  possible  viewing  geometries.  As  Appendix  II  points  out,  A  can  be 
increased  by  using  slit  apertures  to  reduce  Ax.  In  this  way,  increased 
number  density  response  can  be  obtained.  To  illustrate  the  PSI  operating 
characteristics  indicated  by  Equations  49  and  50,  assume  NL  =  10  and 
6  =*  30  ym.  (N^  =  10  has  been  found  to  yield  reasonable  results  in  previous 
applications  and  6  =  30  ym  will  yield  a  size  range  of  1.5  —  30  ym  for  on 
axis  viewing  and  a  15-300  ym  size  range  for  the  case  of  90°  off  axis  observa¬ 
tion  and  an  F/3.5  receiver.  See  Table  I  in  Appendix  I.)  Figure  9  plots  A 
for  the  two  cases  of  interest.  The  figure  shows,  for  example,  that  when 
=  9  that  the  PSI  will  measure  between  1  and  2%  of  the  total  number  of 
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Fig.  8.  Exploded  View  of  Beam  Intersection  Point  Showing  How 
Receiver  Lens  Orientation  Limits  PSI  Sample  Volume. 


-  -  -  OFF  AXIS 

OBSERVATION 


A 


Fig.  9.  Example  of  acceptance  ratio  as  a  function  of 

number  density  for  on  axis  and  off  axis  viewing 
geometry.  Z  values  correspond  to  Z  positions 
(cm)  in  a  50  \im  diameter  moriodisperse  spray  with 
mass  flow  of  667  lbs/hr.  (84.11  gm/sec)  . 
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particles  detected  and  that  90°  off  axis  observation  does  not  give  that 
much  better  response  for  than  the  on  axis  geometry.  Also  plotted  in 
figure  9  are  the  corresponding  Z  positions  in  the  model  spray  described 
for  a  monodisperse  size  distribution,  (Note,  however,  that  in  case  2 
spatial  resolution  will  be  significantly  better  than  for  case  1.) 


Past  experience  has  shown  that  when  A  reaches  a  value  of  0.1  that 
data  acquisition  with  the  PSI  is  too  slow  to  obtain  a  statistically 
significant  data  set  in  a  reasonable  length  of  time.  If  A  =  0.1  is 
chosen  as  the  limiting  criteria  for  data  acquisition,  then  equations 
49  and  50  can  be  solved  to  specify  the  required  and  6  values  to  obtain 
a  reasonable  optical  geometry  for  an  estimated  limiting • number  density. 
Solving  Equations  49  and  50  for  these  cases  yields 


N 


=  3,78  *  10' 
N  6S3 

max  L 


case  1 


_  6,04  *  10 

N  '  N,V 

I  max  L 


case  2 


(51) 


(52) 


Suppose  the  spray  produced  p^  values  up  to  10  particles  cc  .  =  10 

is  approximately  as  small  as  desirable  for  successful  operation  of  the  PSI. 
Equation  51  in  that  case  yields  a  5  value  of  7.23  ym  for  the  on  axis  case 
while  Equation  52  shows  that  the  smallest  acceptable  fringe  period  in  case  2 
is  8.45  ]Jm.  This  means  that  for  p^  =  10  cc“^,  the  size  range  achievable 

for  the  on  axis  case  is  about  0.36  to  7.23  ym.  In  the  90°  off  axis  case, 
Table  I  in  Appendix  I  shows  that  an  F/12.6  scattered  light  receiver  would 
yield  a  size  range  of  approximately  30-300  ym  for  a  p ^  value  of  10^  cc“l. 
While  these  calculations  seem  to  indicate  that  the  PSI  would  function  for 
number  densities  of  the  order  of  10^  cc"l  over  the  required  size  range  for  a 
90°  off  axis  viewing  geometry,  the  reader  should  be  cautioned  that  UTSI  has 
not  attempted  such  measurements  (90°  off  axis)  on  number  densities  this  high 
and  cannot  be  certain  that  such  measurements  are  possible.  However,  measure¬ 
ments  have  been  obtained  for  hygroscopic  smokes  with  reasonable  success  in 
the  on  axis  case  with  fringe  periods  similar  to  those  calculated  for  case  1. 
These  results  are  described  in  Appendix  IIL  Table  III.l  summarizes  the 
results  of  these  computations. 


3.1.2  Optical  Transmission  Limitations 

Optical  transmission  is  a  fundamental  limitation  imposed  by  particle 
number  density  of  PSI  operation.  In  this  area,  the  spray  density  and  droplet 
scattering  cross-section  must  be  sufficiently  small  that  adequate  beam 
transmission  can  be  obtained  through  the  spray.  During  the  work  reported 
In  Appendix  III,  It  was  found  (but  not  reported)  that  when  the  optical 
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PSI  CHARACTERISTICS  REQUIRED 
FOR  RESPONSE  TO  pN  =  106  cc“ 
FOR  AN  ACCEPTANCE  RATIO  OF  0 
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Uncertainty) 


transmittance  fell  below  about  0.8,  the  PSI  data  rate  fell  to  an 
unusably  low  value.  Thus,  in  this  section,  PSI  operational  limitations 
are  examined  in  light  of  this  transmittance  limit.  The  Beer-Bouaher 
transmission  law  may  be  written  as 

T  =  exp(-apNL)  (53) 

where  a  is  the  mean  extinction  cross-section  and  L  is  the  transmission 
path  length.  For  a  polydispersion,  O  can  be  written  as 

—  f  2  7T 

0=jQEWD2l-j—  C54) 

o 

where  QE  is  the  extinction  efficiency  and  dp^j/p^  is  the  numeric  size 
distribution.  The  extinction  efficiency  is  defined  as  4a/uD^  where  O  is 
the  extinction  cross-section. ^  When  D  »  X  (the  optical  wavelength)  QE  -  2 
independent  of  D.^  This  is  assumed  to  be  the  case  for  fuel  spray.  Then 
Equation  54  reduces  to 


®  -*  2  ^2  (55) 

where  QE  =  2  has  been  substituted  into  Equation  54  and  p2  understood  to 
represent  the  second  moment  of  the  size  distribution. 


’1 


1.42 


10 


max. tran. 


y„  L 
M2ym 


cc 


-1 


C56) 


2 

where  y2ym  is  the  second  moment  of  the  size  distribution  expressed  in  ]im 
and  L  is  in  cm. 


For  the  monodisperse  size  distribution  described  by  Equation  33, 
Equation  56  reduces  to 


pN(°) 


5.68  •  10' 


cc 


-1 


max. tran. 


(57) 


Comparison  of  Equations  33  and  57  shows  that  two  competing  effects  exist 
t  to  limit  the  PSI  operation  -  either  the  number  density  may  be  too  high 

to  obtain  a  useable  data  rate  or  the  transmission  path  through  the  spray 
may  be  too  long  even  at  reduced  number  densities.  An  estimate  of  the 
maximum  value  of  Z  in  the  monodisperse  spray  can  be  obtained  by  letting 
L  =  2Ztan0  and  equating  Equations  33  and  57.  The  result  it 

Z(0)  =  38.7  cm  (58) 

max 
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For  the  single  mode  log-normal  distribution  with  Dg  =  50  ym,  Equations 
22  and  56  show  that  Pn C 1 ^max , tran<  is  given  by 


pNa> 


5.68  •  103  “2  ln  °g 

- - -  e  ° 


(59) 


max. tran. 


Using  Equations  33,  35,  57,  and  59,  Z  (1)  thus  can  be  written  as 

max 


v°> 

PN0) 


,  Mil  2  lA« 

Z(0)  e 


max, tran. 


(60a) 


-2.5  In  a 


Z  (1)  =  Z  (0)e 
max  max 


(60b) 


Hence,  for  the  single  mode  log-normal  distribution  the  transmission  path 
decreases  because  the  particle  number  density  increases.  For  example,  if 
lnCg  =  0.5,  then  Zmax(l)  is  about  50%  of  ^^(0).  For  the  bimodal  distri¬ 
bution  specified  by  Equations  36-39  and  a  small  particle  fraction  of  92%, 
the  second  moment  is  computed  to  be 

U2pm  ‘  9‘68  '  wW  <61> 


fmax. tran. 
bution: 


as  expressed  in  Equation  56  thus  becomes  for  the  bimodal  distri- 


P„(2) 


1.46  •  10 


max. tran. 


(62) 


Using  Equation  44  with  Equation  62  yields 

Z  (2)  ^  22.6  cm  (63) 

max 

These  examples  clearly  show  that  the  size  distribution  characteristics  of  the 
spray  can  significantly  affect  the  operational  performance  of  the  PSI.  As 
a  result,  specification  of  number  density  operating  characteristics  for  the 
PSI  must  be  dependent  on  the  spray  to  be  measured. 


3.2  Velocity  Limitations 

The* PSI  signal  processing  system  imposes  velocity  limits  which  depend 
on  fringe  period  used  in  the  optical  system.  The  minimum  velocity  measurable 
by  the  PSI  can  be  derived  from  the  equation 
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V  -<  =  <Sv_  . 

~z  min  Dmin 


(64) 


where  is  the  minimum  Doppler  frequency  acceptable  by  the  data 

acquisition  system,  UTSI  has  supplied  an  optical  system  which  produces 
a  minimum  fringe  period  of  6  ym.  vDmin  for  t*ie  signal  processing 

electronics  is  limited  by  the  high  pass  filter  and  is  30  KHz,  Hence, 

V  .  =  0.18  m/sec,  (65) 

~z  min 

The  maximum  Doppler  frequency  obtainable  with  the  signal  processing 

electronics  is  approximately  3.5  MHz  and  is  limited  by  acceptable  signal 

distortions  introduced  by  the  full  wave  rectifier  as  the  signal  frequency 

increases.  This  leads  to  a  maximum  measurable  velocity  V  for  the  small 

•  max 

particle  size  range  (6=6  ym)  of 

V  =21  m/sec.  (66) 

~z  max  . 

The  largest  fringe  period  for  the  PSI  system  is  100  ym.  In  this  case, 
the  maximum  and  minimum  velocities  correspond  to 

V  ,  =  3  m/sec.  (67) 

-vz  min 


V  =  350  m/sec. 

~z  max 


(68) 


However,  as  figure  9  shows,  the  maximum  number  density  will  correspond  to 
approximately  10^  cc~l  for  an  intermediate  fringe  period  of  30  ym.  For  a 
100  ym  fringe  period  the  maximum  number  density  will  correspond  to  about 
370  cc~l  as  estimated  using  Equation  51.  Table  III.l  summarizes  the  velocity 
limitations.  It  should  be  noted  that  if  only  velocity  and  number  density 
are  to  be  measured  then  the  maximum  measurable  Doppler  frequency  is  about 
35  MHz (upper  frequency  response  of  the  Doppler  processor)  leading  to  maximum 
velocities  of  210  m/sec  (6=6  ym)  and  3500  m/sec  (6  =  100  ym) . 


3.3  Particle  Size  Range  Limitations 

It  has  already  been  indicated  in  previous  sections  that  the  particle 
size  range  is  a  function  of  viewing  geometry  and  fringe  period.  Appendix 
I  and  reference  5  detail  how  these  factors  directly  affect  the  size  range. 
The  ultimate  limitations  of  the  size  range  depend  on  two  factors: 
acceptable  error  for  the  small  size  measurements  and  the  monotonic  limit 
in  the  visibility  when  the  particle  sizes  become  larger  than  a  fringe 
period  (or  equivalent  fringe  period).  In  the  latter  case,  the  error  in 
diameter  will  be  of  the  same  order  as  that  for  the  visibility.  On  the 
other  hand,  when  the  diameter  to  fringe  period  ratio  is  about  0.1,  the 
uncertainty  in  diameter  is  about  30  times  that  of  the  visibility  (see 
figure  3).  As  a  result,  a  10:1  size  range  is  achievable  with  these  devices 
for  a  30%  uncertainty  in  the  smallest  sizes  and  20:1  (as  extrapolated  from 
figure  3)  when  errors  greater  than  about  60%  are  acceptable. 
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4.0  CALIBRATION  STANDARDS  AND  METHODS 


Calibration  of  the  PSI  system  falls  into  two  categories: 

1)  optical,  and  2)  electronics.  The  electronics  system  calibration 
is  dependent  on  a  straightforward  operational  procedure  which  is 
detailed  in  the  PSI  Electronics  Operations  Manual.^  The  optical 
system  calibration  procedures  are  discussed  here. 


4.1  Optical  Calibration  Procedures 

Three  different  calibration  procedures  should  be  performed  on 
the  optical  system.  These  are: 

1.  Fringe  period  calibration 

2.  Checks  of  visibility  as  a  function  of  particle  size 

3.  Depth  of  field  and  probe  volume  diameter  calibration. 

The  following  sections  describe  how  these  calibrations  have  been 
performed  in  past  tests  and  the  results  of  such  calibrations  are 
described. 


4.2  Fringe  Period  Calibration 

The  electronics  measure  the  Doppler  time  period,  T,  or  Doppler 
frequency  1/t,  Vp,  and  the  deviations  in  frequency  for  a  set  of  measure¬ 
ments.  Equation  2  shows  that  the  fringe  period  6  can  be  written  as 

5  »  TV  .  C69) 

Hence,  for  a  known  velocity  and  measured  T,  6  can  be  computed.  The  known 
velocity  is  typically  generated  by  a  synchronous  motor  which  spins  a 
small  wire  mounted  on  a  wheel  attached  to  the  motor  shaft.  The  radius,  r, 
from  the  center  of  the  shaft  to  the  center  position  of  the  wire  is  measured 
with  a  micrometer  calipers  to  the  nearest  0.001u.  A  typical  value  of  r 
is  2.125  +  0.001”.  Commercial  power  is  regulated  to  a  frequency  of  60  Hz 
with  precision  and  accuracy  typically  better  than  0.1%.  By  placing  the 
wheel  such  that  its  tangent  is  perpendicular  to  the  fringe  planes  the  wire 
velocity  can  be  computed  from 

V  =  2irr60  cm/sec.  (70) 

where  r  is  expressed  in  cm.  By  using  the  PSI  system  to  measure  T,  6  can  be 
computed  from 


6  =  120  TTrT  cm  (71) 

where  T  is  expressed  in  seconds.  0.5%  uncertainty  or  better  in  the 
calibrated  6  values  can  be  achieved  using  this  method  to  determine  the  fringe 
period. 
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4.3  Visibility  as  a  Function  of  Particle  Size 

To  obtain  this  calibration,  UTSI  normally  uses  a  vibrating  orifice 
monodisperse  generator  to  produce  a  stream  of  droplets  with  sizes  in  the 
range  of  60  to  100  micrometers  or  greater.  The  droplet  diameter  can  be 
computed  with  less  than  5%  uncertainty  from  the  fluid  flow  rate  and 
orifice  vibration  frequency.  To  further  reduce  this  uncertainty,  a 
microscope  and  strobe  light  system  is  used  to  photograph  the  droplets  at 
the  same  time  they  were  measured  with  the  PSI.  Because  the  droplets  are 
typically  separated  by  4-5  diameters,  multiple  particles  may  be  present 
in  the  probe  volume  from  the  initial  output  stream.  This  is  overcome  by 
using  a  droplet  charging  and  deflection  system.  By  charging,  say,  every 
64th  drop  and  passing  the  drop  stream  through  a  high  static  electric 
field  (typically  2000  volts/cm)  the  charged  drops  can  be  directed  through 
the  sample  volume  while  the  main  stream  follows  a  different  trajectory. 
Figure  10  shows  photographs  taken  with  the  microscope  system  of  the 
droplet  stream.  The  estimated  particle  diameter  from  flow  rate-frequency 
calculation  and  the  visual  diameter  are  compared  under  each.  Figure  11 
shows  the  particle  deflection  system  and  experimental  arrangement. 


4.4  Probe  Volume  Calibration 

In  this  section,  the  pertinent  parameters  required  for  probe  volume 
calibration  are  described.  The  value  of  the  probe  volume  is  required  if 
particle  number  density  estimates  or  volume  normalized  size  distributions 
are  to  be  obtained.  Two  possible  probe  volumes  should  be  considered. 

These  will  be  referred  to  as  the  "High  Density  Probe  Volume"  (HIRO)  and 
the  "Low  Density  Probe  Volume"  (LORO).  The  HIRO  volume  corresponds  to 
the  volume  defined  by  the  localized  region  where  the  beams  cross.  It  is  a 
function  of  the  incident  beam  intensities,  optical  viewing  system  and 
particle  scattering  cross-section,  and  signal  processing  logic.  Table 
IV  shows  the  basic  geometry  of  some  possible  HIRO  probe  volumes  commonly 
encountered  in  PSI  application  and  the  formulas  used  to  compute  these 
volumes.  The  LORO  probe  volume  is  the  cylindrical  volume  defined  by  the 
cross-sectional  area  of  the  HIRO  volume  normal  to  the  mean  flow  direction 
and  a  height,  h,  defined  by  the  mean  velocity  and  the  total  sample  time, 
i.e.  h  =  vzt.  Table  IV. 1  illustrates  the  LORO  probe  volumes  and  the 
appropriate  formulas  for  computing  the  volumes  associated  with  a  particular 
viewing  geometry.  The  calibration  factor  common  to  both  volumes  is  the 
maximum  number  of  signal  cycles,  NQ,  detectable  by  the  system.  Once 
this  value  is  known,  then  the  HIRO  volume  can  be  directly  estimated.  In 
the  LORO  case,  the  length  of  the  sample  volume  for  a  given  maximum  width 
can  be  a  calibration  factor.  Table  IV. 2  illustrates  depth  of  field  calibra¬ 
tion  results  using  the  arrangement  shown  in  figure  11.  The  droplet  stream 
was  positioned  down  the  probe  volume  axis  until  no  more  signals  were  read 
by  the  processor.  The  tabulated  parameter  was  then  changed  and  the 
measurement  repeated.  The  calculated  values  with  respect  to  the  measured 
values  were  within  +  1  cycle  of  that  predicted  using  the  formulas  used  to 
compute  the  probe  volumes.  The  choice  of  a  particular  probe  volume  is 
somewhat  arbitrary.  However,  a  reasonable  "rule  of  thumb"  is  to  use  the 
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MONODISPERSE  CALCULATED  DIAMETER  =  68  /till 
VISUAL  MICROSCOPIC  DIAMETER  =  74  fj m 


MONODISPERSE  CALCULATED  DIAMETER  =  59/flll 
VISUAL  MICROSCOPIC  DIAMETER  =  64/011 


MONODISPERSE  CALCULATED  DIAMETER  =  60 /tm 
VISUAL  MICROSCOPIC  DIAMETER  =  57/011 


MONODISPERSE  CALCULATED  DIAMETER  =  64  U 
VISUAL  MICROSCOPIC  DIAMETER  =  72M 


Figure  10*  Examples  of  microphotographs  of  the  droplet  stream  from  the 
vibrating  orifice  monodisDerse  generator. 
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FIGURE  11.  MONODISPERSE  GENERATOR  WITH  SYNCHRONIZED 
PARTICLE  SEPARATION  AND  STROBOSCOPIC 
MICROSCOPE  OBSERVATION. 
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TABLE  IV 

HIRO  PROBE  VOLUME  SUMMARY 
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NOTE:  N  -  Number  of  interference-fringes 
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TABLE  IV. 2 

DEPTH  OF  FIELD  CALIBRATION 
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I)  PKT  CAIN  -  8.3  *  lCr<  ~  )  (WHERE  V  •  PMT  VOLTAGE  SETTINC) 


HIRO  volume  when  signal  separation  is  much  less  than  signal  duration 
and  LORO  when  signal  separation  is  much  greater  than  signal  duration 


4,5  Histogram  Weighting  Factors 

A  fundamental  use  of  the  probe  volume  calibration  is  for  the 
estimation  of  a  volume  normalized  particle  size  distribution.  This  is 
necessary  to  account  for  the  fact  that  the  size  of  the  probe  volume 
depends  on  particle  scattering  cross-section  (all  other  parameters 
being  constant)  and  the  variation  in  fringe  contrast  throughout  the 
probe  volume.^ 

This  concept  is  straightforward  if  one  imagines  a  large  box  containing 
a  broad  range  of  particle  sizes.  If  a  scattered  light  detector  with  limited 
sensitivity  such  as  the  eye  or  a  lens-phototube  arrangement  is  used  at  a 
fixed  position  to  observe  the  particles,  the  smaller  particles  will  be 
detectable  over  a  smaller  distance  away  from  the  detector  than  the  larger 
particles.  Hence,  if  a  particle  size  distribution  were  determined  by  the 
detector  the  smaller  particles  would  be  counted  in  a  smaller  volume  than 
the  larger  particles  —  the  number  of  small  particles  would  in  effect  be 
undersampled  by  the  ratio  of  volumes  in  which  the  particles  could  be  detected. 
Multiplying  the  small  particle  count  by  this  ratio  would  thus  weight  the 
count  to  correspond  to  the  same  volume  as  the  larger  particles.  It  should 
also  be  borne  in  mind  that  the  probe  volumes  described  thus  far  have  been 
somewhat  arbitrary  (in  the  sense  that  the  threshold  level  of  the  detector 
was  not  specified)  for  purposes  of  modeling  the  function  of  the  instrument. 


4.5.1  Fringe  Contrast  Weighting  Factors 

The  fringe  contrast  weighting  factors  account  for  the  possibility 
that  the  measured  signal  visibility  may  be  lower  than  that  predicted 
on  the  basis  of  its  particle  size  relative  to  a  fringe  period.  This 
reduction  arises  because  of  the  incident  fringe  contrast  reduction  away 
from  the  geometric  center  of  the  probe  volume  (the  interfering  beams 
no  longer  have  equal  amplitudes).  Table  IV. 3  lists  these  weighting 
factors  as  a  function  of  D/6,  and  the  length  of  the  probe  volume, 
m,  relative  to  the  e”^  length,  2bQ/a,  where  a  is  the  angle  between 
the  beams.  When  the  length  of  the  sample  volume  is  determined,  then 
these  weighting  factors  should  be  used  for  correction  up  to  that  length 
only.  For  example,  suppose  the  maximum  Z  axis  trajectory  which  could  be 
observed  was  0.2  times  the  e~2  length.  Then  no  weighting  factors  would 
be  needed  since  Table  IV. 3  shows  all  weighting  factors  for  all  observable 
trajectories  0  <_  m  <_  0.2  are  1,  i.e.  the  histogram  bins  containing  sizes 
corresponding  to  low  visibilities  would  be  known  not  to  contain  particle 
counts  resulting  from  small  particles  with  low  visibilities  due  to  a  Z 
axis  trajectory  corresponding  to  m=l.  This  case  would  occur  for  90° 
viewing.  On  the  other  hand,  if  all  Z  trajectories  could  be  detected 
and  measured  then  the  bins  corresponding  to  low  visibilities  would  be 
multiplied  by  the  appropriate  weighting  factor  to  reduce  an  incorrect 
number  of  counts  resulting  from  point  particles  passing  near  m=l 
values. 
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TABLE  IV. 3 

FRINGE  CONTRAST  WEIGHTING  FACTORS 
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4.5.2  Particle  Size  Weighting  Factors 

These  weighting  factors  can  be  extremely  complicated  to  compute. 
Experience  has  shown  that  reasonable  results  are  obtained  using  the 
following  approximate  results. 


4. 5. 2.1  On  Axis  Measurements 


* 

In  this  case, 
conditions  using: 

the  weighting  factors  are  computed  for  HIRO  and  LORO 

tf 

HIRO: 

W.  =  WrG  D  2/G.D  2 
i  f  o  o  i  i 

(72) 

LORO: 

W.  =  W,G  %  /g\, 
i  f  o  o  i  i 

(73) 

where  W^  is  the  weighting  factor  for  the  ith  histogram  bin,  is  the 
largest  particle  in  the  ith  histogram  increment,  DQ  is  the  largest  particle 
in  the  nlargest  resolvable  size’1  histogram  bin,  and  G0,  G £  are  the 
scattering  gain  function  for  the  DQ,  th  particle  size  and  W£  is  the 
fringe  contrast  weighting  factor.  These  functions  may  be  computed  from 
the  equation 

G  =  J  2(3’e)  -  J  2(3’)  +  J  2C3'e)  -  J  2(3')  (74) 

i  o  o  i  i 

Bf  =  7rD±AF  (75) 

where  JQ,  are  zeroth  and  first  order  Bessel  functions  of  the  first  kind, 
e  is  the  fraction  of  the  collection  aperture  radius  obscured  by  a  circular 
beam  stop  and  F  is  the  receiver  F  number.  Reference  4  provides  additional 
discussion  and  insight  into  the  computation  of  the  weighting  factors  and 
detailed  tabulation  of  G±  for  the  PSI  system  were  provided  upon  system 
delivery. 


4. 5. 2.2  Off  Axis  Measurements 

The  computations  in  this  case  (both  for  Class  I  and  Class  II  particles) 
are  generally  more  complicated.  In  this  case,  both  simple  models  and  a 
sophisticated  Mie  scatter  program  developed  at  UTSI  for  this  purpose  should 
be  used  as  deemed  necessary  from  calibration  measurements  with  the  mono- 
disperse  generator.  A  discussion  of  some  aspects  of  this  program  are 
described  in  Appendix  IV.  In  the  special  case  of  droplets  greater  than 
about  10  micrometers  in  diameter  and  observation  angles  less  than  30°, 

Gi  -  GQ  for  all  viewing  angles  and  F  numbers.  Equations  72  and  73  apply 
with  G0/Gi  =  1.  Reference  4  provides  a  general  discussion  of  the  theory 
behind  the  development  of  the  probe  volume  as  a  function  of  particle  size 
and  the  resulting  weighting  factors. 
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4.6  Glass  Bead  Calibration 


Glass  beads  offer  a  convenient  means  for  in  house  calibration 
tests.  Such  particles  should  be  understood  to  be  poor  measures  of 
system  performance  because  of  non-spheristy  and  inhomogenities . 
However,  these  particles  can  be  used  as  an  indication  of  proper  system 
performance.  UTSI  uses  such  calibration  beads  in  4  size  ranges 
(ranging  from  roughly  60  ym  to  several  hundred  microns)  and  a  device 
to  disperse  the  beads  through  the  sample  volume  similar  to  that  shown 
in  figure  12.  These  may  be  purchased  from  Duke  Standards  Inc.  or  the 
Cataphote  Corporation.  (See  the  list  of  references  for  appropriate 
addresses. ) 


5.0  DATA  SOFTWARE  PROGRAMS 

Normally  the  data  reduction  programs  used  in  the  UTSI  PSI  are 
recorded  on  floppy  disks  and  are  used  to  manipulate  data  recorded  in 
memory.  Since  program  requirements  tend  to  vary  with  applications, 
the  development  of  a  specialized  program  for  BRL  is  offered  as  an 
option.  A  program  commonly  used  by  UTSI  is  included  in  order  to 
illustrate  system  capability.  It  uses  a  21  bin  size  histogram  and 
computes  numeric  and  mass  probability  density  functions-  It  is 
possible  (and  has  been  used  in  previous  research)  for  forming  histograms 
of  up  to  100  bins  or  to  distribute  the  bins  logarithmically  in  size. 


5.1  Current  UTSI  Data  Reduction  Program 

In  this  section,  examples  of  current  UTSI  PSI  software  programs 
and  their  interpretation  are  presented.  Section  5.3  provides  data 
examples  from  a  program  used  in  a  number  of  obscurant  tests.  Section 
5.3  also  shows  the  results  from  a  program  which  totals  or  summarizes 
data  from  a  number  of  individual  data  files. 


5.1.1  Interpreting  PSI  Data  Printouts 

The  following  is  intended  to  instruct  the  reader  in  interpreting 
data  printouts  from  the  Particle  Sizing  Interferometer.  In  examining 
a  data  sheet,  the  reader  will  observe  that  the  information  is  presented 
in  eight  blocks  with  a  space  between  each  block.  The  first  block  of 
information  which  begins  !,Long  Fringe  Count"  specifies  the  logic 
parameters  which  were  used  to  control  the  signal  processor  during  data 
acquisition,  to  specify  the  size  of  the  sample  which  was  recorded  in 
memory,  and  to  indicate  a  run  number  which  is  also  operator  generated. 
The  second  block  contains  three  lines  for  operator  generated  messages 
which  pertain  to  experimental  conditions  or  instrument  parameters  not 
covered  in  the  first  block.  The  third  block  of  data  which  begins  "Data 
Read  From  Diskette  File  ..."  contains  the  following  information.  The 
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Fig.  12.  Method  for  Glass  Bead  Dispersal  Through  PSI  Sample  Volume. 
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diskette  file  number  identifies  the  location  of  the  complete  data  set 
recorded  on  the  magnetic  disk.  The  second  line  provides  the  mean  particle 
velocity  in  m/s,  the  standard  deviation  in  m/s  and  the  turbulence 
intensity  (velocity  distribution  standard  deviation  relative  to  the  mean 
velocity)  in  percent.  The  line  containing  "frequency”  gives  the  mean 
Doppler  frequency  from  which  the  mean  velocity  was  computed.  "Aperiodicity 
Accepted  and  aperiodicity  Failures"  indicate  the  number  of  signals  in 
the  total  sample  size  which  the  control  logic  indicated  was  either  that 
resulting  from  a  single  particle  (Aperiodicity  Accepted)  or  those  resulting 
from  two  or  more  particles  (Aperiodicity  Failures) •  High  pass  failures 
indicate  the  number  of  signals  which  had  frequencies  that  were  above  the 
high  pass  filter  limit  set  in  block  one.  Acquisition  time  is  the  time  in 
seconds  required  to  obtain  the  sample  size  being  analyzed.  The  next  line 
shows  IRIG-B  time  at  which  the  sample  acquisition  was  initiated.  The 
next  block  indicates  the  number  of  measurements  which  passed  tests  imposed 
on  the  signal  in  order  that  it  be  accepted  for  a  size  measurement.  The 
first  two  parameters  indicate  the  number  of  particles  which  generated 
signals  that  were  beyond  the  resolution  capability  of  the  instrument 
(either  too  small  or  too  large) .  The  "Pedestal  Under  Range"  shows  the 
number  of  signals  which  were  detected  but  too  small  to  measure.  "Ped 
Amp  Saturation  indicates  the  number  of  signals  which  produced  scatter 
magnitudes  large  enough  to  drive  the  system  amplifiers  into  nonlinear 
operation.  "Doppler  Under  Range"  shows  the  number  of  signals  with  Doppler 
signal  components  too  small  for  measurement.  The  next  block  of  data  is 
a  summary  of  the  statistical  parameters  for  the  measured  size  distribution. 
The  "Mean  Diameter"  is  the  geometric  mean  diameter.  "Sigma  G"  is  the 
geometric  standard  deviation  and  the  resulting  moments  are  computed  using 
the  generally  accepted  equations  from  statistics.  The  sixth  block  computes 
particle  number  density  and  aerosol  mass  concentration  using  the  model 
described  in  Appendix  II.  Inputs  required  are  relative  slit  width  Kappa, 
and  observed  average  N^.  The  material  density  is  computed  for  orthophos- 
phoric  acid  with  relative  humidities  an  input.  The  seventh  block  labeled 
Visibility  Calibration  Coefficient"  contains  the  exponents  and  coefficients 
for  a  least  squares  curve  fit  correction  for  input  data  which  has  been 
degraded  due  to  long  transmission  line  losses.  For  cable  lengths  less  than 
50m  between  processor  and  interferometer,  this  correction  has  not  been 
found  necessary.  The  eighth  block  is  the  computed  histogram  information 
from  the  number  of  measurements  indicated  in  the  size  accepted  statement. 

The  Bin  column  identifies  a  particular  particle  size  increment  with  a 
bin  number.  The  "Diameter"  column  lists  the  largest  particle  diameter 
associated  with  a  given  bin.  The  "Probability"  column  lists  the 
accumulative  probability  as  the  bin  number  increases.  "WTPOP"  lists  the 
weighted  population  associated  with  the  measurements.  The  weighted 
population  is  the  population  which  is  computed  to  exist  for  a  common  volume 
for  the  entire  particle  size  range.  It  is  computed  in  order  to  accommodate 
the  fact  that  the  sampling  volume  of  the  PSI  depends  on  particle  size. 

Hence,  large  particles  may  be  observed  over  a  larger  volume  than  the 
smaller  particles.  The  "Weight"  column  indicates  the  normalization  factors 
used  for  each  bin  increment  in  order  to  normalize  the  measured  population 
to  a  common  volume.  The  column  listed  "VIS"  indicates  the  cut  point  in 
the  visibility  function  for  a  particle  bin  as  computed  from  Equation  8. 
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The  column  listed  "POP"  is  the  actual  measured  particle  population. 

The  page  immediately  following  this  data  shows  the  probability  density 
distribution  and  is  intended  only  to  give  the  reader  some  idea  of  the 
shape  of  the  distribution.  The  amplitudes  of  the  probability  density 
plot  are  relative.  The  largest  amplitude  is  chosen  to  fit  a  convenient 
scale  on  the  paper.  The  first  column  indicates  the  largest  particle 
size  for  a  given  probability  density  amplitude.  The  second  column 
gives  the  accumulated  probability  density. 
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5.2  DEFINITIONS  OF  COMPUTED  VALUES 


ACCEPTANCE  RATIO  —  The  nunber  of  signals  accepted  for 
neasurenent  with  respect  to  the  nunber  of  signals  detected  which 
satisfy  the  counting  constraints  required  for  entry  into  the 
Histogran  Counter. 

ACQUISITION  TIME  —  The  tine  which  elapsed  during  data 
acquistion.  Does  not  include  data  reduction  tine. 

APERIODICITY  ACCEPTED  —  The  nunber  of  velocity 
neasurenents  which  passed  the  aperiodicity  test. 

APERIODICITY  FAILURES  —  The  nunber  of  velocity  readings 
which  failed  the  aperiodicity  test. 

ARITHMETIC  MEAN  —  The  arithnetic  nean  of  the  particle 
size  distribution. 


where  fi=weighted  fi/bin 

Di=largest  particle  dianeter 


in  a  bin. 

N  =total  weighted  bins 


BIH  —  The  bin  nunber  in  a  size  histogran. 

COEFFICIENT  —  A  coefficient  of  the  visibility  versus 
frequency  calibration. 

DIAH  —  Particle  dianeter.  The  units  are  nicron  (i.e., 
flicroneters)  if  the  fringe  Period  is  in  nicrons. 

DQPP  ---  The  neasured  nean  Doppler  conponent  of  the  signal. 

DOPPLER  UNDER  RANGE  —  The  nunber  of  neasurenents  for  which 
the  Doppler  conponent  of  the  signal  was  too  snail  to  neasure 
accurately. 

DR  —  The  integrator  range  (0,1, 2, 3)  for  the  Doppler  signal 
channel,  where  Range  0  =  Tine  constant  of  IE-6, Range  1  =  Tine 

constant  of  IE-5, Range  2  =  Tine  constant  of  IE-4, Range  3  =  Tine 

constant  of  IE-3. 


DIFF  1  —  Nuneric  frequency  of  occurance  per  unit  size 

inteval  in  percent. 
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EXP  —  The  exponent  of  frequency  in  a  visibility  versus 
frequency  calibration. 

FIRST  MOMENT  —  The  first  moment  of  the  particle  size 
distribution. 


V 


1 


l 

N 


n 

£ 

1=1 


fiDi 


where  fi=weighted  8/bin 

Di=largest  particle  diameter 
in  a  bin. 

N=total  weighted  bins 


FOURTH  MOMENT  —  The  fourth  moment  of  the  particle  size 
distribution. 


y 


4 


1 

N 


n 

£ 

1=1 


fiDi 


where  fi=weighted  ft/bin 

Di=largest  particle  diameter 
in  a  bin 

N=total  weighted  bins 


FREQUENCY  ---  Doppler  signal  frequency. 


GAMMA  —  A  dimensionless  constant  proportional  to  the 
Humber  Density.  It  is  a  function  of  the  average  number  of  cycles 
in  a  signal  and  the  fringe  period. 


GEOMETRIC  MEAN  —  The  Geometric  mean  is  defined  as: 


D 

g 

,i 

-  e*P<jJ 

n 

£ 

i=l 

f, InD . 
i  l 

HIGH 

PASS 

FAILURES 

— 

The  number  of  measurements  whose 

frequency 

was 

less  than 

the 

high  pass  filter  setting  in  the 

parameter 

list. 

KAPPA  —  An  experimental  constant  which  is  the  ratio  of 
the  x  dimension  to  the  y  dimension  in  the  probe  volume. 


K1*  MOMENT  —  u  =  -  T  f  n  k 

k  “  i  i 


LCC  - 

65535. 

Z  HASS 
inteval. 


-  The  long  clock  count.  An  integer  between  1  and 


- —  Hass 
fiDi3 

EfiDi3 


frequency  of  occurance  per  unit  size 

where  fi=veighted  H/bin 

Di=largest  particle 
diameter  in  a  bin 
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HASS  CONCENTRATION  —  Aerosol  Hass  contained  in  a  unit  voluhe  of 
aerosol . 

HASS  MEAN  —  Ratio  of  the  Fourth  to  the  Third  Monent  in 
the  particle  size  distribution. 

MATERIAL  DENSITY  —  Mass  per  unit  volune  of  a  naterial. 


MEAN  VELOCITY  —  The  nean  particle  velocity. 


NUMBER 

volune 

=  (■ 


DENSITY  — 

27rGamma*  3 
2 

NL  5 


Nunber  of  aerosol  particles  per  unit 

where  NL=long  count+precount+2 
5  r delta  *  fringe  period 


OVT  —  The  tiner  over-range  counter.  This  is  a  software 
counter  and  is  in  error  if  the  hardware  tiHer  overflowed  nore 
than  once  between  two  successive  data  readings. 


PED  The  neasured  nean  pedestal  coHponent  of  the  signal. 


PEDESTAL  UNDER  RANGE  -  The  nunber  of  readings  for  which 

the  pedestal  Tonponent  of  the  signal  was  too  snail  to  neasure 
accurately. 


POP  —  The  population  (i.e.,  the  nunber  of  counts)  in  a  bin 
of  the  particle  size  histogran. 

PR  —  The  integrator  range  (0,1 ,2,3)  for  the  pedestal 
signal  channel. 


PROBABILITY  ---  The  cunulative  probability  corresponding 
to  the  i^  consecutive  bin  beginning  fron  the  snallest  particle 

size  bin. 


RUN  NUMBER  —  The  nunber  used  by  the  operator  to  identify 
a  data  set. 


SCC  — -  The  short  dock,  count. 


SECOND  MOMENT 
distribution. 

1  n 
^2  =  N  Z 


i=l 


—  The  second  nonent  of  the  particle  size 

fiDi2 


SICHA  —  The  standard  deviation  of  a  neaurenent  (i.e.,  the 
square  root  of  the  variance). 
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SIGMA  D  —  The  standard  deviation  of  the  geonetric  wean  of 
the  particle  size  distribution. 

/  ,  f .ln2(D  /D  )  v  h v 

««-”>{(  »  s)  } 

SIZE  —  Particle  diameter  measured  in  the  units  in  which 
the  fringe  period  is  expressed. 

SIZE  ACCEPTED  —  The  nunber  of  measurements  which  passed 
all  tests  applied  to  the  size  data. 


THIRD  MOMENT  —  The  third  moment  of  the  particle  size 
distribution. 

1  n  3 

P3  "  N  fi°i 

TIME  —  Tine  in  seconds  Measured  from  that  tine  at  which 
the  operator  struck  the  space  bar  to  take  data. 

TIME  AV£RAG£D  MOMENTS  AND  MASS  CONCENTRATION  —  Tine 
averages  of  the  ensemble  of  the  Measurements  entered  into  the 
conputations. 

TIME  AVERAGED  HASS  CONCENTRATION  — 


6=1 


n 


Z  At 
6=1 


6 


TIME  AVERAGED  MOMENTS  — 


6-1 

n 

E  At 
6=1  B 


where  V>j.  =K  “  nonent 

n«=na55  concentration/trial 
At=acquisition  tine/trial 

TOTALIZED  HISTOGRAM  ---  Ensenble  average  of  all  histograns 
entered  into  conputations. 
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VEL  —  Measured  particle  velocity. 

VIS  ---  Measured  particle  visibility. 

VISIBILITY  CALIBRATION  CONSTANTS  —  The  coefficients  and 
corresponding  exponents  in  a  visibility  versus  frequency 
calibration. 

VOLUME/AREA  MEAN  —  Mean  particle  of  the  volune  to  area 
ratio. 


VOLUMETRIC  MEAN  —  Mean  particle  volune. 

-to" 

HEIGHT  —  The  weighting  factor  for  a  bin  in  the  particle 
size  histogran. 

UTPOP  —  The  weighted  population  (i.e.,  Measured  population 
tines  the  weighting  factor)  of  a  bin  in  the  particle  size 
histogran  acquistion.  Does  not  include  data  reduction  tine. 


*  NOTE:  Some  of  these  values  are  available  on  a  total  data 
printout  which  is  not  present  in  the  data  example 
summaries. 
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5.3  PSI  Data  Examples 


U  T  S  I  G  A  S 

D  I  A  G  N 

0 

S  T  I  C  S  DIVISION 

LONG  FRINGE  COUNT 

12 

SHORT  FRINGE  COUNT 

9 

# 

PRECOUNT 

o 

s-' 

4 

CLOCK  FREQUENCY 

70.00 

MHZ 

HIGH  PASS  FILTER- 

.  2000 

MHZ 

FRINGE  PERIOD 

6.000 

MAX I MUM  APER I ODI CITY 

1 . 0000 

7. 

DC  OFFSET 

4.  900 

VOLTS 

- 

€ 

DC  OFFSET 

4 . 900 

PED  ADJUST  OFFSET  42.0 

TIMER  RATE 

1000.0 

HZ  I RIG  CLOCK  ON 

SYSTEM  VOLTAGE  GAIN 

40.0 

HISTOGRAM  NOT  PRINTED 

SAMPLE  SIZE 

3000 

RUN  NUMBER 

1 

PARTICLE  SI 

ZING 

M 

ODE 

AUGUST  16,1980  DATA 

FROM  PS I 

#2 

ON  PAD 

TRIAL  #25  XM825  WP 

SMOKE  WEEK  III  EGLIN 

AFB 

PMT-600  SLIT- 100 

LPF-2.0MHZ 

ZCD-3.7 

DATA  READ  FROM  DISKETTE  FILE  TRIAL25  1 


MEAN  VELOCITY-  2.672  SIGMA-  .4670 


FREQUENCY  =  .4453  MHZ 

APERIODICITY  ACCEPTED  -  1320  -  44.0 

APERI GDI CITY  FAILURES  =  1680  =  56.0 

HIGH  PASS  FAILURES  =  0  0.00 

ACQUISITION  TTMF  =  95.14  - 

HOURS  =  21  MINUTES  =  46  SECONDS 


7. 

X 

7. 


17. 


SIZE  ACCEPTED  =  828 

VISIBILITY  TOO  SMALL  -  3 

VISIBILITY  TOO  BIG  =  43 

PEDESTAL  UNDER  RANGE  =  425 

DOPPLER  UNDER  RANGE  =  67 

GEOMETRIC  MEAN  =  .593 

FIRST  MOMENT  =  1.09 

THIRD  MOMENT  =  10.6 


=  27.6  7. 

=  .100  7. 

=  1.43  7. 

=  14.2  7. 

=  2 . 23  7. 

SIGMA  6  «  1.98 

SECOND  MOMENT  =  2.95 

FOURTH  MOMENT  -  42.0 


ARITHMETIC  MEAN  =  1.09 

VOLUME/AREA  MEAN  =  3.60 

VOLUMETRIC  MEAN  =  2.20 

MASS  MEAN  -  3.95 


7. 


GAMMA  =  1.54  ACCEPTANCE  RATIO  = 

NUMBER  DENSITY  (N/CC)  =  1 . 228E+05 
MASS  CONCENTRATION  (GM/CUBIC  METER)  = 
ALPHA  EXTINCTION  =  .6113  '  SIGMA 


.440  KAPPA  =  .360 

MATERIAL  DENSITY  =  1.4 

.991  WHEN  NL  =  1 

BAR  =  4.931 


l 


VISIBILITY  CALIBRATION  CONSTANTS 
EXP  -  0  •  COEFFICIENT  =  8.77333E-01 

EXP  -  1  '  COEFFICIENT  -  3.47290E-01 

EXP  =  2  COEFFICIENT  — 5.79147E-01 
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I o  co 


RUN  NUMBER®  1 

AUGUST  16,1980  DATA  FROM  PS I  #2  ON  PAD 
TRIAL  #25  XM825  WP 

SMOKE  WEEK  III  EGLIN  AFB 

PMT =600  SL I T= 1 00  LPF-2.0MHZ  ZCD»3.7 


NUMBER  DISTRIBUTION 


BIN 

DIAM 

DIFF  7.  PROBABILITY 

WTPOP 

WEIGHT 

VIS 

POP 

1 

.30 

.  646 

.  646 

2. 002E+03 

45.5 

.  997 

44 

2 

.60 

5.873E-02 

.  705 

182. 

45.  5 

.  988 

4 

3 

.90 

2. 033E-02 

.  725 

63. 0 

12.6 

.  972 

5 

4 

1 . 20 

1 . 569E-02 

.741 

48.  6 

6.  08 

.951 

8 

5 

1.50 

1 . 775E-02 

.  758 

55.  0 

3.93 

.925 

14 

6 

1.80 

1 . 355E-02 

.772 

42.  0 

2.80 

.  893 

15 

7 

2.  10 

1 . 897E-02 

.  791 

CO  © 

JO  .  o 

2.  10 

.  856 

28 

S 

2.40 

2.502E-02 

.816 

77.5 

1.65 

.815 

47 

.  9 

2.70 

2. 465E-02 

.  84 1 

76.4 

1.34 

.  770 

57 

10 

3.00 

3.  331E-02 

.  874 

103. 

1.11 

.  722 

93 

11 

3.  30 

2. 420E-02 

.  898 

75.0 

1 . 000 

.  670 

75 

12 

3.60 

2.  689E-02 

.  925 

oo  o 

VO.  J 

.  896 

.617 

93 

13 

3.90 

2. 804E-02 

.  953 

86.  9 

.790 

.  562 

110 

14 

4.20 

1 . 999E-02 

.  973 

62.0 

.  765 

.  506 

81 

15 

4.50 

1 . 047E-02 

.  984 

32.4 

.  676 

.449 

48 

16 

4 . 80 

6. 165E-03 

.  990 

19.  1 

.579 

.  393 

33 

17 

5.  10 

4. 654E-03 

.994 

14.4 

.577 

.  337 

25 

13 

5.40 

2, 74  J  F— 03 

.  997 

8.  50 

.531 

.  283 

1  6 

19 

5.70 

1 . 850E-03 

oejQ 

5.73 

.273 

.231 

21 

20 

6.00 

7. 305E-04 

1 . 000 

2.26 

ooo 

.  181 

8 

21 

6.30 

2. 836E-04 

1 . 000 

.879 

.  293 

0.00 

3 

MASS 

DISTRIBUTION 

BIN 

DIAM 

7.  MASS  PROBABILITY 

WTPOP 

WEIGHT 

VIS 

POP 

1 

.30 

1 . 639E-03 

1 . 639E-03 

2. 002E+03 

45.5 

.997 

44 

2 

.  60 

1 . 1 92E-03 

2. 83 IE— 03 

182. 

45.  5 

.  988 

4 

3 

.90 

1 . 393E-03 

4 . 223E-03 

63.0 

12.6  . 

.972 

5 

4 

1.20 

2. 548E-03 

6. 772E-03 

48.6 

6.  08 

.951 

8 

5 

1.50 

5. 630E-03 

1 . 240E-02 

55 . 0 

3.93 

.925 

14 

6 

1.80 

7. 427E-03 

1 . 983E-02 

42.  0 

2.80 

.  893 

15 

7 

2.  10 

1 . 651E-02 

3. 634E-02 

58 . 8 

2.  10 

.856 

28 

8 

2.40 

3. 250E-02 

6. 884E-02 

77.5 

1.65 

.815 

47 

9 

2.70 

4.558E-02 

.  1 14 

76.4 

1.34 

.770 

57 

10 

3.00 

8. 45 IE— 02 

.  199 

103. 

1.11 

.722 

93 

11 

3.30 

S. 172E-02 

.  281 

75.0 

1 . 000 

.  670 

75 

12 

3.60 

.  1  IS 

.  399 

83.3 

.896 

.617 

93 

13 

3.90 

.  156 

.  555 

86.9 

.790 

.562 

110 

14 

4.20 

.  1 39 

.  694 

62.  0 

.  765 

.  506 

81 

15 

4.50 

8. 965E-02 

.  784 

32.4 

.676 

.449 

48 

16 

4.  80 

6.407E-02 

.  848 

19.  1 

.579 

.  393 

33 

17 

5.  10 

5. 802E-02 

.  906 

14.4 

.577 

.  337 

25 

18 

5.  40 

.4 . 056E-02 

.  946 

8.  50 

.531 

.  283 

16 

19 

5.  70 

*3.  219E-02 

.  979 

5.73 

.273 

.231 

21 

20 

6.00 

1 . 483E-02 

.  993 

2.26 

OC'O 

.  181 

8 

21 

6.30 

6. 664E-03 

1 . 000 

.  879 

.  293 

0.00 

3 
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RUN  NUMBER*  1 

AUGUST  16,1980  DATA  FROM  PS  I  It  2  ON  PAD 
TRIAL  #25  XM825  WP 

SMOKE  WEEK  III  EGLIN  AFB 

PM T -600  SLIT- 100  LPF-2.0MHZ  ZCD-3.7 


N 


.  300 
,  600 


U  M 
646 
S73E 


E  R 


-02 


.  900 

2. 033E-02 

1 , 20 

1 . 569E-02 

1 , 50 

1 . 775E-02 

1 . 80 

1 . 355E-02 

2.  10 

1 . 897E-02 

2.40 

2. 502E-02 

2.70 

2. 465E-02 

3.  00 

3. 331E-02 

3.  30 

2. 420E-02 

3.60 

2. 689E-02 

3.90 

2. 804E-02 

4.20 

1 . 999E-02 

4.50 

1 . 047E-02 

4.80 

6. 165E-03 

5.  10 

4. 654E-03 

5.40 

2. 741E-03 

5.70 

1 . 850E-03 

6 . 00 

7.305E  -C’4 

6.30 

2.  836E--04 

I 

1 

1 

1 

I 

1 

( 

1 

MASS 

.  300 

1 . 639E-03 

.  600 

1 . 1 92E-03 

.  900 

1 . 393E-03 

1.20 

2. 548E-03 

1 . 50 

5. 630E-03 

1.80 

7. 427E-03 

2.  10 

1 . 65 1 E— 02 

2.40 

3. 250E-02 

2.70 

4. 558E-02 

3.00 

8. 451E-02 

3.30 

8. 172E-02 

3.  60 

.  118 

3.90 

.  156 

4.20 

.  139 

4.50 

S. 965E-02 

4.80 

6 . 407E-02 

5.  10 

5. 802E-02 

5.  40 

4. 056E-02 

5.  70 

3. 219E-02 

6.00 

1 . 483E-02 

6.  30 

6. 664E-03 

I  C  S  I  Z  E  DISTRIBUTI  0  N 

x* ************************************** *********** 

x#### 

X* 

X* 

X* 

X* 

X* 

X* 

X* 

X#* 

X* 

X** 

X** 

X* 

X 

X 

X 

X 

X 

\J 

4\ 

X. 

DISTRIBUTION 

X 

X 

X 

X 

X* 

X**- 

x#*-***##*-** 

X  *•*#*■  *■*■******•»•* 

X#**#*#*-*-******************* 

X#*#*##  *##*#*■****  ********** 

X************************************** 

X************************************************** 

X******************************************** 

X**************************** 

X******************** 

x****************** 

x************* 

x********** 

X**** 

X** 
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TOTALED  DISTRIBUTION  H  I  S  T  0  G  R  A  M 


RUN  NUMBERS  FOR  TOTALED  HISTOGRAM 
2  3  4  5  6 

AUGUST  16,1980  DATA  FROM  PS I  #2  ON  PAD 
TRIAL  ft 24  5“  ZUNI  WP 

SMOKE  WEEK  III  EGLIN  AFB 

PMT=600  SL I T= 1 00  LPF=2.0MHZ  ZCD=3.7 


TOTAL  SAMPLE  SIZE  =  6154 

APERIODICITY  ACCEPTED  =  3525 
TOTAL  ACQUISITION  TIME  =  226.0 

TOTAL  SIZE  ACCEPTED  =  2903 


GEOMETRIC  MEAN  = 

.  387 

SIGMA  G  = 

1 . 59 

FIRST  MOMENT  = 

.  555 

SECOND  MOMENT  = 

.912 

THIRD  MOMENT  = 

2.82 

FOURTH  MOMENT  = 

10.4 

ARITHMETIC  MEAN  = 

.  555 

VOLUME /AREA  MEAN  = 

3.  10 

VOLUMETRIC  MEAN  = 

1 . 41 

MASS  MEAN  = 

3.69 

.  360 
1.43 
18 


VISIBILITY  CAL I  BRAT I  ON  CONSTANTS 
EXP  =  0  COEFFICIENT  =  8.77333E-01 

EXP  =  1  COEFFICIENT  =  3.47290E-01 

EXP  =  2  COEFFICIENT  =-5. 79147E-01 


GAMMA  «  .937  ACCEPTANCE  RATIO  =  .573  KAPPA  = 

NUMBER  DENSITY  (M/SC)  -  2.777CK>4  MATERIAL  DENSITY  ** 
MASS  CONCENTRATION  (GM/CUBIC  METER)  =  5.S64E-02  WHEN  NL  = 
ALPHA  EXTINCTION  =  .6262  SIGMA  BAR  =  1.322 


**TIME  AVERAGED  STATISTICS** 

GEOMETRIC  MEAN  =  .389 

FIRST  MOMENT  =  .559 

THIRD  MOMENT  =  2.84 

ARITHMETIC  MEAN  =  '  .559 

VOLUME/AREA  MEAN  =  3.09 

VOLUMETRIC  MEAN  =  1.42 

MASS  MFAN  =  3.69 


SIGMA  G  -  1. 

SECOND  MOMENT  =  .9 

FOURTH  MOMENT  =  10.5 


MASS  CONCENTRATION  (GM/CUBIC  METER)  -  .108  WHEN  NL  =  18 

ALPHA  EXTINCTION  COEFF.  SIGMA  BAR 

.631135  1.33207 
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RUN  NUMBERS  FOR  TOTALED  HISTOGRAM 
2  3  4  5  6, 

AUGUST  16.1980  DATA  FROM  PS I  #2  ON  PAD 
TRIAL  424  5"  ZUNI  WP 

SMOKE  WEEK  III  EGLIN  AFB 
PMT=600  SLIT=1C»0  LPF=2.0MHZ  ZCD=3.7 


NUMBER  D I STR I BUT  I ON 


BIN 

DIAM 

DIFF  7. 

PROBABILITY 

WTPOP 

WEIGHT 

VIS 

POP 

1 

.30 

.  S44 

.  844 

2. 621E+04 

45.5 

.997 

576 

2 

.60 

3. 955E-02 

a  fc'C'**' 

1 . 228E+03 

45.5 

.  988 

27 

3 

.90 

1 . 663E-02 

.900 

517. 

12.6 

.972 

41 

4 

1.20 

1 . 155E-02 

.912 

359. 

6.08 

.  951 

59 

5 

1 . 50 

1 . 240E-02 

.924 

ooc; 

•_» 

3.  93 

.925 

98 

6 

1 .  SO 

8. 745E-03 

.  933 

272. 

2.80 

.  893 

97 

7 

2.  10 

8. 993E-03 

.942 

279. 

2.  10 

.  856 

133 

8 

» 

2.40 

9. 403E-03 

.  951 

292. 

1.65 

.815 

177 

9 

2.70 

8. 2S4E-03 

.  959 

257. 

1.34 

.  770 

192 

10 

3.  00 

7. 6 84 E- 03 

.  967 

239. 

1.11 

.  722 

215 

1 1 

3.30 

7. 856E-03 

.  975 

244. 

1 . 000 

.  670 

244 

12 

3.60 

7.  616E-03 

a  V  O  3 

237. 

.  896 

.617 

264 

13 

3.90 

6.  512E-03 

.  9S9 

202. 

.  790 

.562 

256 

14 

4.20 

4 . 508 E -03 

.  994 

140. 

.765 

.  506 

183 

15 

4.50 

3.  025E-03 

.  997 

94.0 

.676 

.449 

139 

16 

4.80 

1 . 622E-03 

.  998 

50.4 

.579 

.  393 

87 

17 

5.  10 

1 . 003E-03 

.  999 

31 . 2 

.577 

.  337 

54 

18 

5.  4U 

9h> 8t£-04 

1 . 000 

15.4 

.531 

OO'Tj 

29 

19 

5.70 

1 . 406E-04 

1-.  000 

4.37 

.273 

a  231 

16 

20 

6.00 

8.201E-05 

1 . 000 

2.55 

.  283 

.181 

9 

21 

6.30 

5. 660E-05 

1 . 000 

1.76 

.  293 

0.00 

6 

MASS 

DISTRIBUTION 

BIN 

DIAM 

7.  MASS  1 

PROBABILITY 

WTPOP 

WEIGHT 

VIS 

POP 

1 

.30 

8.  073E--03 

8. 07 6 E  08 

2. 621E+04 

45.5 

.  997 

576 

2 

.60 

3.02SE-03 

1 . 1 10E-02 

1 . 228E+03 

45.  5 

.988 

27 

3 

.90 

4. 297 E- 03 

1 . 540E-02 

517. 

12.6 

.972 

41 

4 

1.20 

7 . 072E-03 

2. 24 7 E— 02 

359. 

6.08 

.951 

59 

5 

1.50 

1 . 4S3E-02 

3. 730E-02 

cj  trj 

3.93 

.  925 

98 

6 

l.SO 

1 . 807E-02 

5. 537E-02 

272 . 

2.80 

.  893 

97 

7 

2.  10 

2. 95 IE— 02 

8. 488E-02 

279. 

2.  10 

.  856 

133 

8 

2.40 

4. 606E-02 

.  131 

292. 

1.65 

.815 

177 

9 

2.70 

5. 778E-02 

.  189 

257. 

1.34 

.770 

192 

10 

3.00 

7. 352E-02 

.  262 

239. 

1.11 

.722 

215 

11 

3.30 

.  100 

.  362 

244. 

1.000 

.  670 

244 

12 

3.60 

.  1 26 

.  488 

237. 

.896 

.617 

264 

13 

3.90 

.  137 

.  625 

202. 

.  790 

.  562 

256 

14 

4.20 

.  1  IS 

.743 

140. 

.765 

.  506 

183 

15 

4.50 

9. 769E-02 

.841 

94.0 

.  676 

.  449 

139 

16 

4.  SO 

6. 356E-02 

.905 

50.4 

.579 

•  393 

87 

17 

5.  10 

4. 716E-02 

.  952 

31 . 2 

.577 

00-7 

54 

18 

5.40 

'2. 767E-02 

.  979 

15.4 

.531 

.  283 

29 

19 

5.70 

9. 229E-03 

.  989 

4.37 

.  273 

.231 

16 

20 

6.00 

6. 277E-03 

.  995 

2. 55 

.  283 

.  181 

9 

21 

6.30 

5. 015E-03 

1 . 000 

1 . 76 

.  293 

0 . 00 

6 

56 

RIJN  NUMBERS  FOR  TOTALED  HISTOGRAM 
2  3  4  5  6 

AUGUST  16,1930  DATA  FROM  PS I  42  ON  PAD 
TRIAL  424  5"  ZUNI  WP 

SMOKE  WEEK  III  EGLIN  AFB 

PMT=600  SL I  T=  1 00  LPF=2.0MHZ  ZCD=3.7 


N  U  M  E  R 

I  C  SIZE  D  I  S  T  R  I  B  U  T  I  0  N 

.  300 

.844 

X*#**##***#**####*****#-*-*##  *■  *#*###  ##*#####•##  ##** 

.  600 

3. 955E-02 

Xbb 

.  900 

1 . 663E-02 

X 

1 . 20 

1 . 155E-02 

X 

1 . 50 

1 . 240E-02 

X 

1 .  SO 

8. 745E-03 

X 

2.  10 

8. 993E-03 

X 

2.40 

9. 403E-03 

X 

2.70 

8. 284E-03 

X 

3.  00 

7. 684E-03 

X 

3.30 

7 . 856E-03 

X 

3.  60 

7. 616E-03 

X 

3.  90 

6.512E-03 

X 

4.20 

4 . 508E-03 

X 

4.50 

3. 025E-03 

X 

4.80 

1 . 622E-03 

X 

5.  10 

1 . 003E-03 

X 

5.40 

4. 958E-04 

X 

5.70 

1 . 406F-04 

X 

6.00 

8. 20 1 E— 05 

X 

6.30 

5. 660E-05 

X 

.300 
.  600 
.900 
1.20 

1 . 50 
1.80 
2.  10 

2.40 

2.70 
3 . 00 
.3.30 
3.60 
3.90 
4.20 

4.50 

4.  SO 

5.  10 

5.40 

5.70 
6.00 
6.30 


MASS 
3. 073E-03 
3. 028E-03 
4. 297E-03 
7 . 072E-03 
1 . 433E-02 
1 . 307E-02 
2. 95 IE -02 
4. 606E-02 
5. 77SE-02 
7 . 352E-02 
.  1 00 
.  126 
.  137 
.  118 

9. 769E-02 
6 . 356E- 02 
4. 716E-02 
2. 767E-02 
9. 229E-03 
6. 27 7 E— 03 
5. 015E-03 


DISTRIBUTION 

X** 

X* 

X* 

Xbb 

Xbbbbb 

XBBBBBBBBBB 

Xbbbbbbbbbbbbbbbbbbbbb 

X  B  a-  B  -a-  B  B  *  B  *  B  B  B  B  B  B  B  B  *  B  B  B  B  B  B  B  B  *  B  *  B  *  B  B  B  B  B  B  B  B  B  B  B  B  *  *  B 
XBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

Xbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbb 

Xbbbbbbbbbbbbbbbbbbbbbbb 

Xbbbbbbbbbbbbbbbbb 

Xbbbbbbbbbb 

Xbbb 

X-BB 

XB 
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A  SIMPLE  MODEL  FOR  THE  VISIBILITY  OF  LARGE 
SPHERES  OBSERVED  WITH  A  LASER  VELOCIMETER 


W.  M.  Farmer 

The  University  of  Tennessee  Space  Institute 
Gas  Diagnostics  Research  Division 
Tullahoma ,  Tennessee 


ABSTRACT 

A  model  is  developed  to  express  the  signal  visibility  for  large 
refracting  or  reflecting  spheres  observed  with  a  fringe  type  laser 
velocimeter.  The  model  yields  a  closed  form  solution  for  the  visibility 
when  scattering  is  observed  off  axis  in  the  near  forward,  90°,  or  near 
backward  direction.  Results  show  that  particle  size  range  can  be  adjusted 
by  aperture  adjustments  at  the  scattered  light  receiver.  Comparison  of 
model  predictions  with  recently  published  numerical  computations  shows 
good  agreement. 

INTRODUCTION 

Considerable  effort  has  been  given  to  the  measurement  of  particle 

size  using  the  optical  system  of  differential  Doppler  laser  velocimeters 

(see  for  example  references  1-7).  Three  basic  approaches  to  this  problem 

can  be  identified:  1)  correlation  of  particle  velocity  with  aerodynamic 

particle  diameter  in  a  known  velocity  field, ^  2)  correlation  of  scattered 

2 

irradiance  with  particle  diameter,  and  3)  correlation  of  signal  visibility 

3-7 

or  contrast  with  particle  diameter.  Much  of  the  research  involving  the 


60 


latter  method  has  been  concerned  with  particle  diameters  comparable  to 
the  wavelength  illuminating  it.  An  approximate  model  which  relates  particle 
diameter  to  signal  visibility  appears  to  correlate  experimental  data  with 
theoretical  predictions  surprisingly  well  when  the  constraints  imposed  by 

g 

the  model  are  carefully  followed.  This  model  requires  that  the  scattered 
light  be  observed  with  a  large  receiver  aperture  centered  on  the  symmetry 
axis  of  the  transmitting  optical  system.  Measurement  of  spherical  particles 
with  this  model  requires  that  the  particle  diameter  be  less  than  an  illumi¬ 
nating  fringe  period  if  unambiguous  measurements  are  to  be  obtained.  Early 
research  with  particles  much  greater  than  a  wavelength  in  diameter  showed 

that  receiver  alignment  was  very  critical  if  agreement  with  the  approximate 

4 

model  was  to  be  obtained.  It  was  also  found  that  when  the  particle  diameter 
was  greater  than  a  fringe  period  and  the  scattered  light  observed  at  some 
angle  off  the  axis  of  symmetry,  the  measured  visibility  was  usually  much 
higher  than  predicted  by  the  approximate  model.  Thus  it  was  concluded  that 

2 

a  visibility  measurement  could  not  be  used  for  particle  diameter  estimates. 

However  recently  published  calculations  have  shown  this  conclusion  to  be 
9 

incorrect.  Bachalo  has  shown  that  for  particles  much  greater  than  a 
wavelength  in  diameter  and  at  viewing  angles  where  the  reflected  or  re¬ 
fracted  component  of  the  scattered  light  is  greater  than  the  diffracted 
component,  particle  diameter  can  be  related  to  a  visibility  function  which 
is  monotonic  over  a  significant  size  range.  This  work  shows  that  the 
visibility  function  is  monotonic  when  particle  diameter  is  greater  than 
a  fringe  period  and  that  the  visibility  function  depends  on  the  receiving 
optics  F  number,  particle  index-of-refraction,  and  viewing  angle.  The 
analysis  performed  by  Bachalo  consisted  of  numerical  integrations  of  the 
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classical  electromagnetic  field  equations  describing  refraction  and  re¬ 
flection  by  homogeneous  spheres  illuminated  by  plane  waves.  As  such,  this 
work  is  acceptably  rigorous.  However,  because  of  the  approach  taken  and 
the  numerical  integration, it  is  difficult  to  discern  or  physically  visualize 
the  scaling  laws  which  relate  visibility  to  particle  index-of-ref raction , 
receiver  F  number  or  viewing  angle.  The  purpose  of  this  paper  is  to 
develop  a  simple  model  which  elucidates  the  results  obtained  by  Bachalo 
and  yields  simple  scaling  laws  which  can  be  used  to  configure  optical 
systems  for  measuring  a  broad  range  of  particle  sizes  for  a  fixed  fringe 
period. 


SCATTERED  IRRAD LANCE  APPROXIMATION  FOR  REFRACTING  OR  REFLECTING  SPHERES 
ILLUMINATED  BY  TWO  COHERENT  BEAMS. 

A  classical  approximation  for  the  electric  fields  scattered  by  a 
large  (relative  to  the  wavelength)  homogenous  sphere  involves  separating 
the  scattered  field  E  ,  into  three  contributions 


ES  =  ED  +  Et  +  ER 


Cl 


where  E^  is  the  field  component  due  to  diffraction,  E^_  is  that  due  to 

refraction  and  E  is  that  due  to  reflection.  Hodkinson’s  work  can  be 
K 

manipulated  to  show  that  the  mean  value  of  E  will  be  the  predominant 
field  component  in  the  forward  scatter  direction  (rather  than  diffraction) 
when  the  scatter  angle,  3  is  such  that: 
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(2 


*  1/3 


m  +  1 


/ 


m 


where  X  is  the  illuminating  wavelength,  D  is  the  particle  diameter  and 
m  is  the  index  of  refraction  (assumed  real) .  For  the  near  axis  backscatter 
the  reflected  field  is  predominant  when 


(3 


3  «Y 


where  y  is  the  rainbow  angles  and  nth  order  refraction  in  the  back  direction 
is  ignored. 

For  this  analysis  the  geometry  shown  in  figure  1  will  be  assumed. 

The  centerlines  of  two  coherent  beams  cross  at  the  center  of  an  XYZ  co¬ 
ordinate  system.  The  beam  centerlines  are  separated  by  an  angle  a. 

The  Z  axis  is  parallel  to  the  bisector  between  the  beams.  Y  is  perpendicular 
to  Z  and  lies  in  the  plane  of  the  beam  centerlines  and  X  is  perpendicular 
to  YZ.  A  scattered  light  collection  aperture  is  centered  at  some  viewing 
angle  3  and  subtends  a  collection  angle  defined  by  the  F  number  of  the 
receiver.  Into  this  coordinate  system  a  large  homogeneous  sphere  of  radius 
r^  is  placed  at  some  position  X  =  0,  Y  and  Z  =  0  as  shown  in  figure  2. 

The  sphere  in  the  limit  D>>X  may  be  imagined  to  function  as  a  spherical 
lens.  Those  rays  having  angles  of  incidence  near  zero  will  be  predominantly 
refracted  by  two  refractions  (ray  entrance  and  ray  exit) .  These  rays  will 
appear  to  focus  at  some  distance,  f,  from  the  center  of  the  sphere.  For 
angles  of  incidence  less  than  about  20°,  f  is  nearly  independent  of  the 
angle  of  incidence.  Differences  in  reflection  due  to  different  values 
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of  the  Fresnel  reflection  coefficients  are  also  minimal.  Rays  contained 
within  this  angle  of  incidence  will  all  appear  to  diverge  from  f  and  form 
a  spherical  wave  with  radius  of  curvature  measured  from  f .  As  the  angle 
of  incidence  increases  f  is  reduced  and  spherical  aberration  becomes  signi¬ 
ficant.  Furthermore  rays  undergoing  more  than  two  refractions  may  be  ob¬ 
served.  It  will  be  assumed  that  the  refracted  light  observations  are  made 
with  an  observation  aperture  and  observation  angle  which  is  small  enough  to 
preclude  these  effects  thus  allowing  f  to  be  uniquely  specified.  Experi¬ 
mentally  these  assumptions  are  found  not  to  be  too  restrictive.  Consider 
figure  3a  which  is  a  photograph  of  a  114  micrometre  diameter  crown  glass 
(index  of  ref raction  . 5)  sphere  illuminated  by  two  coherent  beams.  The 

central  portion  of  the  sphere  images  the  interference  fringes  produced  by 
the  incident  beams.  A  large  portion  of  the  sphere  cross-section  is  shadow¬ 
ed  indicating  that  rays  refracted  by  this  portion  of  the  sphere  did  not 
pass  through  the  F/2  camera  lens  used  to  obtain  the  photograph.  If  it 
is  assumed  that  rays  reaching  the  photograph  are  only  twice  refracted  by 
the  sphere,  then  a  ray  trace  for  this  geometry  will  show  that  the  angle 
of  incidence  could  not  have  been  more  than  about  21  ,  which  is  sufficiently 
small  for  the  approximations  to  be  used  here.  With  these  assumptions  it 
is  clear  that  if  two  angularly  separated  plane  waves  illuminate  the  sphere, 
two  focus  points  will  be  obtained.  These  foci  produce  two  spherical  waves 
which  produce  nonlocalized  interference  fringes.  The  fringe  contours  are 
identical  in  shape  to  those  observed  in  Young’s  experiment  using  two  cir¬ 
cular  point  sources.  As  is  well  known,  the  fringe  contours  are  hyperboles 
and  are  beautifully  illustrated  in  figure  5  of  Bachalo’s  work.  Figure  3b 
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shows  the  resulting  fringe  contours  for  the  sphere  shown  in  figure  3a.  The 
model  to  be  used  in  this  analysis  contains  assumptions  radically  different 
from  those  made  originally  for  a  visibility  analysis  where  diffraction  is 


3 

predominant.  In  that  analysis  it  was  assumed  that  the  spherical  waves 


scattered  by  the  particle  illuminated  by  two  beams  have  a  common  center. 

In  the  case  where  scatter  is  predominantly  from  the  edge  of  the  particle 
and  the  angle  between  the  illuminating  beams  is  small  this  is  a  reasonable 
assumption.  When  refraction  or  reflection  is  the  dominant  mechanism  through 
which  light  reaches  the  collecting  aperture,  this  model  shows  that  such  an 
assumption  is  no  longer  valid. 

With  the  assumption  that  only  twice  refracted  light  reaches  the  col¬ 
lecting  aperture,  the  electric  fields  originating  from  the  refracted  foci 
may  be  written  by  assuming  spherical  waves  and  the  geometry  shown  in 
figure  2.  These  fields  are 


(5 


where  the  E  (i=l,2)  are  the  incident  field  amplitudes  for  beams  1  and  2. 
oi 


K  are  the  wave  vector  magnitudes 


(6 
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and  (0^  are  the  Doppler  shifted  optical  frequencies  given  by: 


Wtl  =  Wo  +  l  •  <5tl“  *01> 


(7 


Wt2  Wo  +  ~  (~t2  ~02) 


(8 


V  is  the  particle  velocity.  K_^  is  the  incident  wave  vector  for  the  ith 
beam.  It  will  be  assumed  that  the  receiver  aperture  is  small  enough  to 
approximate  K  as  K^.  The  geometry  in  figure  2  shows  that  r^  and  r^  can 
be  written  as 


r2  =  r  2  +  (X  +  X')2  +  (Y  +  Y'  +  A/2)2 
loo  o 


(9 


r  2  =  r  2  +  (X  +  X')2  +  (Y  +  Y'  -  A/2)2 
ZOO  o 


(10 


where  the  Xf ,  YT  coordinate  system  is  centered  in  the  scattered  light 
collection  aperture.  The  coordinates  of  a  point  in  a  circular  aperture 
coordinate  system  are  given  by 
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X'  =  psine 


(11 


Y'  =  pcos9 


(12 


The  total  field  at  the  aperture  point  defined  by  r^  and  Is 


E. 


Etl  +  Et2 


(13 


The  irradiance  resulting  from  Et  is  approximately  given  by 


(14 


where  Z'  is  the  impedance  of  the  medium  in  which  the  observations  are 
made  and  c  is  the  speed  of  light  in  vacuum.  Using  equation  13  in  equation 
14  results  in  the  well  known  form  for  the  scattered  irradiance  from  two 
coherent  beams. 
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lt  -  8ir  '{lEtl!2+  lEt2!2  +  2Re(EtlEt*2» 


(15 


2  ^ 

For  simplicity  assume  that  E^  =  E^  =  Eq.  | E ti  I  and  |  Et2 1  are 

thus  constants  and 


(16 


(17 


The  last  term  in  brackets  in  equation  15  is  the  so  called  interference 
term.  Use  of  equations  4,  5,  7  and  8  allows  the  interference  term  to  be 
written  as 


2Re(EtiEt2)  - 


2  I E 

1  r 


rlr2 


cos(Kt(r1-r2)  +  V 


(<5t2-  Sti 


(18 


+  ^~02  ~01^ 


to  obtain  r. -  r2  use  is  made  of  equations  9  and  10  to  compute 
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(19 


r 


2 

1 


2 


2 (Y  +  pcos9)A 
o 


for  the  case  when  r2»A  and  r2  =  r  [l  +  %  tan2$| 

equation  19  reduces  to 


r 


1 


A(Yq+  pcosB) 
ro  D*  +  ^  tan^jT] 


(20 


2  2  2  2 

where  tan  3  =  (X  +  Y  )/r  arises  from  the  second  order  terms  in  the 
o  o  o 

Taylor’s  series  expansion  of  r^  and  r2«  It  is  important  to  note  that  since 
spherical  waves  are  assumed  to  originate  from  point  sources  in  the  X,Y  co-ordin¬ 
ate  system  any  results  involving  r  ^ ,  r^are  invariant  to  simple  rotation  about 
the  Y  axis  (3>0,  X  =0) .  This  means  that  if  Xf,Y!  is  kept  perpendicular  to  r^ 

and  the  origin  of  x’,Yf  kept  in  the  X  ,Yq  plane,  then  the  only  angular  dependence 

2 

on  3  should  arise  from  a  lengthening  of  r^  as  3  increases  •  The  (1-f^tan  g)  factor 
accounts  for  this  increase.  Hence,  no  angular  correction  is  necessary  to  account 
for  aperture  ellipticity  as  might  be  expected  if  the  X* ,Yf  plane  were  held  in  the 
Xq,Yo  plane.  It  is  straightforward  to  show  that  the  terms  in  equation  18  resulting 
from  the  Doppler  shift  can  be  written  as 
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(21 


(*tl 


-  5t2>  = 


V  K  A/R 

y  t 


Y  *<K02  -  K01)  =  VyKo2sin(a/2)  (22 

R  is  the  distance  from  the  particle  centered  coordinate  system  to  a  point 
of  interest  in  the  aperture.  In  previous  calculations  the  term  expressed 
by  equation  21  was  neglected  on  the  basis  that  V  •  (K  was 

proportional  to  AA/A  where  AX  is  the  difference  in  optical 

wavelengths  resulting  from  the  Doppler  shift.  AX/X  is  of  the  order 

of  10  ^  or  less  in  nearly  all  cases  of  interest  and  hence  can  be  reasonably 
neglected.  In  most  cases  A/R  is  also  small  although  cases  can  be  en¬ 
visioned  where  it  might  be  significant.  This  analysis  will  assume  that 
this  term  is  negligible.  Equations  20  and  22  result  in  the  interference 
term  being  written  as 


2Ee(EtlEt2  > 


A(Y  +  pcos0) 

- 2 - —  +  2ttY/S) 

r  [I  +  tan  0 


(23 


6  =  X/2sin(a/2) 


(24 


Y  =  V  t 

y 


(25 
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6  is  the  fringe  period  generated  by  the  beams  striking  the  particle. 


IRRADIANCE  COLLECTED  BY  A  CIRCULAR  APERTURE 

The  photodetector  used  to  detect  light  scattered  by  the  particle 
is  assumed  to  respond  to  the  total  irradiance,  I^,  incident  on  the  scat¬ 
tered  light  collection  aperture.  For  a  circular  collection  aperture, 

IR  is  given  by 


I  -  fa  /2Tr  I  pdpd0 
r  o  o  t 


(26 


Using  equations  17  and  23  in  equation  26  results  in  an  integral  with  a 
closed  form  solution  given  by 


ca2I  [l  +  2J,(K  Aa/r')  cos  (2ttY/6+KAY  /r ' )1 
_ O  L  1  o  o  o  o  J 

4Z’r2  K  Aa/r' 

o  o  o 


(27 


where  =  rQ  [l  +  ^tan2^]  and  ls  a  first  order  Bessel 

function  of  the  first  kind. 
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VISIBILITY  FUNCTION  FOR  A  CIRCULAR  APERTURE 


It  is  straightforward  to  show  that  the  relative  amplitude  of  the 
interference  term  in  equation  27  is  identical  to  the  signal  visibility 
function,  V,  as  defined  in  previous  research.  In  this  case  equation  27 
shows  that  V  is  given  by 


V  =  2J,  (KAa/r ? )/KAa/r 
l  o  o 


(28 


The  visibility  function  in  equation  28  is  identical  in  form  to  that  cal¬ 
culated  previously  for  circular  diffracting  particles  and  a  sufficiently 
large  aperture  centered  on  the  bisector  between  the  incident  beams.  However, 
comparison  shows  there  are  important  differences  between  the  two  cases. 
Equation  28  shows  that  in  the  case  of  refraction  the  visibility  depends 
directly  on  the  angular  extent  a/r^  of  the  collection  aperture,  its  angular 
position  and  the  separation,  A  ,  of  the  two  point  sources  generated  by 
incident  beams  focused  by  the  particle.  As  will  be  shown,  A  depends  on  both 
the  diameter  of  the  particle  and  its  index-of-ref raction.  It  is  interesting 
to  note  that  a  phase  term  K  yQ/r0  is  predicted  for  the  Doppler  signal 
portion  of  1^.  This  phase  term  depends  on  the  viewing  angle,  y Q/r0  in  the 
plane  of  the  incident  beams,  the  particle  size,  and  index-of-ref raction . 

These  dependences  have  not  been  calculated  by  the  more  rigorous  solution 
to  the  problem  posed  by  Bachalo.  (There  was  no  apparent  need  to  do  so.) 
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Similar  effects  using  Mie  scattering  theory  have  also  been  predicted  for 
scattered  light  observations  involving  particles  comparable  to  a  wavelength 
in  diameter. 

VISIBILITY  FOR  NEAR-AXIS  FORWARD  SCATTER 

The  geometry  shown  in  figure  4  can  be  used  to  express  A  as  a  function 
of  beam  intersection  angle,  a  ,  and  the  focal  lengths  developed  by  the 
sphere  for  the  individual  beams.  Inspection  of  figure  4  shows  that 

A/2f  =  sin (a/ 2)  (29 

A  ■  f(2sina/2)  (30 

f  may  be  obtained  by  simple  ray  traces  through  the  sphere  with  an  applica¬ 
tion  of  Snell's  law  assuming  that  the  angle  of  incidence  0i«  1*  Appendix 
1,  using  this  approach,  shows  that  f  can  be  written  as 

f  _  Pm 

4(m-l)  (31 

where  m  is  the  particle  index  of  refraction,  D  is  its  diameter,  and  it  is 
assumed  that  the  refractive  index  of  the  medium  surrounding  the  particle 
is  1. 
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The  F  number  of  the  scattered  light  collection  aperture  is  defined  as 


F  =  r  /2a 
o 


(32 


Equations  24,  30,  31,  and  32  may  be  used  to  express  the  argument  of  the 


visibility  function  in  28  as 


K  Aa/r’ 
o  o 


tt(D/ 6)m 


(33 


r% 

4F[l  +  ^tan  3]  (m  -  1) 


Thus  the  receiver  aperture  and  index  of  Tefxacrtion  of  the  particle  scale 
the  range  over  which  a  particle  size  measurement  can  be  obtained  from  a 
visibility  measurement.  This  result  agrees  with  conclusions  drawn  by  Bachalo 
from  his  numerical  integration  studies  of  the  problem.  Resolution  of  the  visi¬ 
bility  measuring  electronics  limits  the  smallest  measurable  value  of  equation 
33  to  about  0.1  tt  ,  The  visibility  is  a  mono  tonic  function  of  its  argument 
for  values  up  to  about  1,0  TT  .  Since  the  D/6  (ultimately  D  is  measured  in 
terms  of  a  known  6  )  is  usually  the  parameter  to  be  determined,  the  size 

range  limits  for  a  given  value  of  F  and  m  are 


(D/6)  .  -  0.4F  fo"-1—  [l  +  h  tan2$] 

min  m 


(34 
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(35 


These  equations  show  that  with  a  fixed  fringe  period,  6,  and  particles  all 
having  the  same  index-of-ref raction,  the  size  range  covered  is  easily  varied 
by  adjusting  the  F  number  of  the  receiving  aperture. 

These  results  suggest  that  when  D/6  is  known,  the  index  of  refraction 
of  the  particle  can  be  measured.  For  example,  suppose  liquid  droplets  of 
a  known  size  were  produced  with  a  monodisperse  generator  and  the  measured 
visibility  corresponded  to  an  argument  r|.  Solution  of  equation  33  for  the 
index-of-refraction  yields 


4Fr) 


(36 


m 


4Fr|  -  ttD/6 


Equations  30  -  32  shows  that  the  phase  shift,  ^  ,  in  the  signal  for  this 
case  is  given  by 


=  KAy  /r 
J  o  o 


(37 


*  *  2*  warn  6 


(38 
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The  relative  phase  of  the  signal  is  seen  to  depend  on  the  aperture  position 
3^  only  and  not  on  the  aperture  size.  It  is  linearly  proportional  to  D/6. 
Figure  5  shows  the  phase  shift  effect  as  D/6  increases  from  1  to  10  for 
an  F/10  aperture  located  at  3^=  0.1  and  an  index  of  refraction  of  4/3.  This 
result  suggests  that  particle  size  could  be  obtained  by  measuring  the  phase 
shift,  ip,  rather  than  the  visibility  V. 

VISIBILITY  FOR  NEAR  AXIS  BACKSCATTER 

Previous  numerical  studies  by  Bachalo  suggested  that  this  approach 

could  also  be  used  to  measure  D/6  from  off  axis  backscatter  measurements 

of  highly  reflecting  or  absorbing  particles.  Assume  that  the  observation 

angle  is  far  removed  from  any  glory  or  rainbow  angles.  The  previous  results 

can  be  modified  in  a  straightforward  fashion  to  show  this.  In  this  case 

the  particle  is  no  longer  assumed  to  be  a  refractive  lens  but  a  spherical 

mirror.  The  focal  length  of  a  spherical  mirror  is  half  the  radius  of 

curvature.  Equation  30  then  shows  that  A  for  backscatter,  A  ,  is  given 

B 

by 


D 


4 


(2sin (ot/2) ) 


(40 


The  argument  for  the  visibility  in  the  backscatter  case  may  thus  be 
written  as: 


K^Aa/r^  =  7r(D/6)/4F[l  +  %  tan^gj 


(41 


The  phase  shift  in  the  Doppler  signal  for  backscatter,  iJj  can  be  computed 

D 

from  equation  37  and  40  as 


=  2tt (D/5)3  /A 

&  y 


(42 


Equation  41  shows  that  a  much  greater  range  in  D/6  can  be  covered  in 
backscatter  measurements  since  the  size  range  is  independent  of  the  index 
of  refraction.  For  example  when  near  axis  forward  and  backscatter  measure¬ 
ments  are  compared  for  water  drops  (m  =  4/3),  the  maximum  value  of  D/6  which 
can  be  measured  in  the  forward  scatter  direction  is  F  while  in  the  back¬ 
scatter  direction  it  is  4F.  This  result  suggests  that  backscatter  measure¬ 
ments  over  the  same  size  range  as  forward  scatter  may  give  accurate  results 
even  for  highly  refracting  particles  since  the  collection  aperture  may  be 
significantly  reduced  thereby  eliminating  rays  resulting  from  multiple 
refractions  toward  the  backward  direction. 

VISIBILITY  FOR  90°  SIDE  SCATTER 

The  simple  model  used  for  analysis  of  near-axis  forward  and  backscatter 
measurements  may  also  be  used  to  develop  an  approximate  expression  for  the 
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visibility  function  for  light  scattered  90°  relative  to  the  bisector  between 
the  incident  beams.  Light  reaching  the  collection  aperture  in  this  case 
results  predominantly  from  reflection.  For  observations  made  at  90°  with 
respect  to  the  bisector  between  the  beams,  reflected  rays  originate  for  angles 
of  incidence  of  -  a/2  and  —  +  a/2  for  beams  separated  by  an  angle  a  . 

The  lateral  separation  of  the  reflected  light  sources  on  the  sphere  is  thus 
given  by 


Ago  -  Y  (cos(tt/4  -  a/2)  -  cos(tt/4  +  a/2)) 


(43 


Ago  ~  Dsin (tt/4) sina/2 


(44 


Fig.  6  illustrates  the  coordinate  geometry  for  the  estimation  of  AgQ» 

The  collection  aperture  again  can  be  modeled  as  one  observing  two  point  sources 
separated  by  a  distance  Ago  .  The  visibility  function  for  a  circular  aperture 
is  thus  given  by  equation  28.  The  argument  of  the  visibility  function  in 
this  case  is  thus  found  from  equation  44  to  be 


kAa/r 

o 


/2_ 

4F 


ttD/6 


(45 


The  phase  shift  observed  in  the  Doppler  signal  for  90°  side  scatter  is 

found  using  equations  37  and  44  to  be 
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(46 


(D/6)  6y 


These  results  apply  to  the  experimental  case  of  the  highest  possible  spatial 
resolution  which  can  be  achieved  with  an  instrument  of  this  type. 

THE  GENERALIZED  VISIBILITY  FUNCTION 

The  previous  results  for  a  circular  aperture  may  be  generalized  to 

apertures  of  virtually  any  shape.  The  approach  is  similar  to  that  followed 

by  Robinson  and  Chu  in  analyzing  the  visibility  function  for  diffraction 
5 

only.  With  assumptions  made  previously,  the  intensity  distribution  in  the 
plane  of  the  collection  aperture  may  be  written  using  equations  15,  and 
18  -  22  as 


I  =  I  (1  +  cos((kAY '/r’)  +  0)) 
to  o 


(47 


0  =  kAY  /r'  +  2ttY/6 
o  o 


(48 


where  I  includes  all  necessary  constants  and  it  is  assumed  that  the 
waves  incident  on  the  particle  are  equal  in  amplitude.  The  irradiance 
collected  by  the  receiving  aperture  is  given  by 


r 

kA  Y T  \ 

I  =  I 

dA  +  I 

COB 

— T“  +  0 

R  o  J 

P  °  J 

r T  1 

A 

A 

°  / 

(49 
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where  is  the  area  of  the  receiving  lens.  It  is  convenient  to  define  a 
pupil  function  P(X,Y)  such  that 


P(X,YI)  =  1  for  all  x,  y  falling  inside  A 


0  otherwise 


can  be  written  as 
K 


Xo 


t 

|kA  Y' 

p(x;yj)dA  +i 

P(X',Y')cos  - ^  +0 

J  oo  P  OJ 

oo  \  ro 

dA 


(50 


The  integrals  in  equation  50  may  "be  recognized  as  the  Fourier  Transforms 

of  the  pupil  function  evaluated  at  spatial  frequencies  f  ,  and  f  ,  where 

xo  yo 


f  =  kA/r'  (a 

° 


(51 


f  =0 
x 

o 


(b 


with  this  identification,  equation  50  may  be  written  as: 


P(0,f  ) 

IR  =  IoP(0,0)  [l  +  - - 2-  cos(0)] 


P(0,0) 


(52 


where  P(f^,f^)  is  understood  to  represent  the  Fourier  transform  of  the 
pupil  function  P(X,Y).  From  equation  52  it  is  clear  that  the  visibility 
function  is  given  by 
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V  =  P(0,f  )/P(0,0) 


(53 


Equation  53  shows  how  the  visibility  function  varies  as  a  function  collection 
aperture  geometry.  For  example,  a  rectangular  aperture  of  width  2L  and 
uarbitrary  length"  oriented  with  sides  parallel  to  the  interference  fringes 
in  the  reception  plane  requires  a  visibility  function  given  by 


sin(kAL/r  ) 

v  =  - —  (54 

kAL/r 

o 

The  arguments  of  the  visibility  functions  for  spherical  particles  developed 
for  circular  apertures  transpose  directly  with  the  understanding  that  the 
F  number  of  the  receiver  is  now  defined  as  rQ/2L.  Equation  53  also  indicates 
the  possibility  of  aperture  shaping  to  achieve  increased  resolution  at  the 
small  or  large  ends  of  the  particle  size  range.  For  example,  a  square 
aperture  oriented  with  its  diagonal  perpendicular  to  the  planes  defined  by 
the  sample  space  interference  fringes  requires  a  visibility  function  given 
by 

sin^(kAL  /?/2r  ) 

V  =  - -  (55 

(kAL  /I  /r  ) 2 
o 

In  this  case  the  large  size  limits  of  the  visibility  function  are  extended 
over  those  required  by  the  one  dimensional  aperture  described  by  equation  54 
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because  the  visibility  is  monotonic  over  a  larger  range  (the  first  zero 
in  the  visibility  is  extended  by  a  factor  of  /2  ) .  The  small  size  limit 
is  also  increased  but  because  the  visibility  is  squared  the  small  size 
resolution  limit  can  be  decreased  from  about  0.1  TT  in  the  limiting  argument 
of  the  visibility  to  about  0.075  TT  to  obtain  the  same  limiting  value  of 
visibility.  This  represents  an  approximate  reduction  in  the  smallest 
resolvable  particle  size  of  about  25%. 


COMPARISON  WITH  NUMERICAL  COMPUTATIONS 

Bachalo's  approach  to  the  solution  of  this  problem  is  to  use  more 
exact  expressions  for  the  refracted  waves  and  to  numerically  integrate  the 
resulting  expressions  to  obtain  the  visibility.  As  such  this  method  should 
yield  more  precise  results  than  those  of  the  approximate  model  developed 
here.  Thus,  it  is  of  interest  to  compare  results  from  the  two  approaches 
to  determine  how  well  the  approximate  model  compares  with  the  numerical 
calculations . 

Table  I  summarizes  the  results  of  the  approximate  solution  for  circular 
apertures  and  near  forward,  90°,  and  near  backseat ter  estimates.  When  the 
approximate  model  results  are  compared  with  Bachalo’s  calculations  it  is 
immediately  seen  that  both  agree  in  predicting  that  the  same  functional 
form  will  result  as  the  particle  index  of  refraction  or  aperture  size  is 
varied.  A  reasonable  point  of  comparison  for  the  two  results  is  the  pre¬ 
diction  of  the  value  of  D/6  for  the  first  zero  in  the  visibility  function. 
Table  II  shows  the  results  of  comparing  the  two  predictions  for  the  index- 
of-refrac tion  and  F  number  values  used  in  Bachalo’s  numerical  calculations. 
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In  view  of  the  assumptions  made  for  the  approximate  model  evaluation, 
the  agreement  is  quite  good.  The  comparison  shows  the  essential  correct¬ 
ness  of  the  scaling  relationship  effects  on  the  visibility  resulting  from 
changes  in  viewing  angle,  collection  aperture,  or  index-of-ref raction. 

Note  also  that  as  the  index-of-ref raction  increases,  the  approximation 
becomes  better, 

SUMMARY  AND  CONCLUSIONS 

A  simple  model  has  been  developed  which  gives  in  closed  form  solution 
a  prediction  of  the  visibility  dependence  on  particle  index-of-ref raction, 
receiver  F  number,  and  viewing  angle  for  large  refracting  particles.  The 
analysis  also  shows  the  effects  of  these  parameters  on  the  phase  of  the 
Doppler  signal  and  how  receiving  apertures  can  be  shaped  to  obtain  visi¬ 
bility  functions  which  are  consistent  with  a  desired  size  range  for  a  given 
fringe  period.  Comparison  of  the  model  results  with  recently  published 
results  obtained  by  numerical  integration  shows  reasonable  agreement. 

The  ability  to  adjust  the  measurable  size  range  by  simply  adjusting 
an  iris  In  front  of  the  receiving  lens  is  a  significant  operational  simpli¬ 
fication  over  past  paraxial  interferometers  which  were  designed  to  produce 
very  large  fringe  periods  for  large  particle  measurements. 

Results  in  this  analysis  help  to  qualitatively  explain  the  alignment 
sensitivity  problem  described  in  past  work  for  the  measurement  of  large 

g 

refracting  particles  paraxially.  When  optical  alignment  is  such  that 
refractive  amplitudes  are  greater  than  diffracted  amplitudes,  the  visibility 
can  be  expected  to  increase.  Hence,  paraxial  observations  of  particles 
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producing  diffracted  components  which  are  effectively  blocked  by  beam  stops 
can  lead  to  visibilities  of  the  magnitude  predicted  here. 

The  method  developed  here  may  also  be  used  for  predicting  visibility 
functions  for  non-spherical  particles.  The  parameter  measured  by  the 
visibility  function  is  A.  Hence,  as  the  particle  shape  becomes  non- 
spherical,  geometric  ray  traces  can  be  used  to  relate  6  to  a  particle  size 
parameter. 
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COMPARISON  OF  APPROXIMATE  AND  NUMERICAL  PREDICTIONS  FOR  D/6 
VALUES  AT  FIRST  ZERO  IN  VISIBILITY  FUNCTION 
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(1)  Private  communication  with  Bachalo  indicates  that  the  published  values  are  incorrect. 
They  are  not  consistent  as  reported  in  the  publication  .  9  The  numerical  values 
published  here  were  computed  using  a  program  developed  at  UTSI. 
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In  this  appendix  a  ray  trace  is  developed  to  find  the  value  of  f 
for  small  angles  of  incidence  (0^$  20  )  and  spherical  particles.  Consider 
the  two  rays  in  figure  7.  The  assumption  is  made  that  Snell’s  law  of  re¬ 
fraction  can  be  written  as 


SIN0  =  mSIN0  (a 

1  K 


6i  “  n6R 
0±  <20° 


(b 


(56 


Elementary  geometrical  arguments  show  that  the  angles  of  interest  can  be 
written  as 

3'  =  2(0.-  0R1)  (57 

Y’  =  20r1-  0±  (58 

Using  the  approximation  for  Snell fs  law  (equation  56b)  equation  57  and 
58  may  be  written  as 
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(60 


Y’ 


=  9±  M 
i  l  m 


Inspection  of  figure  7  shows  that  p  and  q  can  be  written  as 


p  =  rsiny' 


(61 


q  =  p/tanfi’ 


(62 


q  ■  rsiny' /tan$* 


(63 


For  small  angles 


siny'  -  y1 


(64 


tanfi'  - 


(65 


cosy '  -  1 


(66 


Using  equations  64  and  65  in  equation  63  yields 

_  r  [2  -  ml 

q  "  2  i^Tj 


(67 


Equation  66  allows  f  to  be  written  as 


f  -  r  +  q 


(68 
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Substituting  equation  67  into  equation  68  produces  the  desires  relationship: 


f 


(69 
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Figure  1.  Probe  Volume  Co-ordinate  System  Geometry  for 
Interferometric  Particle  Size  Analysis. 
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Figure  2.  Co-ordinate  System  Geometry  for  Visibility 
Estimate  for  Refracting  Spheres. 
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Figure  3.  114  micrometre  glass  sphere  in  a  fringe  pattern  and 

the  resulting  scattered  field  pattern. 
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Focal  Point  Separation 
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Figure  5.  Example  Signal  Phase  and  Visibility  for  a  Large  Refracting  Sphere. 
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Figure  6*  Ray  Trace  Co-ordinates  for  the  Estimation  of  f 
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A  MODEL  FOR  ESTIMATING  PARTICLE 
NUMBER  DENSITY  USING  DIFFERENTIAL 
DOPPLER  LASER  VELOCIMETERS 

W.  M.  Fanner 

ABSTRACT 

A  model  is  developed  which  may  be  used  to  predict  particle  number 
density  using  data  from  currently  available  differential  Doppler  laser 
velocimeter  optical  systems  and  burst  counter  signal  processors.  The 
results  may  also  be  used  to  predict  relative  system  performance  as  the 
LDV  system  parameters  used  in  the  model  change. 

INTRODUCTION 

The  estimation  of  particulate  number  density  in  a  moving  fluid  is 

a  particularly  difficult  experimental  problem  subject  to  numerous  sources 

of  error.  For  example,  mechanical  sampling  systems  such  as  cascade 

impactors  must  sample  isokinetically  and  at  a  fixed  flow  rate.  These 

kinds  of  measurements  yield  number  density  estimates  from  computations 

involving  weight  measurements  for  each  stage  of  the  impactor  corresponding 

12  .  . 

to  "equivalent"  aerodynamic  particle  diameters.  ’  The  resulting  estimate 
is  thus  not  a  direct  numeric  count  but  rather  an  "equivalent"  number 
density.  Because  mechanical  samplers  often  introduce  unacceptable  per¬ 
turbations  in  the  flows  being  measured,  optical  and  other  more  exotic 
techniques  have  been  developed,  each  having  its  own  peculiarities  and 
limitations.  These  methods  have  ranged  from  multiple  wavelength 
transmissometer  schemes  to  sophisticated  multiple  scatter  angle  techniques 
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involving  highly  sophisticated  inversion  schemes  and  considerable 
computational  effort.  Single  particle  optical  counters  which  determine 
a  particle  size  from  a  scattered  light  measurement  and  maintain  an 
accurate  particle  count  rate  for  a  fixed  sample  flow  rate  are  among  the 
most  accurate  devices  for  estimating  numeric  number  densities.  They 
are,  however,  cumbersome  to  use  except  in  the  most  benign  environments 
because  of  the  need  to  draw  a  fluid  sample  through  their  optically 
sensitive  region.  Furthermore,  when  number  densities  become  large  , 
it  is  necessary  to  dilute  the  concentration  in  order  to  reduce  the 
probability  that  more  than  one  particle  will  exist  in  the  sample  volume 
during  a  measurement  interval.  In  hostile  environments,  such  as 
encountered  in  combustion  or  chemically  active  flows  or  in  flows  with 
particles  that  may  be  affected  by  the  sampling  system,  it  is  not  practical 
to  use  these  kinds  of  optical  counters.  Differential  Doppler  laser 
velocimeter  systems  are  appropriate  for  use  in  these  kinds  of  particle 
climates.  The  sampling  volumes  of  these  devices  may  be  modeled  as 
periodic  structures  resulting  from  sets  of  localized  interference  fringes. 
The  number  of  cycles  and  signal  shape  in  a  signal  generated  by  a  particle 
scattering  light  from  this  fringe  structure  may  be  related  both  to  the 
size  and  the  trajectory  of  the  particle.  As  a  result,  signal  processing 
electronics  has  been  developed  which  can  test  signals  in  terms  which  can 
at  least  be  heuristically  related  to  trajectory  through  the  sample  volume 
or  the  presence  of  two  or  more  particles  generating  the  signal.  Thus, 
these  optical  systems  do  not  rely  on  plumbing  or  dilution  to  restrict 
particle  flow  through  the  sampling  volume.  Rather,  it  is  common  practice 
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to  focus  these  devices  into  particle  laden  flows  and  to  let  the  particle 

flux  arrive  as  the  flow  velocity,  concentration,  and  spatial  distribution 

**v,  * 

dictate.  The  signal  processing  electronics  and  signal  test,  logic  then 
sample  the  particle  flux  in  a  manner  analogous  to  the  functions  imposed 
by  the  sampling  and  dilution  plumbing  of  single  particle  optical  counters. 
The  purpose  of  this  paper  is  to  present  a  model  that  shows  how  data 
obtained  with  a  differential  Doppler  laser  velocimeter  (LV)  may  be  used 
to  estimate  particle  number  density  in  otherwise  inaccessible  flow  fields. 
The  method  utilizes  data  which  is  readily  available  from  most  existing  LV 
systems  and  should  therefore  be  straightforward  to  implement  on  a  wide 
variety  of  experiments  where  these  devices  are  in  use  and  number  density 
data  are  of  interest.  In  addition  to  providing  a  number  density  estimate 
from  LV  data,  the  method  may  also  be  used  to  address  a  problem  often 
encountered  in  the  application  of  the  LV.  The  problem  is  to  estimate  how 
often  data  can  be  acquired  for  a  specified  number  density  given  the 
constraints  of  the  optical  and  signal  processing  system. 

Because  the  approach  taken  in  this  paper  is  somewhat  different  than 
normally  encountered  in  dealing  with  number  density  limitations  on  LV 
system  performance,  a  review  is  first  given  of  the  results  of  a  typical 
LV  system  measurement  of  a  particle  laden  flow.  Next,  statistics  normally 
associated  with  this  process  are  examined  in  the  various  limits  which 
might  be  encountered  experimentally.  Factors  which  affect  the  applicability 
of  these  statistics  are  examined.  A  model  is  then  proposed  to  circumvent 
the  difficulty  imposed  by  these  factors.  Finally,  experimental  results 
are  presented  which  suggest  that  the  proposed  model  can  be  effectively 
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used  to  estimate  particle  number  density. 


TYPICAL  LV  EXPERIMENT 

A  typical  LV  experiment  might  involve  an  optical  system  of  the 
type  shown  in  figure  1.  Two  coherent  beams  are  made  to  cross  and  focus 
at  some  position  in  a  particle  laden  flow.  The  scattered  light  collector 
will  be  assumed  to  be  sufficiently  large  that  individual  particles  can 
be  detected  as  they  traverse  the  beam  cross  region.  This  region  is 
commonly  called  the  probe  volume.  For  TEMq  laser  beams  and  angles 
between  the  beams  of  about  10°  or  less,  this  volume  will  appear  to  be 
an  ellipsoid  of  revolution.  In  most  experiments,  the  flow  velocity  is 
not  constant  but  varies  randomly  about  some  mean  value.  Normally,  the 
particles  in  the  flow  may  be  distributed  in  particle  size,  shape,  and 
composition.  Depending  on  the  drag  coefficients  which  may  be  ascribed 
to  the  ensemble  of  particles  measured  and  the  fluid  velocity  parameters, 
the  particles  may  follow  the  fluid  velocity  or  "slip”  relative  to  it. 
Normally,  it  is  assumed  that  the  particles  follow  the  fluid  velocity 
well  within  the  uncertainties  acceptable  for  the  experiment.  In  a 
typical  experiment,  the  photodetector  is  set  to  a  fixed  gain  level  which 
is  commensurate  with  amplifier  linearity  used  in  the  signal  processor 
electronics.  A  signal  threshold  level  is  set  for  a  value  that  minimizes 
or  eliminates  the  time  spent  by  the  signal  processor  attempting  to 
measure  sources  of  non  flow  related  noise.  If  it  is  assumed  that  the' 
threshold  is  set  such  that  the  processor  does  not  spend  time  attempting 
to  measure  noise,  then  for  a  sample  time  At,  with  valid  measurements 
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(valid  in  the  sense  that  all  signal  processor  requirements  are  met) 
the  following  factors  account  for  the  value  of  At. 

1.  The  amount  of  time  required  for  particles  to  traverse  the 
sample  volume  -At^ 

2.  The  amount  of  time  during  which  no  particles  are  present  in 
the  sample  volume  “At^ 

3.  The  amount  of  time  due  to  2  or  more  particles  are  simul¬ 
taneously  in  the  sample  volume  -At^ 

4.  The  amount  of  time  due  to  single  particles  producing  signals 
with  too  few  cycles  to  be  measureable  (normally  associated 
with  very  small  particles  or  particles  passing  near  the  edge 
of  the  sample  volume)  -At^ 

5.  The  amount  of  reset  time  required  by  the  signal  processor  and 
data  acquisition  system  after  each  measurement  attempt  -AtR 

Thus  At  may  be  expressed  as: 

At  =  At^  +  At2  +  At j  +  AtR  (1) 

At^  may  be  expressed  as: 

n 

At-  -  £  Nt  ,T .  (a 

1  i=l  Li  ± 

or  (2) 

Atj  =  AtP(l)  (b 


where  N  is  the  number  of  cycles  generated  in  the  ith  signal  and  T 

til 

is  the  signal  time  period.  P(l)  is  the  probability  of  1  particle  being 
present  in  the  sample  volume.  may  be  expressed  as: 


AtR  "  (NE"1)TR 


(3) 
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where  T  is  the  signal  processor  reset  time  and  N_  is  the  total  number 
K  E 

of  attempted  measurements  made  by  the  signal  processor.  No  similar 
simple  relationships  exist  for  At ^ >  At^,  and  At^.  These  times  can  be 
estimated  using  statistical  estimators.  For  example,  At2  may  be  written 
as: 

At2  -  AtP(O)  (4) 

where  P(0)  is  the  probability  that  zero  particles  are  present  in  the 
sample  volume.  Similarly,  At^  might  be  written  as: 

At3  *  At  P(n  >2)  (5) 

where  P(n  >2)  is  the  probability  that  two  or  more  particles  are  present 
in  the  sample  volume.  At^  might  be  computed  using  P(l)  and  a  weighting 
factor  to  account  for  the  fact  the  sample  volumes  over  which  these 
measurements  can  be  made  is  probably  much  smaller  than  that  corresponding 
to  Atj  or  Aty 

Experimentally,  it  is  not  unusual  to  find  N  »  N  and  T  «N  ,T,  . 

E  R  Li  i 

It  is  thus  possible  that  At^  can  make  a  significant  contribution  to  At. 

It  is  often  observed  that  in  flows  with  moderate  number  densities  of 
particles  that  At^/At  <<  1,  and  that  N^/N^  (N^  being  the  number  of  signals 
generating  At^)  generally  varies  between  0.1  -  0.8  and  seems  to  depend  on 
number  density.  This  suggests  that  a  statistical  model  might  relate 
particle  number  density  to  N^/N^.  The  next  section  examines  a  statistical 
model  commonly  used  to  analyze  LV  optical  systems. 


106 


COMMON  STATISTICAL  ANALYSIS  OF  FINDING  n  PARTICLES  IN  THE  LV  SAMPLE  VOLUME 

It  is  generally  assumed  that  the  probability  of  finding  n  particles 

in  the  LV  sample  volume  is  fully  equivalent  to  that  posed  by  Chandrasekhar 

for  N  particles  randomly  distributed  in  a  volume  V  which  is  sampled  by  a 
4  5 

small  volume  v.  *  In  this  case,  Chandrasekhar,  starting  from  the 
Bernoulli  distribution,  derives  the  Poisson  distribution  for  simultaneously 
obtaining  n  particles  in  v  as 


P(n) 


(Np)n  -Np 
n! 


(6) 


where  p  is  the  probability  of  finding  a  particular  particle  in  v.  If  the 
particles  are  randomly  distributed  throughout  V  then  by  definition  of  a 
randomly  variable ^ 

p  =  v/V  (7) 


The  value  Np  is  usually  written  as  p^v  where  is  the  number  of  particles 
per  unit  volume.  P(n)  may  then  be  written  in  the  form  normally  seen  for 
equation  6 


P(n) 


(pNv>  pnv 

-  e 

n! 


(8) 


For  later  reference  it  is  to  be  noted  the  probability  of  obtaining  2  or 
more  particles  in  the  sample  volume  may  be  written  as 


P<2)  =  1 


~PNV  -pNV 

e  -  PNve 


(9) 
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Note  that 


lim 


P 


oo 

N 


lim 

V° 


P(2)  =  1 


P(2)  =  0 


(10) 


(ID 


In  the  case  where  p^v  >>  1,  the  Poisson  distribution  approaches  a  Gaussian 
distribution.  Chandrasekhar  shows  that  P(n)  can  be  expressed  as 


P  (n) 


x  ~(n  -  PNv)2/2pNv 

-  e 

V  2ttpnv 


(12) 


These  equations  are  relatively  straight forward  and  present  no 
interpretational  difficulties  until  attempts  are  made  to  apply  them  to 
experimental  LV  data.  The  fundamental  experimental  parameter  in  equations 
8  and  12  is 


PNv  =  N  v/V  (i3) 

v  and  V  are  well  defined  quantities  in  the  mathematics.  Experimentally, 
this  may  not  be  so.  Consider  the  fact  that  v  may  be  a  time  dependent 
variable.  For  example,  v  for  a  single  particle  can  be  defined  as 

v  =  v  a  t  (14) 

pv 

where  v  is  the  particle  velocity  component  normal  to  the 
sample  volume  cross  sectional  area  a  and  t  is  the  particle 
transit  time  across  the  sample  volume.  If  account  is 
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taken  of  the  fact  that  the  particles  sampled  may  have  a  speed  distri¬ 
bution  as  well  as  a  size  distribution,  then  the  sample  volume  for  a 
particle  having  the  ith  velocity  and  jth  particle  size  is 


v 


ij 


V.a  4 

~i  pvj 


t .  . 
iJ 


(15) 


a  and  t  have  j  subscripts  because  it  is  presently  understood  that  an 
LV  sample  volume  depends  on  the  particle  size  (or  more  specifically  the 
scattering  cross-section)  producing  the  signal. ^  Hence,  it  is  to 
be  expected  that  v  will  vary  from  one  particle  to  the  next.  A  second 
difficulty  arises  in  attempting  to  define  the  total  volume  of  aerosol 
sampled,  V.  For  example,  if  V  were  defined  in  a  fashion  similar  to 
equation  14,  then 


-ru 


a  v 


( ^  At,  4-  a  At9  +  a  AtQ  +  a  At,  )f  (v)f  (a)dadv 
pv1  1  p  v  2  2  pv3  3  pv^  4' 


where  f (v)  and  f (a)  are  the  probability  density  functions  for  the  velocity 
and  extinction  cross-sections  for  the  measured  ensemble.  The  resolution 
of  equation  16  depends  not  only  on  assuming  functions  for  f  (v)  and  f  (cr) 
but  also  in  defining  a  when  no  particles  are  present  or  when  2  or 
more  occupy  the  sample  volume. 

Equation  15  shows  that  the  sampling  volume  varies  particle  to 
particle  as  a  function  of  both  particle  size  and  speed.  Except  for 
artificially  contrived  cases,  it  is  not  constant  in  value.  Equation  16 
illustrates  the  problem  associated  with  attempting  to  estimate  V  from 
experimental  data  and  models  for  the  velocity  and  size  distributions. 
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Erdmann  and  Gellert  have  pointed  out  that  when  the  velocity 
is  distributed  the  Poisson  distribution  describing  the  arrival  rate 
of  particles  crossing  the  sample  volume  must  be  replaced  by  the  compound 
Poisson  distribution.5  In  their  case,  the  compound  Poisson  distribution 
for  some  distribution  of  velocity  f (v)  is  given  by 

Q(n,XT)  -  f  {Ut) "exp (-XT)  .  j  f  (j)dj  (1?) 

v  n ! 

where  X  is  the  mean  particle  arrival  rate  and  T  is  the  sample  time  of 
interest.  Feller  indicates  how  compound  Poisson  distributions  lead  to 

.  g 

solutions  for  Brownian  motion  and  "randomized  random11  walk  problems.  For 
the  application  under  consideration  here,  the  compound  Poisson  process  yields 

f  (p  v)  n  -PMV 

Q(n,PNv)  =  \{  -  e  }  <f>(v)dv  (18) 

*  n! 

v 

where  v  is  the  sample  volume  and  $(v)  is  the  sample  volume  distribution 
which  may  depend  on  both  the  particle  scattering  cross-section  and  velocity 
distribution.  Hence,  Q(n,PNv)  will  be  a  strong  function  of  experimental 
conditions  and  instrumentation  response.  To  formally  solve  equation  18,  it 
is  necessary  to  assume  a  functional  dependence  for  <j>(v).  In  this  work, 
rather  than  directly  postulating  a  4>(v)  function,  a  somewhat  different 
approach  will  be  taken  which  has  been  found  to  yield  reasonably  good 
agreement  with  experiment  and  provides  some  physical  insight  into  the 
measurement  process.  In  this  approach,  it  is  first  assumed  that  at  some 
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instant  in  time  the  particle  ensemble  to  be  measured  is  randomly 
distributed  in  the  total  measured  volume.  As  the  sample  volume  is 
moved  through  the  ensemble,  it  will  be  assumed  to  vary  (because  of 
velocity  and  particle  size)  in  a  fashion  that  can  be  described  by  a 
random  walk  process.  Such  a  sample  volume  will  be  called  a  "Random 
Walk  Sample  Space."  In  the  following  section,  the  probability  of  two 
or  more  particles  existing  within  some  distance  of  each  other  for  each 
of  three  orthogonal  coordinates  in  the  random  walk  sample  space  is 
computed.  It  will  be  assumed  that  the  logic  circuits  of  the  LV  signal 
processor  are  sufficiently  accurate  to  reject  all  signals  resulting  from 
two  or  more  particles.  Obviously,  the  accuracy  with  which  this  can 
be  accomplished  is  a  major  limitation  to  the  accuracy  of  this  model. 

(Our  experience  has  been  that  such  circuits  are  fooled  by  broadband 
Gaussian  noise  in  less  than  1%  of  the  measurement  attempts  which  might 
correspond  to  a  multiple  particle  signal).  The  resulting  probability 
distribution  for  particle  separation  in  three  dimensions  and  assuming 
accurate  multiple  particle  detection  lead  to  the  definition  of  a  signal 
acceptance  ratio,  A.  "A"  specifies  the  ratio  of  number  of  particles 
measured  to  those  detected.  "A"  is  then  taken  as  the  sampling  efficiency 
of  the  LV  and  is  a  quantitative  measure  of  how  well  the  LV  performs  as 
a  function  of  number  density.  Conversely,  if  the  LV  optical  system 
parameters  and  signal  processor  characteristics  are  carefully  specified, 
then  a  measurement  of  "A"  can  be  used  to  estimate  particle  number  density. 


Ill 


PROBABILITY  OF  TWO  PARTICLES  SIMULTANEOUSLY  IN  THE  RANDOM  WALK 
SAMPLE  SPACE 

The  development  of  the  probability  density  distribution  for 

two  particles  simultaneously  in  a  random  walk  sample  space  is  based 

on  the  well  known  random  walk  analysis  commonly  found  in  text  books 

6  8 

on  statistics.  *  Because  the  functional  dependences  of  the  resulting 
distribution  are  not  immediately  evident,  the  steps  leading  to  the 
result  are  briefly  sketched. 

Because  the  particles  are  assumed  to  be  randomly  distributed  in 
space  and  a  random  walk  is  performed  to  reach  one,  the  probability 
P(x,y,z)  of  a  particle  existing  at  some  x,y,z  position  in  a  rectangular 
coordinate  system  may  be  written  as 

P (x,y , z)  *  P (x)  P (y)  P(z)  (19) 

where  P(x)  etc.,  is  the  probability  that  the  particle  position  is  at 
the  coordinate  x.  Consider  figure  2,  a  dimension  in  the  x,  y,  z  coordinate 
system,  is  divided  into  a  set  of  cells  of  dimension  G.  The  occupied 
cell  is  approached  by  advancing  or  retreating  in  unit  cells  with  equal 
probability,  i.e.  the  probability  of  advancing  a  unit  cell  is  1/2. 

Let  some  arbitrary  cell  be  chosen  as  the  origin  and  count  cells  to  the 
left  until  a  cell  is  found  which  is  occupied.  Let  this  number  of  cells 
be  p.  Next,  let  q  be  the  number  of  cells  counted  from  the  origin  before 
the  next  successive  occupied  cell  is  found.  Define  the  total  number  of 
cells  involved  as  given  by  n  while  the  number  of  cells  between  the 


particles  is  m,  then 


P  +  q  =  n  (18) 

p  -  q  =  m  (19) 


The  probability  that  a  sequence  n  leads  to  a  value  of  m, 
by  the  Bernoulli  equation 


P  is  given 
n  j  in 


p 

n,m 


(20) 


Applying  Stirling's  formula  for  n!  and  manipulating  variables,  there 
results  the  Laplace  de  Moii-ve  limit^ 


P 

n,m 


exp (-m  /2n) 


|  m  |  «  n 


(21) 


It  is  convenient  to  make  the  following  identities 


x  =  me 


(22) 


X  = 


n£ 

2 


(23) 


where  x  is  the  spatial  displacement  between  particles  and  X  is  the 
arithmetic  average  of  the  distance  between  the  particles.  Substituting 
Equations  22  and  23  into  Equation  21  results  in 
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Let  P  (x,X)dx  be  the  probability  of  a  distance  x  to  x  +  dx  existing 


between  particles  when  the  average  distance  is  X, 
is  approximately  constant.  Then  P'(x,X)  dx  is  P^  m 
number  of  values  of  m  occurring  in  dx  for  fixed  n, 


For  small  dx,  P 

n,m 

multiplied  by  the 


P(x,X)dx  =  P 

n,m 


dx 

£ 


(25) 


A  factor  of  ^  enters  since  for  fixed  values  of  n,  values  of  m  are  either 
all  odd  or  all  even.  Equation  24  can  now  be  expressed  using  Equation  25 
as 


P(x,X)dx 


1 

(TrAeX)*4 


dx 


(26) 


It  is  assumed  that  the  arithmetic  average  distance  between  particles  along 
a  specified  direction  for  a  given  number  density  is  given  by 


X  = 


(27) 


The  cell  length  e  is  chosen  to  be  the  fringe  period  6  of  the  LV  system 
since  this  is  the  basic  resolution  element  through  which  the  logic 
circuits  in  the  LV  processor  must  decide  if  two  or  more  particles 
occupy  the  sample  volume.  With  these  identifications 


.  4tt5  . 


exp 


x 


P(x,X)dx 


dx 


(28) 


It  is  convenient  to  define  a  as 


1/3 

a  =  PN  M5 


(29) 


With  these  results  and  definitions,  the  probability  of  finding 
2  or  more  particles  in  some  increment  of  space  AxAyAz  when  the  particles 
are  separated  by  some  average  distance  X  is  given  by 


P (AxAyAz) 


Ax  Ay  Az 


P (x,X)  P(y,X)  P (&,X)  dxdydz  (30) 


o  o  o 


Because  the  integrals  are  separable  and  identical  in  form,  differing 
only  in  the  limits  of  interest,  only  one  is  explicitly  evaluated.  Using 
Equations  28  and  29  for  this  evaluation  gives 


(31) 


which  readily  yields 


P(x,X)  dx  =  erf  (-yT  Ax) 


(32) 


An  accurate  approximation  for  a  closed  form  solution  to  Equation  32  may 


lib 


Q 

be  obtained  by  using  an  approximation  due  to  Menzel 


erf(b)  =  -  exp  (-4b2 /it)  j 


(33) 


Therefore,  the  solution  to  the  integral  is 


Ax 


o 


(34) 


In  the  limit  as  Ax  -*•  the  integral  has  a  value  of  Since  in  the 
limit  it  is  required  that  the  particles  be  somewhere  in  the  box  AxAyAz, 
the  probability  of  finding  the  particles  separated  by  some  distance  Ax  is 
normalized  to  the  value  for  infinite  separation. 


p(Ax) 


(x,X) 


o 


dx 


P (x,X)  dx 


p(Ax) 


Ax 


P(x,X) 


o 


dx 


(a 


(35) 


(b 


p (AxAyAz) 


P (AxAyAz) 

P  (°°  ,co,0°) 


(c 


llo 


Using  Equations  34  and  35  in  Equation  30  and  manipulating  variables, 
a  functional  relationship  for  p (AxAyAz)  is  obtained 


p (AxAyAz) 


(36) 


p (AxAyAz)  can  be  simplified  to  a  somewhat  more  convenient  form 


p (AxAyAz) 


=  8^  exp 


X/3 

— —  (  Ax2  4*  Ay2  4*  Az2  )1  x 

4tt6  J 


sinh 


(37) 


p (AxAyAz)  represents  the  probability  that  2  or  more  particles  will  be 
found  in  a  volume  AxAyAz*  By  using  values  of  Ax,  Ay  and  Az  appropriate 
to  the  LV  probe  volume  or  sample  space,  an  estimate  can  be  made  of  how 
often  two  or  more  particles  may  be  expected  in  the  probe  volume.  Note 
that  Equation  36  or  37  shows  that 


lim  p (AxAyAz)  =  1 


(a 


(38) 


lim  p (AxAyAz)  =  0 

V  0 


(b 
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which  is  a  result  predicted  by  Poisson  statistics  and  should  also 
be  reasonably  predicted  for  these  statistics.  It  should  also  be 
noted  that  while  Ax,  Ay,  Az  appear  separately  in  equation  37  and  not 
directly  as  a  volume,  equation  37  does  not  predict  that  p (AxAyAz) 
is  shape  dependent,  p (AxAyAz)  is  perfectly  symmetric  in  Ax,  Ay, 

Az.  As  a  result,  there  is  no  preferred  shape  dependence  for  any 
one  dimension  or  equivalently  for  a  particular  volume  shape. 


THE  ACCEPTANCE  RATIO 

The  results  from  the  previous  section  can  now  be  used  to  define 
an  experimental  parameter  called  the  acceptance  ratio*  A.  A  is 
defined  as 

A  =  dumber  of  Signals  Accepted  for  Measurement  .  . 

Number  of  Signals  Detected  ^  ' 

where  "number  of  signals  detected"  Is  understood  to  mean  those  which 
generate  sufficient  signal  cycles  to  be  measured  -  both  those  resulting 
from  single  particles  and  those  resulting  from  2  or  more  particles. 

If  it  is  assumed  that  the  logic  circuits  in  a  burst  type  LV  signal 
processor  are  totally  effective  in  rejecting  all  multiple  particle 
signals,  then  if  N  is  the  total  number  of  signals  detected,  the  number 
of  signals  rejected  is  Np (AxAyAz)  where  Ax,  Ay,  Az  are  appropriate 
probe  volume  or  sample  space  dimensions.  With  these  assumptions  and 
definitions.  Equation  can  be  used  to  express  A  as 


A  =  1  -  p(AxAyAz) 

(a 

A  =  P(0)  +  P(l) 

(b 

(40) 


where  P(0)  and  P(l)  are  the  probabilities  of  finding  no  particles  and 
a  single  particle  in  p(AxAyAz).  In  order  to  explicitly  express  A  in 
terms  of  LV  system  parameters,  it  is  convenient  to  make  the  following 
identifications.  Assume  an  LV  system  similar  to  that  in  Fig.  1  capable 
of  measuring  one  velocity  component,  then  let  Ay  be  the  dimension 
parallel  to  the  velocity  component  measured.  Express  Ay  as 


Ay  =  N  6  (41 ) 

Li 

where  N  is  taken  to  be  the  average  number  of  observable  cycles  generated 
Li 

in  the  LV  signal.  Let  Ax  be  the  dimension  perpendicular  to  Ay  and  the 
bisector  between  the  beams  (i.e.  in  the  fringe  model  of  the  LV  it  is 
parallel  to  the  fringe  planes).  Write  Ax  as 

Ax  =  kN  6  (42 ) 

Li 

where  k  is  a  constant  which  specifies  the  relative  size  of  Ax  and  Ay. 

In  practice  Ax  _<  Ay.  This  may  result,  for  example,  from  logic  constraints 
in  the  signal  processor  or  slit  apertures  placed  in  the  receiver  optics 
to  limit  the  size  of  the  probe  volume.  Az  is  assumed  parallel  to  the 
bisector  between  the  beams  and  may  be  expressed  in  terms  of  Ax  and  Ay, 

For  most  LV  systems  currently  in  use,  Az  »  Ay.  For  example,  if  Az 
were  made  to  correspond  to  the  length  of  the  transmitted  probe  volume 


lly 


i 


-2 

corresponding  to  the  average  e  intensity  contour,  then 


Az 


*y 

a 


(43) 


where  a  is  the  angle  between  the  beams.  A  typical  value  of  a  is  less 

1  / 

than  0.1  making  Az  =  40Ay.  When  3Az2/2tt6  is  computed  for  the  range 
of  values  commonly  encountered  in  most  LV  systems  and  applications , (i.e. 
Az>1006,  >  5cc  \  5  >  6  pm)  it  is  found  that 


slnh 


V* 

2k& 


exp 


Va 

PB 


1  i 


(AA) 


and 


■/. 

2tt6 


>  2 


(A5) 


It  follows  for  many  applications  that  A  will  be  independent  of  Az.  With 
these  definitions  and  approximations,  A  can  be  explicitly  expressed  as 


1  - 


2  exp(  (1  +  k2))  £sinh  (y)  sinh  (k2y)j 

Vs 


D  N2  6 
L 

2tt 


(A6) 


(A7) 
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Equations  46  and  47  show  that  if  N  ,  6,  and  k  (all  experimental  constants 

Li 

for  a  given  LV  system)  are  known,  then  a  determination  of  A  can  be  used 

to  compute  p  .  Conversely,  specification  of  p  with  a  required  k,  N  , 

N  N  Li 

and  5  provides  an  estimation  of  the  acceptance  ratio.  If  A  is  taken  as 
a  measure  of  system  performance,  then  examination  of  the  response  of  A 
for  variation  in  system  parameters  will  show  how  to  optimize  the  system 
for  expected  operational  conditions. 

LV  SYSTEM  PERFORMANCE 

Figure  3  plots  A  as  a  function  of  y  f°r  selected  values  of  k 

between  1.0  and  0.2.  As  the  figure  shows,  for  fixed  and  6, 

significant  increases  in  acceptable  p^  for  given  values  of  A  can  be 

achieved  by  reduction  of  the  probe  volume  cross-section  through  apertures 

(a  fact  well  known  to  anyone  who  has  used  an  LV) .  Figure  4  graphically 

illustrates  this  dependence  by  plotting  y  as  a  function  of  k  for  selected 

values  of  A.  The  figure  shows,  for  example,  that  if  A  =  0.1,  then 

reducing  k  from  1  to  0.2  when  N  and  6  are  held  constant  results  in  an 

L 

effective  increase  of  acceptable  p^  by  a  factor  of  nearly  8,000. 

In  applying  an  LV  to  a  turbulent  flow  measurement,  it  may  be 
necessary  to  add  particles  to  the  flow  to  satisfy  sampling  rate  require¬ 
ments.  However,  Fig.  3  shows  that  as  the  particle  number  density  (i.e.  y) 
increases,  the  acceptance  ratio  decreases.  Eventually,  the  number 
density  will  reach  a  level  where  the  data  rate  actually  decreases  with 
increasing  number  density  because  the  acceptance  ratio  has  become  so 
small.  As  the  number  density  increases,  a  point  will  be  reached  where 
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numerous  particles  always  exist  in  the  sample  space  and  either  it 
becomes  expedient  to  switch  to  a  local  oscillator  optical  system 
or  to  use  signal  processors  which  optimally  function  with  continuous 
type  signals.  The  above  analysis  can  be  used  to  estimate  the  maximum 
number  density  which  can  continue  to  increase  the  data  rate  when  all 
other  system  parameters  are  fixed.  Note  that  the  data  rate  T\  can  be 
written  as 


n  -  NA 


(48) 


where  N  is  the  rate  at  which  particles  pass  through  the  sample  volume. 

Assume  that  the  time  scales  involved  are  such  that  A  can  be  assumed 

constant  and  that  the  flow  velocity,  V,  is  constrained  to  the  dimension 

normal  to  the  cross-sectional  area  a  of  the  sample  space,  then  N 

pv 

can  be  written  as 


Px,  V  a 
~  pv 


(49) 


combining  Equations  48  and  49,  taking  the  derivative  with  respect  to 
P^,  and  making  use  of  Equation  47,  there  results 


(50) 

/ 


X]  is  an  extreme  when  the  bracketed  term  in  Equation  50  is  zero.  Using 

Equation  48  and  performing  the  indicated  algebra  in  Equation  50  ,  it  is 

3  n 

found  that  « -  is  an  extreme  when 

3pN 


IO 
^4, 


Y  ,  , 

_  6  (1  -  coth  (Y))  +  k2  (1  -  coth  QOr)) 

A  X  (1  -  coth  (Y))  +  k2  (1  -  coth  (kzy) )  -  1 
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Figure  5  plots  A  as  a  function  of  Y  for  various  values  of  k  for  optimum 
data  rate.  The  figure  shows  the  interesting  result  that  after  k  reaches 
about  0.5,  A  becomes  nearly  constant.  This  means  that  the  acceptance 
ratio  or  equivalently  the  LV  sampling  efficiency  cannot  be  improved  by , 
for  example,  further  reduction  of  the  sample  space  cross-section  via 
apertures.  However,  further  reduction  of  k  will  provide  an  increase  in 
the  data  rate  for  increasing  number  density. 

For  the  case  where  A  is  chosen  to  yield  a  maximum  value  for  f| 
as  a  function  of  p^,  Y  is  assumed  constant.  Solving  Equation  47  for 
with  Y  constant  yields 

.  I27iY)l  (52) 

N  Nt6  6> 

Equation  5  2  shows  this  extreme  value  of  p^  is  much  more  drastically 
affected  by  the  number  of  signal  cycles  generated  by  a  scattering  particle 
than  by  the  absolute  magnitude  of  the  LV  fringe  period.  Hence,  in  seeking 
to  optimize  data  rates  for  a  particular  LV  geometry ,  it  is  most  expedient 
to  reduce  the  number  of  cycles  in  the  signal  rather  than  to  change  the 
fringe  period. 
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NUMBER  DENSITY  DETERMINATION 


In  this  section,  the  results  of  experimental  applications  of 
the  expression  derived  for  A  in  the  previous  section  are  described. 

It  is  very  difficult  to  produce  a  direct  quantitative  comparison  of 
experiment  with  theory  since  there  are  no  generally  accepted  calibration 
standards  for  number  density.  In  this  section,  the  results  of  three 
experimental  comparisons  are  described.  The  results,  while  inferential, 
suggest  that  the  previously  computed  expression  for  the  acceptance  ratio 
does  yield  reasonably  good  results. 

The  first  set  of  results  reported  here  were  obtained  with  a 
particle  sizing  interferometer  (PSI)  in  a  laboratory  chamber  test  and 
in  an  outdoor  field  test.  Details  of  these  tests  may  be  found  else¬ 
where.^  Only  a  brief  sketch  of  the  experiments  is  presented  here. 

The  optical  geometry  for  the  PSI  is  shown  in  figure  6.  The  system 
was  arranged  to  measure  a  particle  size  range  of  0.3  -  6.0  micrometres. 
The  smallest  particle  size  detectable  with  the  system  was  estimated  to  be 
approximately  0.2  micrometres  for  water  droplets.  This  PSI  used  a  5 
milliwatt  HeNe  laser  and  F/2.2  receiving  lens.  Data  obtained  with  the 
PSI  signal  processor  and  data  acquisition  system  in  these  experiments 
were  recorded  on  magnetic  disks  for  hard  copy  presentation  and  data 
storage.  The  PSI  measured  particle  velocity,  signal  magnitude  (which  is 
also  a  measure  of  particle  size),  the  PSI  size  parameter-visibility,  and 
the  acceptance  ratio. 

The  laboratory  measurements  were  conducted  in  a  sealed  box  with 
an  atmosphere  initially  at  ambient  temperature  and  humidity  conditions. 


The  PSI  system  was  placed  directly  inside  the  chamber.  Filter  samples 
and  samples  obtained  with  an  Andersen  Cascade  Impactor  were  obtained 
via  samples  drawn  from  wall  ports. 

Phosphorus  and  HC  smokes  were  generated  by  igniting  small  samples 
at  one  end  of  the  chamber.  Fog  oil,  while  not  hygroscopic,  was  easy  to 
generate  using  a  small  commercially  available  fog  oil  generator  and  it 
also  was  measured  using  both  PSA  systems.  Hass  concentration  as  a 
function  of  time  was  obtained  using  filter  samples  which  trapped  a 
portion  of  the  smoke  withdrawn  from  the  chamber  at  a  fixed  flowrate. 

Figure  7  shows  a  comparison  of  concentration  estimates  obtained 
mechanically  by  weight  and  that  obtained  using  PSI  acceptance  ratio 
data.  The  concentration,  C,  was  computed  from 


C  - 


7  pVn 

o 


(53) 


where  p  is  the  material  density,  and  is  the  third  moment  of  the  size 
distribution.  A  major  uncertainty  in  the  comparison  is  the  material 
density  of  the  smokes  measured.  Orthophosphoric  acid  (density  =  1.28  gm/ 
at  82%  RH)  is  assumed  for  red  phosphorus  and  for  fog  oil  a  density  of 
0.9  gm/cc  is  assumed. 

As  Figure  7  shows,  agreement  between  the  concentration  estimate 
using  the  acceptance  ratio  data  and  that  determined  by  weighing  is 
surprisingly  good. 

It  is  of  interest  to  compare  the  number  density  determined  by  the 
acceptance  ratio  using  equation  46,  i.e.  for  a  random  walk  sample  space 


for  the  data  in  figure  7  and  that  obtained  using  commonly  used  Poisson 
or  Gaussian  statistics.  Poisson  statistics  yield  an  acceptance  ratio 
given  by 

~PNV 

A(P)  =  (1  +  PNv)  e  (54) 

Gaussian  statistics  (p^v  »  1)  yield 

exp (- (2-p  v) 2/2p  v) ) 

A (G)  =  1 . . . -  - - —  (55) 

Y2ttpnv 

The  same  parameters  entering  into  the  computation  of  Y  were  used  to  compute 
V  and  subsequently  p^  using  A(P)  and  A(G).  Figure  8  shows  that  A(P)  and 
A(G)  yield  number  density  values  which  are  significantly  in  excess  of 
those  estimated  via  the  random  walk  sample  space.  If  these  p^  values 
were  used  with  the  PSI  size  distribution  data,  mass  concentration  would 
have  been  significantly  overestimated  for  the  data  shown  in  figure  7. 

It  may  be  argued  that  there  may  be  compensating  errors  in  the  PSI  data, 
i.e.  it  overestimates  y^  anc*  P^  is  underestimated.  However,  comparison  of 
PSI  size  data  with  that  from  other  optical  counters  shows  that  third  moments 
in  the  measured  size  distribution  agree  within  10%.  Compensating  errors 
would  require  that  y^  be  roughly  a  factor  of  4  too  large  which  seems 
unlikely  in  lieu  of  agreement  between  particle  size  analyzers. 

The  same  PSI  system  was  also  used  in  a  series  of  field  measurements 
on  the  same  kind  of  particles  as  those  examined  in  the  laboratory  chamber 
tests.  In  the  field  measurements,  the  PSI  was  placed  within  3m  of  a 
commercially  available  optical  particle  counter  (Cliraet) .  This  device 


dilutes  a  sample  of  the  aerosol  in  order  to  count  individual  particles. 
It  is  claimed  that  by  dilution,  this  device  can  determine  particle 
number  densities  approaching  lO^cc  The  field  tests  involved 
measurements  of  spatially  inhomogeneous  clouds  over  times  of  5  -  10 
minutes.  Figure  9  shows  typical  results  of  the  two  independent 
number  density  determinations  in  14  separate  test  trials.  A  linear 
regression  of  these  data  shows  that  the  two  measurements  are  related  by 

Pn(CLIMET)  =  0.2pn(PSI)  +  2.48  *104  (56) 

with  a  correlation  coefficient  of  0.58.  The  slope  obtained  for  the 
linear  regression  is  in  surprisingly  good  agreement  with  the  fact  that 
80  -  90%  of  the  PSI  measured  size  distribution  was  beyond  the  size  range 
capability  of  the  Climet.  In  view  of  the  fact  that  the  two  devices 
1)  examined  a  spatially  nonhomogeneous  cloud  with  a  separation  of  about 
3m,  2)  covered  different  size  ranges,  and  3)  one  device  required  aerosol 
dilution,  agreement  between  the  two  sets  of  data  is  surprisingly  good. 

The  third  comparison  is  actually  an  examination  of  the  internal 
consistency  of  a  set  of  measurements.  These  data  were  obtained  in  a 
smelter  operating  at  2760°  C.  The  particles  measured  are  believed  to 
be  slag  particles  from  the  walls  of  the  smelter.  The  number  density  for 
these  measurements  has  been  estimated  using  the  random  walk  sample  space 
acceptance  ratio  and  from  the  equation 

PN  ’  N/?apvAt  (57> 

va  At  will  be  called  the  standard  volume. 

-  pv 


N  is  assumed  to  be  the  number  of  measurements  satisfying  the 


requirement  that  the  signal  have  a  minimum  number  of  cycles.  V  is 

the  average  speed  measured  for  the  ensemble  of  measurements  and  At 

is  the  ensemble  acquisition  time.  Water  cooled  probes  were  used  to 

limit  the  probe  volume  length  to  approximately  2.5  cm.  The  fringe 

period  of  the  system  was  36.6  micrometres  and  signal  observation 

with  an  oscilloscope  showed  that  the  maximum  number  of  cycles  in 

the  probe  volume  was  approximately  12.  The  signal  processing 

electronics  required  the  signal  to  have  a  minimum  of  6  cycles. 

Assuming  a  rectangular  cross-section,  the  maximum  value  of  a^  is 

-2  2 

computed  to  be  approximately  9.3  •  10  cm  .  For  the  optical  geometry 
in  these  experiments,  K  =  1,  and  equation  46  reduces  to 


A  =  e 


-2Y 


(58) 


Using  equation  47  and  solving  equation  58  for  p^,  there  results 

In3 (A-1)  (59) 

A  value  of  8  (as  determined  by  oscilloscope  observation)  was  assumed  for  N 
Table  I  lists  the  experimental  values  used  in  the  computation  of  using 

equations  59  and  47.  Note  that  by  comparison  with  the  first  set  of  experi 
ments,  the  A  values  are  quite  large.  The  number  of  particles  producing 
signals  with  too  few  cycles  to  process  was  typically  25%  more  than  N. 

Comparison  of  p^  computed  from  equation  57  with  that  computed  using 
equation  .59  shows  that  when  A  >  0,88  agreement  between  the  two  estimates 


differs  by  no  more  than  about  15%,  In  these  cases,  N  was  very  close  to  all 

particles  detected,  the  turbulence  intensity  (velocity  distribution  standard 

deviation/mean  velocity)  was  about  20%,  the  observed  size  range  was  relatively 

narrow,  and  relatively  few  signals  were  due  to  multiple  particles.  As  A 

decreases,  the  observed  size  distribution  becomes  broader  (with  several 

predominant  mode  sizes)  and  At  becomes  less  representative  of  the  time  required 

to  sample  N.  Hence,  the  greatest  divergences  in  the  estimates  should  occur 

when  the  true  is  the  largest.  The  standard  volume  calculation  shows  roughly 

the  same  particle  concentration  for  all  values  of  A.  Experimentally,  this  is 

known  not  to  be  the  case.  When  the  random  walk  value  of  is  roughly  greater 

than  20  cc  \  the  smelter  was  much  cooler  and  the  observed  particle  concentration 

was  expected  to  be  much  higher  than  for  conditions  where  the  temperature  was 

near  2700°  C.  Finally,  ultraviolet  transmission  measurements  were  found  to  yield 

consistent  number  densities  predicted  using  the  random  walk  sample  space. 

Note  in  examining  the  values  in  Table  I  that  very  small  changes  in  A  correspond 

to  rather  large  changes  in  p^.  This  is  consistent  with  the  plot  for  K  =  1  in 

figure  3.  What  is  surprising  to  this  author  is  how  well  the  model  apparently 

predicts  p._  when  p._  is  small.  It  is  of  interest  to  note  that  when  Poisson 
N  N 

statistics  are  used  to  estimate  p*T  for  the  measured  values  of  A,  that  number 

N 

densities  of  400  -  800  cc  ^  are  predicted  for  the  data.  These  values  were  found 
to  be  inconsistent  with  the  UV  transmission  measurements. 


129 


SUMMARY  AND  CONCLUSIONS 


A  random  walk  analysis  has  been  used  to  predict  how  often  two  or 
more  particles  will  appear  in  a  randomly  varying  sample  volume  of  an 
LV  system.  From  this  model,  the  acceptance  ratio  is  defined  and  applied 
to  the  evaluation  of  fringe  type  LV  system  performance  as  a  function  of 
particle  number  density  and  the  estimation  of  particle  number  density 
when  the  acceptance  ratio  is  measured.  This  model  leads  to  the  following 
conclusions: 

1*  For  coaxial  optical  systems,  the  acceptance  ratio  is  a  function 
of  the  probe  volume  cross-section  geometry  and  is  essentially 
independent  of  the  optical  system  depth-of-f ield . 

2-  It  is  found  that  as  number  density  increases,  the  acceptance 

ratio  decreases.  Thus,  data  rate  for  an  LV  cannot  be  increased 
indefinitely  by  increasing  p  .  For  a  given  optical  system 
geometry,  there  exists  an  optimum  acceptance  ratio  beyond  which 
the  data  rate  will  decrease. 

3.  When  the  number  density  is  constant,  it  is  found  that  the 
acceptance  ratio  increases  as  k,  the  probe  volume  cross-section 
ratio,  decreases. 

4.  When  k  reaches  a  value  of  about  0.5,  the  acceptance  ratio  will 
not  significantly  improve  as  the  number  density  is  increased 

and  k  is  decreased.  The  data  rate  can  be  increased  by  decreasing 
k  as  number  density  increases. 

5.  Comparison  of  number  density  estimates  via  the  random  walk 
acceptance  ratio  with  other  experimental  methods  such  as  optical 


particle  counters  and  mass  concentration  measurements  shows 
that  the  acceptance  ratio  method  yields  reasonable  results. 

6.  The  acceptance  ratio  was  not  found  to  yield  number  densities 
consistent  with  experiment  when  constant  sample  volumes  were 
used  with  Poisson  or  Gaussian  statistics  which  attempt  to 
estimate  the  probability  of  two  or  more  particles  in  the  sample 
volume  on  the  basis  of  mean  numbers  of  particles  in  the  sample 
volume . 

7.  The  acceptance  ratio  method  can  be  used  to  estimate  number 
density  with  conventional  LV  fringe  type  optical  systems  and 
burst  signal  processors  with  modest  modification  and  system 
calibration  to  determine  k. 

The  limitations  and  potential  modifications  of  this  model  should  be 

clearly  understood.  The  model  does  not  directly  account  for  number  density 

as  a  function  of  particle  size.  N  will  ultimately  be  a  function  of  particle 

Li 

size.  Using  an  average  N  determined  experimentally  thus  weights  the  number 
density  estimate  toward  those  sizes  most  frequently  giving  a  measureable 
signal.  In  light  of  experimental  results,  this  appears  to  be  a  reasonable 
approach.  There  is,  however,  considerable  room  for  model  refinement  in 
this  area.  For  example,  this  technique  assumes  all  multiple  particle 
signals  are  detected.  This  is  known  to  not  be  the  case  when  A  is  small 
and  the  model  accuracy  can  be  expected  to  hinge  directly  on  signal  processor 
efficiency  in  detecting  such  signals.  It  should  also  be  realized  that 
when  the  particle  number  density  becomes  so  high  that  optical  transmission 
becomes  limited,  the  LV  system  cannot  function  and  the  model  must  fail. 


Thus,  future  refinements  must  also  include  transmission  effects. 
Experimentally,  it  has  been  determined  that  the  instrument  and  model 
appear  to  work  well  for  transmittances  as  low  as  80%. 
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Figure  1.  Generation  of  free  space  interference  fringes 
with  a  laser  Doppler  velocimeter  system. 
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ACCEPTANCE  RATIO 


Fig.  3.  A  as  a  Function  of  y  With  K  as  a  Parameter. 
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Fig.  4.  y  as  a  Function  of  K  With  A  as  a  Parameter. 


157 


V 


ACCEPTANCE  RATIO 


K-J 


K>4 

K-.6 


K*.8 


X-t 


Fig. 5  Optimum  A  for  Maximum  Data  Rate. 
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FIGURE  6 .  Schematic  of  the  PSI  optical  system  used 
in  the  laboratory  and  field  measurements. 
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MASS  CONCENTRATION  (  gm/m3) 


Fig.  7.  Comparison  of  mass  concentrations  as  determined  by 

filter  weights  and  PSI  measurements  for  red  phosphorus 
and  vaporized  condensed  fog  oil  smokes. 


NUMBER  DENSITY  (ccM) 


Fig.  8.  Comparison  of  number  densities  for  the  smoke  data  in 
Fig.  1 .  as  estimated  using  three  statistical  models. 
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Fig*  9 •  Comparison  of  number  densities  as  determined  by  a  single 

particle  counter  and  that  using  a  random  walk  sample  space 
model. 
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OPTICAL  PARTICLE  SIZE  MEASUREMENTS 


OF  HYGROSCOPIC  SMOKES  IN 
LABORATORY  AND  FIELD  ENVIRONMENTS 

W.  M.  Fanner  ,  R.  D.  Morris  ,  F.  A.  Schwartz 


ABSTRACT 

Optical  Particle  Size  Analyzer  (PSA)  measurements  of  hygroscopic 
smokes  are  reported  from  experiments  conducted  in  laboratory  chambers 
and  in  field  tests*  Three  different:  instruments  were  used:  a  particle 
sizing  interferometer,  a  Climet  (model  211),  and  a  Particle  Measurements 
System  CSASP-100-HC .  Field  measurements  were  made  under  relative 
humidity  conditions  ranging  between  65  and  97%.  The  results  indicate 
that  the  size  distributions  are  at  least  bimodal  -  not  log-normal  - 
and  that  mean  particle  diameter  does  not  appear  to  increase  with 
increasing  relative  humidity. 

INTRODUCTION 

Hygroscopic  smokes  form  a  major  portion  of  the  conventional 
obscurants  inventory  of  the  U.  S.  Army.  Additionally,  there  exist 
numerous  other  particle  sources  which  are  thought  to  be  hygroscopic. 
Changes  in  the  shape  of  the  particle  size  distribution  as  relative 
humidity  increases  is  of  fundamental  importance  in  describing  the  optical 
properties  of  these  smokes.  For  example,  predictive  models  hypothesize 
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that  the  particle  sizes  in  the  distribution  increase  as  relative 
humidity  increases.  As  a  result,  optical  transmission  through  such 
smokes  may  significantly  decrease  over  a  broad  range  of  wavelengths. 

In  an  attempt  to  produce  a  solid  base  of  field  test  data  on  the 
performance  of  conventional  hygroscopic  smokes  at  RH  values  greater 
than  60%  and  to  begin  resolution  of  data  differences,  the  U.  S.  Army 
PM  Smoke  Office  conducted  the  High  Humidity  Hygroscopic  Smoke  (H3S) 
field  test  at  the  Edgewood  Area  of  Aberdeen  Proving  Ground  during 
July  1979. 

3 

During  the  H  S  test,  particle  size  distribution  data  were  obtained 
under  conditions  which  allowed  the  results  from  the  PSA  systems  to  be 
reasonably  compared  (i.e.  close  spatial  positions  and  similar  size 
ranges).  These  particle  sizing  systems  were  1)  a  particle  sizing 
interferometer  operated  by  the  University  of  Tennessee  Space  Institute 
and  2)  a  Climet  operated  by  Dugway  Proving  Grounds  (DPG) . 

The  initial  intent  was  to  compare  the  PSA  systems  in  a  small 
laboratory  test  chamber  immediately  before  they  were  moved  to  the  H  S 
test.  In  this  way  specific  differences  in  PSA  data  under  controlled 
conditions  could  be  identified  and  data  could  be  obtained  by  the  same 
instruments  in  both  laboratory  and  field  environments.  Because  of 
shipping  difficulties,  one  of  the  PSA  systems  (a  Climet  operated  by  DPG) 
did  not  arrive  in  time  for  the  laboratory  comparison.  However,  the  second  PSA 
system  (a  Particle  Sizing  Interferometer  operated  by  UTSI)  was  used  in  both 
the  field  test  and  in  a  laboratory  comparison  with  a  Particle  Measuring 
Systems  PSA  (PMS  CSASP—100— HC)  operated  by  the  Research  Division  of  the 
Chemical  Systems  Laboratory.  Direct  comparison  of  laboratory  and  field  test 
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results  can  thus  be  made  only  for  data  obtained  with  the  particle 
sizing  interferometer. 

In  reporting  these  results,  the  objectives  of  this  paper  will 
be  to: 

1.  Compare  hygroscopic  particle  size  distribution  data  obtained 
from  three  (3)  different  optical  PSA  systems  operating  in 
laboratory  and/or  field  conditions. 

2.  Present  results  for  the  measurement  of  size  distributions  °f 
hygroscopic  smokes  under  conditions  of  high  RH. 

3.  Test  the  hypothesis  that  particle  size  increases  with 
increasing  RH. 

4.  Determine  how  well  PSA  concentxatiion  estimates  compare  with 
chemical  impinger  data. 

To  achieve  these  objectives,  the  pretest  PSA  comparison  and 
subsequent  results  are  first  described.  Next,  PSA  results  obtained 
during  H  S  and  examples  of  size  distributions  as  measured  by  the  PSA 
systems  are  given.  The  measured  effects  of  RH  on  mean  particle  size  are 
shown  and  concentration  estimates  obtained  with  the  particle  sizing 
interferometer  are  compared  with  those  obtained  with  a  nearby  chemical 
impinger. 

PRETEST  PSA  COMPARISON 

The  optical  geometry  for  the  PSI  is  shown  in  Figure  1.  The  system 
was  arranged  to  measure  a  particle  size  range  of  less  than  0.3  -  greater 
than  6.0  micrometres.  The  smallest  particle  size  detectable  with  the 
system  was  estimated  to  be  approximately  0.2  micrometres  for  water  droplets. 
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This  PSI  used  a  5  milliwatt  HeNe  laser,  and  F/2.2  receiving  lens. 

A  detailed  description  of  the  theory  and  operation  of  the  PSI  system 
is  given  in  references  2-5.  The  PSI  functions  on  the  principle  that 
particle  sizes  can  be  determined  by  measurement  of  a  signal  shape 
parameter  from  light  scattered  out  of  a  localized  fringe  pattern  called 
the  probe  volume.  Both  particle  size  and  velocity  data  are  obtained 
with  a  PSI.  Data  obtained  with  the  PSI  signal  processor  and  data 
acquisition  system  in  these  experiments  were  recorded  on  magnetic 
disks  for  hard  copy  presentation  and  data  storage.  The  PSI  measured 
particle  velocity,  signal  magnitude  (which  is  also  a  measure  of  particle 
size),  the  PSI  size  parameter-visibility,  and  a  particle  number  density 
parameter.  The  particle  size  histogram  was  divided  into  21  bins  corresponding 
to  0.3  micrometre  size  increments  ranging  from  all  detectable  particles  less 
than  0.3  to  all  greater  than  6.0  micrometres. 

The  PMS  CSASP-100-HC  uses  an  intake  horn  to  shape  an  air  flow  drawn 
through  the  instrument  with  a  small  fan.  The  device  uses  a  5  milliwatt 
HeNe  laser  and  measures  the  scattered  signal  magnitude  observed  with  a 
lens  in  the  forward  scatter  direction.  The  signal  magnitudes  (corresponding 
to  equivalent  latex  sphere  diameters)  are  sorted  into  histogram  bins 
beginning  at  0.3  micrometres  and  increasing  in  0.04  micrometre  increments 
to  1.0  micrometre  where  the  bin  width  increases  to  0.1  micrometre.  The 
last  two  histogram  bins  cover  1.8  -  2.1  and  2.1  -  2.4  micrometres.  For 
this  experiment  the  histogram  bin  counts  were  read  and  recorded  manually 
from  a  panel  display  on  the  PMS  data  acquisition  and  signal  processing  system. 

The . laboratory  measurements  were  conducted  in  a  sealed  box  with  an 
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atmosphere  initially  at  ambient  temperature  and  humidity  conditions. 

Figure  2  schematically  illustrates  the  experimental  arrangement  for  the 
measurements.  The  PSI  system  was  placed  directly  inside  the  chamber. 

Filter  samples  and  samples  obtained  with  an  Andersen  Cascade  Impactor 
were  obtained  via  samples  drawn  from  wall  ports.  Since  the  PMS  particle 
analyzer  measurements  become  uncertain  due  to  the  probability  of  multiple 
particle  signals  when  particle  number  density  is  high,  it  was  necessary 
to  first  dilute  the  smoke  by  allowing  it  to  enter  a  dilution  chamber 
through  a  small  hole  (diameter  approximately  3  mm)  in  the  side  of  the 
main  chamber. 

Phosphorus  and  hexacloroethane  (HC)  smokes  were  generated  by  igniting 
small  samples  at  one  end  of  the  chamber*  Fog  oil,  while  not  hygroscopic, 
was  easy  to  generate  using  a  small  commercially  available  fog  oil 
generator  and  it  also  was  measured  using  both  PSA  systems.  After  initial 
smoke  generation,  the  smoke  was  given  time  to  become  evenly  dispersed  in 
the  chamber.  During  this  time  the  operating  parameters  for  the  PSI  were 
set  and  the  other  instruments  tested  for  operational  readiness.  A  common 
clock  was  used  to  synchronize  time  for  all  instrument  samples.  Transmisso- 
meter  recordings  were  made  continuously  on  strip  charts.  As  soon  as  filter 
samples  were  obtained  they  were  weighed  to  minimize  any  evaporation  losses. 
The  PSI  data  acquisition  system  was  instructed  to  obtain  chronological 
histograms  of  10  counts  each  and  to  record  all  data  on  magnetic  disks. 

The  PMS  system  was  manually  interrogated  at  approximately  10  minute 
intervals.  A  typical  test  might  run  30  to  60  minutes,  after  which  the 
chamber  was  exhausted  in  preparation  for  the  next  test. 
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Transmis some ter  measurements  were  used  to  compute  an  optical  depth 
as  a  function  of  time  after  smoke  initiation.  Mass  concentration  as  a 
function  of  time  was  obtained  using  filter  samples  which  trapped  a  portion 
of  the  smoke  withdrawn  from  the  chamber  at  a  fixed  flowrate.  Since  it 
has  been  assumed  that  these  smokes  follow  log-normal  distributions,  the 
geometric  mean  and  standard  deviations  for  the  PSA  data  were  estimated. 

The  PMS  data  were  plotted  on  log-probability  paper  to  obtain  the  geometric 
mean  diameter  and  logarithmic  geometric  standard  deviation  and  to  determine 
how  closely  the  distribution  approximated  a  log-normal  distribution. 

PSI  data  were  used  to  compute  the  probability  density,  geometric  mean 
diameter,  geometric  standard  deviation  and  particle  number  density.  The 
method  of  estimating  particle  number  density  is  new  and  described  in  a 
recent  UTSI  report.^  By  using  the  PSI  measured  distribution  functions, 
a  mean  particle  volume  was  computed  which  together  with  the  particle 
number  density  estimates  provided  an  estimate  of  the  mass  concentration. 

EXPERIMENTAL  RESULTS  FROM  LABORATORY  MEASUREMENTS 

Table  I  shows  a  comparison  of  typical  data  for  the  three  smokes 
examined:  HC,  red  phosphorus,  and  fog  oil.  Geometric  mean  diameter, 

geometric  standard  deviation  and  third  moments  are  compared  for  the  particle 
size  distributions  determined  by  the  PSI  and  the  PMS  CSASP.  It  should  be 
borne  in  mind  when  comparing  these  data  that  the  PSI  counted  and  recorded 
particles  less  than  0.3  micrometres  in  diameter  while  the  PMS  system  began 
at  0.3  micrometres  and  the  smoke  entered  a  dilution  chamber  through  a 
relatively  small  hole  prior  to  measurement.  There  is  surprisingly  good 
agreement  between  the  two  measurements  when  the  geometric  means  are  compared. 
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The  agreement  is  deceptive  however.  The  PSI  measurements  include  a 
large  number  of  particles  less  than  0.3  micrometres  in  diameter.  This 
becomes  apparent  when  the  logarithmic  geometric  standard  deviations  are 
compared.  The  PSI  values  are  typically  20  -  30%  greater  than  those 
measured  by  the  PMS  system.  This  results  from  the  fact  that  the  PMS 
system  measured  virtually  no  particles  greater  than  1.0  micrometres  in 
diameter.  The  PSI,  on  the  other  hand,  typically  found  3  -  4%  of  the 
distribution  was  greater  than  1.0  micrometre.  This  difference  is 
significant  when  the  third  moments  of  the  distribution  are  computed  in 
order  to  estimate  mass  concentration.  The  PSI  third  moments  are  roughly 
a  factor  of  10  greater  than  the  PMS  values.  A  number  density  estimate 
could  not  be  obtained  from  the  PMS  system  since  the  smokes  were  diluted 
before  measurement.  The  number  density  column  in  Table  I  is  the  estimate 
obtained  with  the  PSI.  In  order  to  indicate  the  fraction  of  particles 
less  than  0.3  micrometres  measured  by  the  PSI,  a  column  labeled  %<_0.3  \xm 
provides  the  probability  density  amplitudes  for  this  size  fraction.  Figure 
3  shows  a  comparison  of  concentration  estimates  obtained  mechanically  by 
weight  and  that  obtained  using  PSI  data.  A  major  uncertainty  in  the 
comparison  is  the  material  density  of  the  smokes  measured.  Orthophosphoric 
acid  (density  =  1.28  gm/cc  at  82%  RH)  is  assumed  for  red  phosphorus  and  for 
fog  oil  a  density  of  0.9  gm/cc  is  assumed.  Agreement  between  the  filter 
and  the  PSI  estimates  is  remarkably  good  for  the  red  phosphorus  and  fog 
oil  when  the  large  differences  in  the  two  methods  to  arrive  at  the 
concentrations  are  considered.  The  PSI  technique  is  purely  optical  and 
subject  to  any  localized  variations  in  number  density  which  may  exist  in 
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the  test  chamber.  The  filter  method  of  determining  the  concentration 
is  purely  mechanical  -  the  smoke  is  pulled  to  the  filter  at  a  fixed 
flowrate  for  a  fixed  time  interval  and  weighed. 

Figures  4  through  6  plot  examples  of  the  particle  size  distributions 
obtained  with  the  PMS  system.  The  deviations  from  the  straight  line  fits 
to  the  data  occur  where  the  particle  counts  are  relatively  small  and 
potential  uncertainties  large.  As  the  figures  show,  the  PMS  size 
distributions  are  represented  reasonably  well  by  log-normal  distributions 
with  mean  diameters  and  logarithmic  standard  deviations  indicated  in  the 
figures. 

Figure  7  shows  representative  probability  density  functions  as 
measured  by  the  PSI  for  the  smokes  and  for  the  background  aerosol  present 
in  the  laboratory  where  the  measurements  were  made.  The  data  clearly 
show  the  contribution  of  the  smokes  to  the  small  particle  background. 

The  strong  multimodal  characteristic  of  the  fog  oil  smoke  is  a  result 
consistently  seen  in  a  large  number  of  different  tests. 

RESULTS  OF  PSA  FIELD  TEST  MEASUREMENTS 

After  the  pretest  laboratory  calibration,  the  PSI  was  moved  to  position 
3 

in  the  H  S  field  test.  It  was  placed  approximately  2  m  off  grid  center, 

1  m  off  ground  level,  and  3  m  from  the  Clime t  PSA  and  a  chemical  impinger. 
Actually  operated  were  two  Climet  PSA  systems.  The  system  data  which  will 
be  reported  here  covered  a  size  range  of  0.3  to  3.0  micrometres  in  six 
unequal  size  increments.  The  second  Climet  system  covered  a  range  of 
0.3  to  13  micrometres  in  six  unequal  increments.  Typically,  the  second 
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system  obtained  measurements  in  only  the  three  smallest  size  intervals 
(the  largest  size  measured  being  5,2  micrometres).  This  data  was 
considered  to  be  too  poor  in  size  resolution  for  the  comparison  with 
PSI  data,  although  this  data  does  provide  some  insight  into  the  Climet 
and  PSI  data  comparisons. 

Unlike  the  PSI,  the  Climet  PSA  determines  particle  size  by  measuring 
the  pulse  magnitude  of  white  light  scattered  in  the  forward  direction  by 
single  particles  withdrawn  from  the  smoke  cloud,  diluted  in  number  density 
and  blown  across  the  incident  light  beam  at  a  fixed  flowrate.  This  PSA 
is  calibrated  with  latex  spheres.  As  a  result,  the  particle  size  measured 
by  the  Climet  must  be  interpreted  as  an  equivalent  latex  sphere  diameter. 
Calculations  given  in  reference  show  t:hat  when  the  particles  are 
nonabsorbing  dielectrics,  the  response  function  of  this  instrument  is 
mono tonic. 

The  results  of  the  particle  size  distribution  measurements  as  a 
function  of  RH  will  be  presented  as  1)  sample  moments  of  the  size  distri¬ 
bution,  2)  mean  particle  diameters  (geometric,  volumetric  and  mass  means) 
and  3)  histograms  of  the  size  distribution  for  selected  trials. 

The  means  and  moments  of  the  particle  size  distribution  are  defined  as 
follows.  Let  the  data  be  grouped  into  n  intervals  with  a  frequency  of 
occurrence  in  the  ith  interval  given  by  f^.  The  total  number  of  measurements 
is  then  given  by 

N  -  If.  (1) 

i=l 
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The  kth  sample  moment  of  the  size  distribution,  is  then  given  by 


1  n  k 
T7=-  Z  f.D? 

^t  i=l  1  1 


(2) 


where  D  is  defined  as  the  largest  diameter  value  in  the  size  increment 
i 

Using  these  definitions  the  following  mean  diameters  can  be  written  as: 


1 .  Numeric  Mean  Diameter 


Di  ■  “i 


(3) 


2.  Volumetric  Mean  Diameter 


»,  -  *3/3 


(4) 


3.  Mass  Mean  Diameter 


D  =  U./U, 

mm  4  3 


(5) 


A  log-normal  size  distribution  is  often  assumed  for  analytical  calculations. 
This  distribution  is  characterized  by  the  geometric  mean  diameter  and  the 
logarithmic  geometric  standard  deviation  0  •  These  are  defined  by 


•  D  =  EXP  7j  Z  f .  InD . 

8  Nt  i=i  1  1 


(6) 


n 


0  =  EXP  tIt  z  f.ln  (D±/D  ) 

S  Nt  i=l  1  1  S 


(7) 
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If  the  size  distribution  can  be  approximated  as  a  Gaussian 


distribution,  then  the  second  central  moment,  m^,  (often  called  the 
"variance")  is  used  as  a  measure  of  the  "spread"  of  the  distribution. 


^2  is  defined  as 


i  " 


2 


fi(Di  -  V 


(8) 


The  measurements  reported  here  are  for  phosphorus  smokes  which 
were  dispersed  as  either  red  phosphorus  grenades  or  white  phosphorus 
wicks  or  wedges. 

Tables  II  and  III  summarize  the  computations  for  the  means  and 
moments  for  all  PSA  size  measurements  obtained  during  a  specified  trial. 
Table  II  additionally  tabulates  the  proportion  of  the  total  particles 
detected  by  the  PSI  which  were  less  than  0.3  micrometres  in  diameter. 

These  particles,  while  too  small  to  resolve  for  measurement,  could  be 
detected  and  were  included  in  the  total  numeric  count. 

Figures  8-12  compare  size  distribution  histograms  obtained  by  the 
PSA  systems  for  examples  of  data  for  each  of  the  different  kinds  of 
munitions.  The  PSI  data  in  these  figures  have  had  the  smallest  size  bin 
(D  <0.3  pm)  artificially  truncated  in  order  to  get  a  more  direct  comparison 
with  the  Climet  data.  The  PSI  data  show  a  reasonable  consistency  in  size 
distribution  shape  for  all  phosphorus  trials.  These  data  show  that  the 
distribution  is  at  least  bimodal  with  one  mode  occurring  at  less  than  0.6 
micrometres  and  the  second  occurring  at  around  3.0  micrometres  in  diameter. 
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The  HC  data  show  a  multimodal  distribution  with  a  relatively  large 
number  being  less  than  0.6  micrometres.  The  Climet  data  generally 
shows  a  strong  size  mode  for  sizes  less  than  0.6  micrometres  and  a 
second  mode  which  appears  between  1.0  and  1.5  micrometres.  The  two 
instruments  are  therefore  in  agreement  that  the  predominant  size  mode 
occurs  below  0.6  micrometres.  However,  the  PSI  data  indicates  that 
the  larger  size  mode  should  occur  at  about  3.0  micrometres  and  that 
significant  numbers  of  particle  sizes  exist  in  the  3-5  micrometre 
range.  The  later  results  also  seem  to  be  verified  by  the  second  PSA 
which  detected  particles  in  the  1-5.2  micrometre  range  (see  figure  8). 

As  a  result  of  the  Climet  not  measuring  a  second  mode  with  sizes  as  large 
as  that  observed  with  the  PSI,  Tables  II  and  III  show  that  significant 
differences  exist  in  the  computed  moments  for  the  two  systems. 

EFFECTS  OF  RELATIVE  HUMIDITY  ON  MEAN  PARTICLE  SIZE 

An  assumption  made  in  the  analysis  of  hygroscopic  smokes  is  that  as 
the  relative  humidity  increases  the  mean  particle  diameter  will  increase 
significantly.1  It  is  of  interest,  therefore,  to  examine  the  PSA  results 
for  the  effects  of  RH  on  mean  particle  diameter.  Figures  13  and  14  plot 
mean  particle  diameter  as  a  function  of  RH  as  measured  by  the  two  PSA 
systems.  The  PSA  data  do  not  show  an  increase  in  mean  particle  diameter 
with  increasing  RH.  In  fact,  the  data  show  a  slight  volumetric  and 
geometric  mean  diameter  decrease  with  increasing  RH.  Only  the  PSI  mass 
mean  appears  to  increase  slightly.  The  geometric  mean  yields  an  average 
value  close  to  the  most  frequently  occurring  size  in  the  distribution. 
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The  strongest  PSI  mode  in  all  the  data  occurred  for  sizes  less  than 
0.3  micrometres.  Therefore,  the  PSI  geometric  mean  is  smaller  than  that 
for  the  Climet  which  begins  counting  only  for  sizes  greater  than  0.3 
micrometres.  The  differences  in  the  mass  and  volumetric  means  are  due 
primarily  to  relative  values  of  the  second  mode  in  the  size  distributions. 

As  a  result,  the  PSI  values  for  volumetric  and  mass  mean  diameters  are 
roughly  a  factor  of  2  greater  than  those  measured  by  the  Climet.  It  is 
interesting  to  note  that  in  the  case  of  the  mass  mean  measurements,  the 
PSA  data  is  a  factor  of  2  to  4  times  greater  than  that  assumed  in  the 
analytical  computations.^  This  assumption  was  based  on  data  obtained 
with  an  Andersen  cascade  irapactor  measuring  smokes  in  a  laboratory  test 
cell.  However,  if  a  log-normal  distribution  is  assumed,  then  the  geometric 
mean  and  standard  deviation  values  yield  mass  mean  diameter  values  which 
are  very  near  those  assumed  for  the  analytical  calculations. 

The  apparent  decrease  in  mean  particle  diameter  with  increasing 
humidity  is  unexpected.  However,  the  PSI  data  may  provide  a  clue  in 
explaining  this  anomalous  behavior.  Figure  15  plots  the  detected  particles 
smaller  than  0.3  micrometres  in  diameter  as  a  function  of  humidity.  The 
trend  shown  by  this  figure  is  for  the  number  of  small  particles  detected 
by  the  PSI  is  to  increase  with  increasing  relative  humidity.  Also  shown  is 
the  equation  obtained  for  a  logarithmic  curve  fit  with  a  correlation 
coefficient  of  0.68.  From  this  data  it  is  hypothesized  that  very  small 
particles  released  by  the  munition  source  (which  otherwise  could  not  be 
detected  by  the  PSI  and  are  evidently  always  present  in  the  size 
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distribution)  grow  to  measurable  size  as  the  humidity  increases. 

As  the  numbers  of  detectable  small  particles  increase  with 
increasing  humidity,  there  is  a  shift  in  the  distribution  which 
apparently  balances  any  detectable  size  increase  in  the  larger  sizes. 

If  a  PSA  system  were  used  which  could  respond  to  the  small  size  tail 
of  the  distribution  then  Figure  15  suggests  that  particle  growth 
would  indeed  be  observed, 

PSA  CONCENTRATION  ESTIMATES 

The  PSA  systems  produce  data  which  can  be  used  with  the  obscurant’s 
material  density  to  obtain  an  estimate  of  mass  concentration.  Pretest 
comparison  indicated  that  the  PSI  could  obtain  concentration  estimates 
which  were  in  reasonable  agreement  with  weighed  samples  withdrawn  from 
the  test  chamber  (see  Figure  3) ,  It  is  of  interest  therefore  to  determine 
how  well  PSI  concentration  estimates  agree  with  those  obtained  by  chemical 
impingers .  It  should  be  expected  that  the  PSI  data  should  yield  concen¬ 
tration  values  which  are  larger  than  the  dosage  values  obtained  with 
the  chemical  impingers.  To  obtain  concentration  from  a  chemical  impinger, 
a  yield  factor  multiplies  the  dosage  value  to  account  for  water  absorbed 
by  the  elemental  phosphorus.  By  ratioing  the  PSI  concentration  and  the 
chemical  impinger  dosage  an  ’’experimental"  yield  factor  for  the  field 
data  is  obtained. 

To  compare  PSI  concentration  estimates  with  dosage  estimates  it  is 
necessary  to  compute  the  trial  time  average  PSI  concentration.  This  is 
given  by 

7T 

<  C(PSI)>  =  -  <  U~N >  (9) 

6  5 
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where  <  >  indicates  a  time  average,  C(PSI)  is  the  concentration  as 
determined  from  PSI  data,  y^  t^le  third  sample  moment  and  N  is  the 
corresponding  number  density.  Number  density  data  from  the  Climet  PSA 
was  not  available  for  this  calculation.  Time  averaged  dosage  data 
were  obtained  for  the  chemical  impinger  located  nearest  the  PSI 
during  the  field  test.  The  values  for  p  were  computed  using  equation 

*  V 

25  in  reference  1.  Values  of  <  C(PSI)>  and  <D  (CI)>,  the  time  average 
dosage  value  for  the  chemical  impinger  nearest  the  PSI,  are  listed 
in  Table  IV  for  the  phosphorus  trials.  Table  IV  lists  the  theoretical 
yield  factor  used  to  compute  the  concentration  from  the  dosage.  An 
experimental  yield  factor  YF^,  is  listed  and  is  computed  from 

<  C(PSI)> 

YF  =  - t -  (10) 

e  <  D  (CI)> 

Comparison  of  the  experimental  with  the  theoretical  yield  factor  shows 
large  differences  with  a  few  points  showing  close  agreement.  In  general, 
the  dosage  values  are  much  closer  to  the  PSA  determined  concentration 
values  than  the  values  obtained  by  using  the  theoretical  yield  factors. 

CONCLUSIONS 

Examination  of  the  data  shown  in  the  previous  sections  leads  to 
the  following  conclusions. 

1.  Size  distributions  measured  by  the  PSA’s  are  usually  bimodal. 
One  mode  occurs  for  sizes  less  than  0.6  micrometres.  The  second  mode 
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appears  for  sizes  between  1.0  -  1.5  micrometres  or  at  about  3.0 
micrometres  depending  on  which  PSA  data  is  used. 

2.  Size  distributions  measured  in  the  laboratory  and  field  by 
the  PSI  have  similar  mode  characteristics. 

3.  Although  the  measured  size  distributions  are  not  log-normal, 
mass  mean  diameters  computed  from  geometric  mean  and  standard  deviations 
assuming  a  log-normal  distribution  are  close  to  previously  assumed  values. 

4.  PSI  size  distribution  data  indicates  that  most  of  the  detected 
particles  are  less  than  0.3  micrometres  in  diameter. 

5.  The  numbers  of  particles  less  than  0.3  micrometres  in  diameter 
appear  to  increase  as  relative  humidity  increases. 

6.  Measured  mean  particle  diameters  (geometric,  volumetric,  or 
mass)  do  not  appear  to  increase  as  relative  humidity  increases. 

7.  Measured  mass  means  are  a  factor  of  2  to  4  times  greater  than 
those  assumed  in  a  previous  analytical  computation. 

8.  Concentration  estimates  made  with  PSI  data  are  generally 

smaller  than  those  obtained  using  a  theoretical  yield  factor  and  dosage  data. 
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Figure  1.  Schematic  of  the  PSI  optical  system  used  in  the  laboratory  and  field  measurements. 


Figure  2.  Schematic  of  the  smoke  chamber  for  laboratory  comparisons. 
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Figure  5.  Pretest  comparison  of  mass  concentration  estimates  obtained 
using  filter  sampling  and  PSI  measurements. 
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Figure  6.  PMS  measured  size  distribution  for  fog  oil. 
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determined  by  the  PS  I. 
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PARTICLE  SIZING  INTERFEROMETER 
PARTICLE  SIZE  DISTRIBUTION  DATA 
OBTAINED  DURING  THE  H3S  TEST 
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Volumetric  mean  diameter  evereged  over  the  trial. 

Kaaa  mean  dlemeter. 

Second  central  moment  or  variance. 

Geometric  standard  deviation  averaged  over  the  trial 
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NUMERICAL  CALCULATIONS  OF  VISIBILITY  USING  MIE  THEORY 
OR  REFRACTION  THEORY  OF  GEOMETRICAL  OPTICS 

Mie  theory  numerical  calculations  of  the  visibility  parameter  for 
the  interference  fringe  Laser  Doppler  Velocimeter  (LDV)  have  been  made 
by  Adrian  and  Early, *  Adrian  and  Orloff^,  Chu  and  Robinson^,  Hong  and 
Jones11,  and  others.  These  calculations  have  typically  considered  only 
a  paraxial  forward  or  back  scatter  configuration.  An  exception  was 
the  work  of  Hong  and  Jones11  which  achieved  a  more  flexible  program  which 
allowed  the  calculation  of  visibility  with  the  photomultiplier  tube  (PMT) 
located  in  any  direction  from  the  probe  volume.  However,  the  published 
calculations  of  Hong  and  Jones  were  only  for  small  rectangular  apertures. 

The  approach  of  Hong  and  Jones  has  been  extended  in  this  research  to  de¬ 
velop  a  computer  code  for  a  more  realistic  configuration:  a  circular 
aperture  of  the  sizes  used  in  actual  experimental  work.  Calculations  of 
visibility  may  be  made  with  or  without  beam  stops  with  the  PMT  aperture  in 
any  direction  from  the  probe  volume.  For  aperture  locations  and  particle 
sizes  where  refraction  theory  is  adequate,  refraction  theory  may  replace 
the  Mie  theory  subroutine  by  following  the  work  of  Bachalo.^ 

Figure  1  shows  the  scattering  geometry  of  the  two  laser  beams  with 
the  z  axis  defined  as  the  bisector  of  the  angle  between  the  beams.  The 
y-axis  is  in  the  plane  of  the  beams  aligned  normal  to  the  z  axis.  Polarization 
of  the  incident  beams  is  in  the x direction . Associated  with  the  (x,y,z)  system 
is  an  (r,  6,  <J>)  system  as  shown.  It  is  assumed  with  that  all  light  which  is 
scattered  into  the  solid  angle  of  the  PMT  aperture  is  delivered  to  the  PMT. 


Following  this  assumption,  the  intensity  resulting  from  the  sum  of  scattered 
fields  is  calculated  in  the  plane  of  the  aperture,  and  integration  of  this 
intensity  over  the  aperture  determine  radiant  power  delivered  to  the  PMT. 

The  vector  _r  locates  the  field  points  in  the  aperture  plane. 

The  scattering  process  of  the  LDV  was  discussed  briefly  by  Chu  and 
Robinson3.  Adapted  versions  of  their  equations  (1)  through  (12)  are  given 
below  as  background  to  the  development  of  an  explicit  expression  for  the 
visibility  which  may  be  evaluated  by  numerical  integrations.  Far  from 
the  scattering  region  at  position  _r  and  time  _t  the  electric  field  scattered 
by  the  particle  is 

Jil, 2  (£>  t)  =  E1>2  (r)  exp  [i  (w  +  t]  (1) 


where  the  subscript  1  or  2  indicates  one  of  the  two  incident  laser  beams 
as  shown  in  Figure  1.  The  Doppler  frequency  shift  is  2»  Born  and 

Wolf^  discuss  an  averaging  procedure  for  the  Poynting  flux  which  results 
in  the  expression 

S  <r,t)  =  h.  Re  ({  h  <r,t)  +  E2  (r,t)  }  X  {Hx  (r,t)  +  H2  (r,  t)  }*)  (2) 

where  2  Cr>t)  are  the  magnetic  fields  scattered  by  the  particle  from  beams 
1  and  2.  Defining 

.S  ~  £a  x  Sg  a,  3  =  1,2  (*  indicates  complex  conjugation)  (3) 

allows  (2)  to  be  written  as 

S  (_r,t)  =  \  Re  (  S^n  (r)  +  $22  (£,)  +  S.12  U)  exP  +S2l(r)  exp  (-iMjt))  (4) 

where  -  A'u^  -  is  the  Doppler  modulation  frequency. 

Use  of  the  subscripts  a  and  3  will  indicate  the  quantity  in  question  to  be 
a  function  of  position  _r  only. 


b2 


3 


The  energy  per  unit  time  that  is  scattered  into  the  solid  angle  of 
the  PMT  aperture  is  given  by 

W(t)  =  Jr2  dfi  e  r  .•  £  (r,t)  (5) 

where  the  integral  is  taken  over  the  solid  angle  of  the  aperture  and 
where  er  is  the  radial  unit  vector.  Substituting  E  q.  (4)  into  Eq.  (5) 


gives 

W(t)  =  %  Re(wll  +  w22  +  w12  exP  U  +  w21  exp  (-ifc^t))  (6) 

with  W  defined  as 

ag 

Wag  =  ^r2  dn  er  (7) 

Eq.  (6)  may  then  be  expressed  as 

W(t)  =  W  ,  +W  cos  W.t  +  4>)  (8) 

dc  ac  d 

where  Wdc  =  h  Re  (Wn  +  W22)  (9) 

Wac  “  h  (  {  Re  («12  +  W21)}  2  +  {1^12  -W21> }  2)  ^  (10) 

$  =  tan  '  Im  ~W21^  (11) 

Re  (W12  +  W21) 


The  visibility  function  V  is  defined  as 

V  =  Wac 
Wdc 

For  a  plane  wave  propagating  in  a  medium  with  intrinsic  impedance  X- 

Hg  =  (er  X  Eg)/Z 

where 


(12a) 


\  > 

Z  =  (op/ (we  +  i  o)J 


(13) 


4 


where  er  is  the  unit  vector  in  the  direction  of  propagation. 

u 

For  a  medium  with  electrical  conductivity  o  =  0,  7  =(p/e)^  is  a  real 
quantity.  Substituting Eq .  (13)  into  Eq.  (3)  gives 


lag 


Ea  x  (er  x  jig)  /z 


(14) 


If  we  require  that  the  aperture  be  located  adequately  far  from  the  probe 
volume, 

Ea  *  er  =  0  at  the  aperture  and  it  follows  that 


Mia 


is)/z 


(15) 


Inserting  Eq.  (15)  into  Eq.  (7)  gives 
W  =  /  r2  dft  (er-er)  (E^  .  E^)/z 


(16) 


With  the  electrical  conductivity  of  the  medium  equal  to  zero,  the  intrinsic 

* 


impedance  is  real,  and  it  may  be  noted  that  Wg 


Wa&  so  W  and  ^22 


are  real  and 


W12  +  W2!  "  2  Re  U12 


W12  W21 


i2  lm  W 


12 


(17) 

(18) 
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Substituting  Eqs.  (17)  and  (18)  into  (10)  and  (11) 
W3c  =  {(Re  W12)2  +  (Im  W12)2}’i 

wdc  -  Cwn  +  w22) 

and  4>  =  tan  /  Im  \ 

\  Re  W12  / 


V  =  2  { (Re  W12)2  +  (IM  W12)2) 


(Wn  +  w22) 


(19a) 

(19b) 


(19c) 


(20) 


In  Mie  theory  or  refraction  theory,  the  formulas  which  give  the  solutions  to 
scattering  of  an  incident  plane  wave  are  written  in  terms  of  a  coordinate 
system  defined  by  the  directions  of  propagation  and  polarization  of  the  in¬ 
cident  wave  with  Z  the  direction  of  propagation  and  X  the  direction  of  polarization. 
In  matrix  form,  the  electric  field  components  of  far-field  Mie  scattering 
in  a  spherical  coordinate  system  (r,  0 > <f>)  associated  with  an  (x,y,z) 
coordinate  system  may  be  expressed  as 


"E0(0,4>y 

iEo  exp  (-ikr) 

"  -COS  4>  S2  (0)“ 

_E4>(0  40. 

*  kr 

sin  4>  Sj  (0)  ^ 

where  E0  is  the  magnitude  of  the  incident  electric  field  of  either  beam. 
S j  (0)  and  S£  (6)  are  the  amplitude  functions  of  Mie  theory. 

Eq.  (21)  may  be  written  as 


"E0(e,*r 


-  c  C  (0,*) 


(22). 


6 


with  e  =  E0  exp  (-ikr)/kr 


(23) 


C  (0  ,<J>) 


-COS  <)>  S2  (&) 

.  sin  4>  Sl  (0) 


(24) 


Now  introduce  an  alternate  coordinate  system  for  each  beam  as  shown  in 
Figure  2.  Considering  the  scattering  of  beam  1  and  2  as  separate  events 
in  their  own  spherical  coordinate  frames  (Sj,  <J>j)  and  (02 * 4>2^  8ives 


"Eqj  (Sj,  <t>i)‘ 
-E<h  (&v  ♦j). 


e  C  (1) 


(25) 


and 


*Eq2  (02>  ^2)  -  e  c  (2) 
_E{j>2  (®2*  ^2^' 

where  C  (1)  =  C  (0i,<t>i)“ 


cos  <f>^  S2  (0^) 


sin 


v2  SI 


(0l). 


and  C  (2)  -  C  (02,4>2)  - 


►a  sa  (V 


sin  <f>2  S2  (O2). 


(26) 


(27) 


(28) 


The  scattered  electric  fields  as  measured  in  the  principal  coordinate 
system  (0,<J>)  are  the  Ea  and  12^  of  Eq.  (15). 


*1  =  E!0  1  +  El* 


and  E2  =  E2q  Q_  +  E2({ri 


(29) 

(30 


i  Jj> 


7 


where  0^  and  £  are  unit  vectors. 

The  magnitude  of  these  components  as  measured  in  the  (0,<|>  )  system  may 
be  obtained  by  orthogonal  transformations  on  the  components  as  measured 
in  the  beam  coordinate  systems  (0 1,<|>)  and  (02, $2) • 


rF  *1 

he' 

=  A 

i 

w 

T> 

i— » 

CD 

H-* 

f— ‘ 

V-/ 

i 

Le6i  (e^+pJ 

(31) 


J20 


LE 


2<|> 


E02  (.QlA])' 

Le<(,2 


where 


all  a12 


-  a21  a22 


bll  b12 
Lb21  b22 


(32) 


(33) 


A  and  B  are  the  real  orthonormal  matrices  of  the  transformation  with 
a22  "  aH>  a21  ^  ”a2i 9  b22  “  bll’  and  b21  **  ~b12  which  results  in 

A+  A  =  B+B  =  I  (34) 

formulas  for  elements  of  A  and  B  were  given  by  Hong  and  Jones  and  are 
also  given  in  the  Appendix.  Sibstituting  Eqs.  (25)  and  (26)  into  (31) 
and  (32)  gives 

Ej  -  Ae  C(l)  (35) 

E2  -  Be  C(2)  (36) 


Matrix  equations  may  now  be  obtained  for  ,Eot. _E0  which  allow  the  calculation 
of  Wag  by  numerical  evaluation  of  Eq.  (16).  This  is  accomplished  by  noting 
that 


* 

Ea  •  E0 


+ 

E0  Ea 


(35) 


where  Ea,&  are  the  vectors  as  given  by  Eqs.  (29)iand  (30)  and  Ea,8  are 

+ 

the  matrices  as  given  by  Eqs.  (31)  and  (32).  Ea  is  the  complex  conjugate 
of  the  transpose  of  Ea. 


for  (a,B)  =  (1,1) 
+ 


* 

-i  •  Ei 


=  Ei  Ei 


^Aec(l))  +  ^Aec(l)^ 

=  e*e  C+(l)  A+A  C(l) 
"  e*e  C+(i)c(l) 


(36) 


since  A  A  =  I( 

For  (a, 3  )  =  (2,2) 


E2  /  |2  =  E+  E2  -  ^EeC(2))+  (pe  C(2)) 

=  e*e  C+(2)  B+B  C  (2) 


=  e*e  C  (2)  C(2) 


(37) 


since  B+B 


For  (a, 3)  =  (1,2) 
* 


E1  •  ^2  =  E2  E1  =  ^3eC(2))+(aec(1)) 

Substituting  Eqs.  (23),  (36),  (37),  and  (38)  into  (16)  gives 
WaB  =  |Ej2 


la  3 


(39) 


■1/2 


ZK 


where 


IaB  =  J* 


fn  =  C  (1)  C(l) 


(40) 

(41) 


22 


C  (2)  C(2) 


f 12  =  C+(2)  B+  A  C  (1) 


Inserting  Eq.  (39)  into  Eq.  (20)  gives 
V  =  2  ((Rel12)2  +  (Iml12)2)  ** 

(h.*  ’22) 


(42) 


1>S 


For  Mie  theory,  it  is  seen  to  be  necessary  to  evaluate  three  integrals: 


Vl2  "  /d  “  Re  f  12  "  fd  52  Fj 


Ira  *12  =  Ira  f^2  “  F2 

(I11  +  I22)  "  fdQ  (fll  +  f22}  =  fdQ  F3 

For  the  actual  evaluation  of  these  integrals,  define  a  rectangular 
coordinate  system  (x! ,  y’)  in  the  plane  of  the  aperture  with  orgin  at 
the  center  of  the  aperture  as  shown  in  Figure  3  with  (P  ',(}>*)  as  the 
corresponding  polar  coordinate  system.  The  numerical  computation  of  the 
visibility  function  over  the  range  of  geometrical  configurations  encountered 
in  practice  requires  a  numerical  integration  formula  which  may  attain 
arbitrarily  high  polynomial  accuracy  by  the  addition  of  points.  Figures 
4  and  5  show  the  location  of  points  for  a  Legendre-Gaussian  scheme. 

Figure  4  has  16  points  and  Figure  5  has  36  points.  This  numerical  integration 
scheme  over  a  circle  in  the  x'-y*  plane  may  be  expressed  as 


If  dx'  dy'  f  <x',y')  =  up2ax  E±>j  wj.j  f  (p£,*j) 


where  Pjnax  is  the  radius  of  the  circle  and  w^j  is  the  weighting  factor 
for  point  (p^,  4>^) . 


The  integrals  of  immediate  interest  are  of  the  form  /  dft  F,  ( 0  1 ,  cf) 1 )  . 


r 


1 


/  dx’  dy’  0^)2 
r2 


Fk(8 


cos0'  =  rG 


r 


fd  £1  Fk  (O'  ,<))')  = 


/dx'dy'  cos20  ’  Fj,  ( 0  '  , e{) '  ) 


1^6 


in 


-i.  *PL:  Ei,j  Wij  Cos2  0ij  F*  <9i’  *j> 


*  ^>2  Eij  Wij  COs2  6i  FK  (6i >*j> 


/dn  FK  (0' » 4> * )  -  it  tan2  Z±^  Wi-j  cos2  0'  FR(e^  *') 

i  =  number  of  circle 

j  -  number  of  point  on  a  particular  circle 

tan  0^  = 


tan  0 1  =  p 

max  max 


tan  0 
tan  0 


i 

i 

* 

max 


£i 

°niax 


0  <  pj^_  <  1 
Pmax 

The  ratio  Pi/Praax  is  specified  by  the  numerical  integration  scheme.  Note 

that  the  F ///  of  the  aperture  is  giVen  by  F ///  *  -^9—  B  } _ _  _J _ 

2  pmax  2(pmx)  2tan  6'^ 


tan  0 

i 

max  = 

1 

2f /# 

°i  = 

tan  1  / 

’fl 

tan 

eie 

\ 

-1  . 
tan  f 

.  pmax 

pi  . 

JL_ 

V 

P 

max 

2 
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/dn  Fk  (O’,  $')  -  jt _  Ey  Wy  cos2  0j  Fk  (0j,  (f'j) 

4(f///)2 

The  procedure  for  numerical  computation  of  the  visibility  function  is 
outlined  in  Figure  6. 

The  original  motivation  for  development  of  the  UTSI  Mie  theory 
visibility  code  was  to  provide  the  capability  for  critical  review  of 
the  published  work  of  other  researchers.  Such  review  is  possible  only 
if  all  details  of  the  scattering  geometry  is  available  for  a  particular 
set  of  data. 

Calculation  of  visibility  for  paraxial  forward  scatter  to  an  aperture 
with  beam  stops  was  reported  by  Roberds. ^  Roberds1  calculation  was  accomp¬ 
lished  using  scalar  diffraction  theory  and  a  numerical  integration  over 
the  aperture  excluding  the  beam  stops.  His  theoretical  and  experimental 
data  for  a  twin  beam  stop  configuration  using  water  drops  was  presented  in 
his  Figure  8.  The  UTSI  Mie  code  recently  calculated  the  visibility  for 
this  experiment  both  with  and  without  the  beam  stops  with  results  as 
shown  in  Figure  7.  As  reported  by  Robinson  and  Chu,  a  real  index  of 
refraction  results  in  an  oscillatory  phenomena  for  visibility  calculated 
by  the  Mie  theory  in  paraxial  forward  or  back  scatter.  For  the  case  under 
consideration  this  phenomenon  is  again  observed.  However,  the  UTSI  Mie 
results  using  beam  stops  are  very  close  to  Roberd's  scalar  calculation. 

Both  these  calculations  are  observed  to  be  slightly  lower  than  the  experimental 
data.  The  Mie  calculation  without  beam  stops  agrees  quite  closely  with 
the  closed  form  for  visibility  as  given  by  Farmer.9  Roberds1  Figure  10 
reported  backscatter  data  for  water  drops  with  fringe  periods  of  38.1ym 


iy6 


and  66.6  ym  with  F  number  of  collecting  aperture  not  specified.  This 
visibility  data  showed  the  oscillatory  phenomenon  that  has  led  many 
researchers  to  question  whether  the  LDV  may  be  used  for  paraxial  back- 
scatter.  A  UTSI  Mie  computation  for  fringe  period  of  50  ym  with 
f/4.0  aperture  also  shows  these  oscillations  in  Figure  8. 

However,  Figure  9  for  fringe  period  of  20  ym  shows  a  distinctly 
more  mono tonic  behavior.  This  is  also  the  case  for  Figure  10  for  fringe 
period  of  10  ym.  Figures  11,  12,  13,  and  14  show  the  result  of  adding 
absorption  (imaginary  index  of  refraction)  to  the  particles  of  Figure  9. 

The  UTSI  Mie  code  has  also  been  used  to  review  the  work  of  Bachalo^ 
for  large  particles  using  an  off-axis  location  with  aperture  centered  in 
the  x-z  plane  which  is  normal  to  the  plane  of  the  beams.  Bachalo  computed 
the  visibility  using  refraction  theory  in  a  numerical  integration  over 
the  aperture  and  also  reported  experimental  data  for  angle  3  of  30°. 

His  experimental  data  for  3  ~  30°  shows  excellent  agreement  with  his 
computer  calculation  as  does  the  Mie  calculation  of  Figure  15.  Mie  and 
refraction  calculations  have  since  been  made  at  UTSI  (figures  16-23)  which 
indicate  that  the  refraction  calculations  are  quite  good  for  angles  of  3 
between  20°  and  60°  and  a  reasonably  broad  range  of  fringe  periods. 
However,  for  angles  of  less  than  20°  and  greater  than  70  the  Mie  theory 
begins  to  deviate  somewhat  from  the  refraction  theory.  It  is  believed 
that  this  is  due  to  diffracted  light  for  the  small  angles  and  due  to 
reflected  light  for  the  large  angles. 
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APPENDIX 


The  equations  which  relate  the  various  coordinate  systems  are 

x  =  r  sin  0  cos  <|>  xl,2  =  x 

y  =  r  sin  0  cos  <|>  Yl,2  =  y  cosY  -  z  sinY 

z  -  r  cos  0  zi,2  =  ±  y  siny  +  z  cosY 

which  give 

sin  0j  2  cos  $1,2  =  sin  0  cos  ^ 

sin  01>2  sin  ^  2  “  sin  0  sin  |  cos  y  ±  sin  y 

cos  0j  2  =  ±  sin  0  sin  4>  sin  y  ±  cos  0  cos  y 

all 

=  (sin  0  cos  Y  i  cos  0  sin  $  sin  y)/sin  0.  ~ 

bu 

a12 

.  =  ±  cos  <j>  sin  y/sin  0.  „ 

b12  lt* 

a22  au  a21  =  -a12  b22  *=  bu  b21  =  b12 
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Figure  1.  Scattering  geometry.  Ki  and  K2  are 

the  propagation  vectors  of  the  incident 
Laser  beams  in  the  (y,z)  plane. 
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Figure  2.  The  beam  coordinate  systems  with 

x  =  x^  =  X2  normal  to  the  (y,z)  plane 
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Figure  3.  The  vector  ro  locates  the  center  of  the 
aperture  ancT  the  (x',y')  plane  is  normal 
to  rQ . 
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Figure  5. 


Location  of  points  for  Legendre-Gaussian 
numerical  integration  using  36  points. 
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Figure  6,  Flow  Chart  for  Calculation  of  Visibility 
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PARTICLE  DIAMETER/FRINGE  PERIOD 

FIGURE  8 


VISIBILITY  VERSUS  PARTICLE  DIAMETER/FRINGE  PERIOD 


CD 


A 

E 

3 

0 

- 

\ 

s — 

c sr 

CM 

§  I 

LU 

_ 1  ^ 

o  >- 

0 

0 

0 

- 0 - 

cz: -  crr 

*-i  o 

1—  SJ  LU 

CJ  E  X 

LU  3  Q  h- 

_J  CD  0 

__J  CM  a?  LU 

a  cn  or  *-« 

a  cd  co  ta  x 
•  .  CD 

1 

0 

0 

D 

- 0 - 

^ii»e 

L  z  -<  o 

0 

0 

0 

— 0 - 

e 

0 

r s 

V 

0 

0 

0 

u 

0 

€ 

0 

0 

—6 - — 

- 0 - 

0 

0 

0 

KTI 

tS  CD  CO  -'t  CM 

^  s  si  s  s  s 

AiniaisiA  . 


210 


0.0  0.2  0.4  0.6  0.8  1.0  1.2 

PARTICLE  DIAMETER/FRINGE  PERIOD 

FIGURE  9 


VISIBILITY  VERSUS  PARTICLE  DIAMETER/FRINGE  PERIOD 


CD 


E 

D 

© 

© 

ft 

ns 

r-i 

0 

(f) 

S 

_J  w 

o  >- 

61 

0 

(] 

c 

COLLECT IN 
.  33 

. 6328um 

80  80  -  0 
MIE  THEOF 

0 

0 

0 

a 

CbJ  CaJ 

^i«ie 

u.s:  -<  § 

© 

© 

u 

0 

( 

tj 

© 

> 

0 

0 

0 

0 

u 

0 

0 

- u - 

o 

- y 

3 

) 

o 

C3 

CD 

(O 

CM 

<3 

■ 

• 

• 

• 

• 

• 

<51 

C3 

AiniaisiA 


211 


0.2  0.4  0.6  0.8  1.0  1.2 

PARTICLE  DIAMETER/FRINGE  PERIOD 


VISIBILITY  VERSUS  PARTICLE  DIAMETER/FRINGE  PERIOD 


25 


CD 


0 

U  3 

i 

I 

P  55 

S  *o 

LU 

A?  ^  >* 

t 

! 

i 

i 

cd  or 

m  ,  o 

©~  ta  lu 

O  -ri  E  X 

+  CD  "o1- 

_l  f\I  0,1  LU 

cp  cn  cn  • 

u  co  to  ta  x 

n  •  *  ® 

) 

«4  bJ 

^0  i  n  D  0 
luqz:  ^  a 

e 

n 

U 

© 

© 

© 

c 

©  ® 

© 

rs 

( 

W - 

0 

0 

0 

u 

) 

) 

0 

cs 

00 

CD 

CM 

• 

rH 

d 

s 

<3 

K 


AiniaisiA 


iXl 


0.2  0.4  0.6  0.8  1.0  1.2 

PARTICLE  DIAMETER/FRINGE  PERIOD 


VISIBILITY  VERSUS  PARTICLE  DIAMETER/FRINGE  PERIOD 


l  1 


CO 


AiniaisiA 


213 


PARTICLE  DIAMETER/FRINGE  PERIOD 


VISIBILITY  VERSUS  PARTICLE  DIAMETER/FRINGE  PERIOD 


CO 


0 

A 

0 

A  .  .  ..3 

*-« 

0 

Z0 

LU  o 

-10 

O  (Vs  5- 

( 

zzrzP  nr 

*-<  •  o 

h-  n  h  uj 

a  -h  x 

UJ  0  H- 

J  +  0®  o 
_J  CJ  011  LU 

ommn  m 
a  m  to  s  z 

•  •  03  f-| 

cJ  ' — (  i\j  * — i  U 

a  a  h  ©( 

ll  x  ^  e 

) 

0 

0 

 ( 

\  . 

1 

( 

-A 

J 

D 

0 

) 

U 

o 

0 

o 

-   A 

A 

0 

0 

0 

—  Z>- 

) 

) 

0 

C3  CD  CO  (\J  ca 

•  •••SB 

C3  C3  C3  C3  EJ 

AiniaisiA 


214 


0.0  0.2  0.4  0.6  0.8  1.0  1.2 

PARTICLE  DIAMETER/FRINGE  PERIOD 


VISIBILITY  VERSUS  PARTICLE  DIAMETER/FRINGE  PERIOD 


28 


CD 


3 

0 

0 

r\ 

i 

i 

- - rsi 

»— < 

a 

CO 

z  CO 

UJ 

O  >* 

- e 

i 

i 

\ 

i 

| 

0 

i 

i 

t 

i 

d 

o 

ft- 

cj  . or 

»— <  «  o 

1—  G3  LU 

U  -h  H  X 

uj  3  a  p- 

_J  +  oo  0 
_J  (V  ^US 

o  00  00  •-» 

a  OO  CD  C3  X 
.  *  OO 

K 

U 

© 

0 

0 

ft 

■XSJ  ^  GT^ 

^  i  o  o  o 

\  O 

U_  X  •<  5 

U 

0 

0 

n 

0 

0 

n 

0 

0 

0 

0 

0 

n 

tr 

< 

G 

n 

- — 

© 

D 

- e — 

0 

) 

0 

(\J 


C3 


CO 

ca 


CD 

ca 


ca 


(\J 

ca 


ca 


ca 


CD 

ca 


CD 

ca 


ca 


ca 


ca 


AiniaisiA 


21^ 


PARTICLE  DIAMETER/FRINGE  PERIOD 

FIGURE  14 


VISIBILITY  VERSUS  PARTICLE  DIAMETER/FRINGE  PERIOD 


AiniaisiA 


216 


0.0  1.0  2.  0  3.  0  4.  0  5.  0  6.  0  7,  0  8.  0 

PARTICLE  DIAMETER/FRINGE  PERIOD 

FIGURE  15 


Visibility  versus  particle  diameter/fringe  period 


30 


SJ 


AiniaisiA 


217 


particle  diameter/fringe  period 


VISIBILITY  VERSUS  PARTICLE  DIAMETER/FRINGE  PERIOD 


AiniaiSIA 


2. 0  3. 0  4. 0  5. 0  6.  0  7. 0  8. 0 

PARTICLE  DIAMETER/FRINGE  PERIOD 


AiniaisiA 


„;9 


2.  B  3.0  4.0  5.0  6.  0 

PARTICLE  DIAMETER/FRINGE  PERIOD 


VISIBILITY  VERSUS  PARTICLE  DIAMETER/FRINGE  PERIOD 


AiniaisiA 


0.0  1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0 

PARTICLE  DIAMETER/FRINGE  PERIOD 


VISIBILITY  VERSUS  PARTICLE  DIAMETER/FRINGE  PERIOD 


AiniaisiA 


221 


PARTICLE  DIAMETER/FRINGE  PERIOD 

FTGIIRE  20 


VISIBILITY  VERSUS  PARTICLE  DIAMETER/FRINGE  PERIOD 


s 


Q 


F3  S 

AimaisiA 


(S3 


Q 


LlI 

CL 

LlI 


QC 


02 

LlI 

H- 

LlI 


_ I 

a 


< 

CL 


222 


FIGURE  21 


VISIBILITY  VERSUS  PARTICLE  01 AMETER/FRINGE  PERIOD 


El 


AiniaisiA 


22  J 


0.  0  1.0  2.0  3.0  4.0  5.0  6.0 

PARTICLE  DIAMETER/FRINGE  PERIOD 


VISIBILITY  VERSUS  PARTICLE  DIAMETER/FRINGE  PERIOD 


5J 


51  S3  53  S3  53 


AiniaisiA 


224 


PARTICLE  DIAMETER/FRINGE  PERIOD 


NOMENCLATURE 


A 


o 

<C> 


D 


D 

g 

D  . 
gi 

D 

mm 

D 

s 

D 

v 

f 


f  . 

i 

F 


a 


L 


L 

c 

m 


M 

N 


Acceptance  ratio 

Particle  cross  sectional  area 

Cross-sectional  area  of  spray 
-2 

e  Gaussian  beam  intensity  radius 
Time  averaged  mass  concentration 
Particle  diameter 
Geometric  mean  diameter 

Geometric  mean  diameter  in  ith  size  distribution  mode 

Mass  mean  diameter 

Sauter  mean  diameter 

Volumetric  mean  diameter 

Fraction  of  particles  in  a  mode  size 

Number  of  measurements  in  the  ith  size  increment 

Receiver  F  number 

ith  size  increment  scatter  gain 

Illuminating  beam  intensity  distribution 

First  order  Bessel  function  of  first  kind 

Depth  of  field  limit  in  visibility  approximation 

Transmittance  path  length 

Cylindrical  particle  length 

Dimensionless  probe  volume  length  in  Z  direction 
Particle  index  of  refraction 
Mass  flow  rate 

Number  of  particles  per  second  crossing  a  plane 


z2S 


NOMENCLATURE  (cont.) 


N, 


Nt 

P(D) 


Average  number  of  signal  cycles 

Total  number  of  measurements  in  a  size  distribution 
Probability  density  distribution 
Extinction  efficiency 
Calibration  wheel  radius 


T 


Transmittance 


v 


Particle  velocity 


V  Signal  visibility 


V 

o 


V 

~z 


W, 

1 


x,y,z 


Ax 


Off  axis  probe  volume 
Receiver  limited  probe  volume 
Transmitter  limited  probe  volume 
Spray  droplet  velocity  normal  to  ag 
ith  size  increment  weighting  factor 

Probe  volume  coordinates,  z  parallel  to  bisector  between  the 
beams,  y  perpendicular  to  z  and  in  plane  of  beams,  x  perpendicul 
to  yz.  Origin  at  the  intersection  of  beam  centerlines 

width  of  PSI  sample  volume  in  x  direction 


Ay  width  of  PSI  sample  volume  in  y  direction 

A z  width  of  PSI  sample  volume  in  z  direction 

Z  Axial  distance  from  spray  nozzle  outlet 

Z(0)  Z  distance  corresponding  a  transmittance  value  of  0.8  and  a 

max  monodisperse  size  distribution 

Z(l)  Z  distance  corresponding  a  transmittance  value  of  0.8  and  a 

max  single  mode  log-normal  distribution 

Z (2)  Z  distance  corresponding  a  transmittance  value  of  0.8  and  a 

max  bimodal  log-normal  distribution 


NOMENCLATURE  (cont.) 


a 

3 

6 

6 


KN 

PN(°) 

PNa) 


max. 


Angle  between  illuminating  beams 

Receiver  orientation  for  off  axis  observation 

Fringe  period 

Equivalent  fringe  period 

Material  density  (specific  gravity) 

Number  density 

Number  density  in  a  monodisperse  size  distribution 

Number  density  in  a  single  mode  log-normal  size  distribution 

Maximum  measurable  number  density  based  on  an  acceptance 
ratio  of  0.1 

Maximum  measurable  number  density  for  a  transmittance  of 


max.tran.  0.8. 


<°> 


e 

G 

G 


gi 


Vi 


V 


D  min 


Maximum  measurable  number  density  for  a  transmittance 
max.tran.  of  0.8  and  a  monodisperse  size  distribution. 

Maximum  measurable  number  density  for  a  transmittance 
max.tran.  of  0.8  and  a  single  mode  log-normal  size  distribution 

Maximum  measurable  number  density  for  a  transmittance 
max.tran.  of  0.8  and  a  bimodal  log-normal  size  distribution 

Droplet  spray  cone  half  angle 

Mean  extinction  cross  section 

Geometric  standard  deviation 

Geometric  standard  deviation  in  the  ith  size  distribution  mode 

Doppler  signal  time  period 

ith  moment  in  a  size  distribution 

Minimum  Doppler  frequency  accepted  by  the  data  acquisition  system 
Wavelength 
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4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 
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